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1. INTRODUCTION AND SUMMARY 


The use of constant-volume balloons (CVB) has proven 
a popular method of measuring atmospheric phenomena. Certain 
questions, regarding the proper application of such balloons 
and the proper interpretation <£f resulting data', have arisen. 
Thus a need has developed for determining the true capabilities 
of the CVB for aeronautical research. This study described in 
this report was designed to satisfy this need. 

The initial phase of the study involved a literature 
survey concerned with CVB. A description of this survey is 
contained in Section 2. 

Examination of the literature revealed the need for a 
more rigorous mathematical treatment of the dynamics of CVB in 
a fluctuating flow field. Accordingly a mathematical model 
was developed which describes the response of the CVB to three- 
dimensional periodic non-homogeneous flow. A description of 
this model is provided in Section 3. 

The mathematical model previously noted was incorporated 
into a digital computer program. Section 4 provides a descrip- 
tion of this program, which is called BALLOON. 

By means of BALLOON, over 84 numerical runs were con- 
ducted with the NASA/MSFC UNIVAC 1108. These runs produced 
both digital and graphical results as discussed in Section 5. 
Because of the considerable bulk involved, most of the computer 
plots are not included. They are available in a separate 
document [1] . 

Analysis of the numerical results obtained, combined 
with a study of the work of other investigations has resulted 
in certain observations and conclusions. Section 6 contains 
such observations and conclusions. Because of limitations in 
the scope of the investigation these conclusions are not con- 
sidered to be completely general with regard to CVB behavior. 
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At the same time they do provide some insight into the nature 
of the problem of properly interpreting data collected by CVB’s 

References cited and included in Section 7. In 

addition there are 16 appendices containing supplementary 
material. Of special Significance is a first-order pertur- 
bation analysis contained in Appendix H. 
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2 . LITERATURE SURVEY 


The survey of literature consisted of three separate 
activities. The first involved an inspection of existing 
balloon literature from the files of the NASA Contracting 
Officer's Representative, Dr. George H. Fichtl, representing 
approximately 600 papers and reports. Such an inspection re- 
vealed 60 pertinent documents. The second activity involved 
utilization of the computerized information retrieval system 
at Redstone Scientific Information Center. The third activity 
consisted of a personal review of journals pertaining to mete- 
orology and atmospheric physics back to 1969. Because Dr. 

Fichtl 's files appeared to adequately cover all earlier years, 
the personal review did not proceed back to any earlier journals. 
The journals review included: 

Journal of Meteorology 

Journal of the Atmospheric Sciences 

Journal of Applied Meteorology 

Quarterly Journal of the Royal Meterological Society 
Journal of Geophysical Research 
BeitrUge zUr Physik der freien Atmospfrure 
Journal of Japanese Meteorology Society 

In the process of reviewing the literature some distinc- 
tion had to be made between different types of balloons. Clearly 
not all literature concerning balloons is pertinent to the study. 
At the same time it was recognized that there exist a number 
of types of balloons which are closely related or nearly equi- 
valent to CVB's. These include tetroons, constant-altitude 
balloons, constant level balloons, transondes, horizontal- 
sounding balloons and certain types of super-pressure balloons. 

A description of each of these types is provided in Appendix A. 
The decision was made to include in the survey literature 
pertaining to all of these types of balloons so long as the 
application of the balloon was consistent with the general 
application of CVB's. 
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The results of the survey revealed four categories of 
articles or reports. The first category is concerned with the 
general theory, design, and operation of constant-volume 
balloons as discussed in subsection 2.1. The second category, 
which is described in subsection 2.2, is concerned with the use 
of constant-volume balloons to determine the mean wind velocity. 
The third category, which is of primary interest in the current 
study, is similar to the second except that, in addition to the 
mean wind velocity, some measurement of turbulent fluctuations 
is recorded. This category is presented in subsection 2.3. 
Within the third category of literature, it is appropriate to 
recognize two subcategories labeled 3-A and 3-B. In subcategory 
3-A, the measurements of turbulent phenomena are presented in 
their simplest form, involving variation with respect to spatial 
position or time. In 3-B such data may be present, but in 
addition, some spectral analysis or correlation of such data 
is provided. Clearly subcategory 3-B is of special interest 
in the current investigation. The fourth category which is 
presented in subsection 2.4, pertains to theoretical studies 
of the motion of constant-volume balloons and related objects 
immersed in a turbulent flow field, and also to the theory of 
atmospheric turbulence. 

2.1 Category 1 Literature - General Theory, Design and 

Operation of Constant-Volume Balloons 

For the current investigation, literature describing 
the general characteristics of constant-volume balloons (CVB) 
was not of primary importance. At the same time such literature 
was of some value in providing background material and in 
developing an understanding of the various uses to which the 
balloons have been put. More than seventy papers, reports, and 
articles [2-75] Were identified in this category. Some of 
these also contained information pertinent to other categories 
as discussed in other subsections. For the sake of brevity 
the discussion of Category 1 literature which follows will be 
limited to material of special significance. 
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The first CVB investigation was reported over sixty- 
five years ago by Ley [2] and concerned the use of "balanced 
pilot balloons." The objectives of the study were to obtain 
a better picture of: 

"(a) Periodic oscillations of stratum as apart from 
variations due to altitude. 

(b) Vertical currents or rising winds. 

(c) Local eddies or other phenomena." 

These objectives correspond quite closely to those of more recent 
studies. Ley in a separate article [3] also reported on an 
automatic valve for maintaining such pilot balloons at a 
constant altitude. 

The first apparent use of CVB's in conjunction with 
radio transmitters was reported by Akerman and Piccard [4] 
shortly prior to World War II. During the fifteen years , follow- 
ing the war a number of articles appeared describing various 
methods of improving balloon performance [5-17] . The initial 
works of Lally [14,16] and Angell [15] are included in this 
group as are some of the works of Hopper [ 7,13,17] and Laby 
[13,17]. 

During the ten-year period from 1961 through 1970, more 
than 40 category 1 articles were written [18-61]. One of the 
most useful of these was a summary of the state-of-the-art in 
1961 written by Angell [18]. The introduction of Mylar as a 
material for the balloon skin was probably the most significant 
development in CVB operations in this period. Both the Ghost 
balloon experiment [26,27,29,30,31,36,39,43,50,53,57,58] and 
the EOLE experiment [28,34,35,38,44,54,57] began operations 
during this period. 

Category 1 literature during the period 1971 - 1975 
consisted of thirteen articles [62-74]. Several of these 
were concerned with the EOLE experiment [62,66-68]. One of 
the most significant articles was concerned with tetrocn drag 

4 

coefficients over a range of Reynolds number from 10 to 

5 

6 x 10 [71]. Typical drag coefficient values were .73. 
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2.2 


Category 2 Literature - The Use of the CVB's to Measure 
the Horizontal Mean Wind 


The literature survey uncovered 68 articles, reports, 
and papers dealing with the use of CVB’s to measure the hori- 
zontal mean wind velocity. Of these, 23 also contained category 
1 material, including the initial work of Ley, and have already 
been noted [2,12,15,17,18,20,21,26,29,34,43,44,48,53,55,57-59, 
63,68,69,72,74]. The remaining 42 papers [75-117] cover a 
twenty-five year span (1950-1974) and describe various programs 
in which CVB's were utilized to track the mean horizontal wind. 
One of these papers [80] was not actually so much concerned 
with measuring the mean wind as with utilizing such a wind to 
send "balloon bombs" from Japan to the United States during 
World War II. A number of investigations involved flights over 
urban areas [82-84,92,101,105,109,114] while others dealt with 
the GHOST program [87,89-91,94-96,100,117] and still others 
with the EOLE experiment [98,99,102,111,112,116,117]. Balloon 
trajectories at various altitudes were analyzed including 500mb 
[85], 200mb [l!6], 50mb [86], and 50km [108]. Clustering pheno- 
mena were studied in some cases [85,97] while the relationship 
between the dispersing of CVB f s and turbulent diffusion in the 
atmosphere was investigated in others [102,104,115]. A number 
of studies were concerned with balloon trajectories in the 
planetary boundary layer [97 , 103 , 109 , lio]. Still others were 
directed toward measurements of rain clouds and thunderstorms 
[93,1X3], 

In essentially all the studies noted the basic assump- 
tion was made that the mean motion of the balloon corresponded 
to the mean motion of the atmosphere. 

2.3 Category 3 Literature - The Use of CVB's to Measure 

Turbulent Fluctuations 

Category 3 literature can be divided into two groups, 
as noted previously. All literature which contained some 
measurement of turbulent fluctuations, either vertical or 
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horizontal, was considered Sub-Category 3-A if it did not provide 
for some form of spectral analysis or correlation of the fluc- 
tuations. Sub-Category 3-B consisted of all literature which 
contained both measurement of fluctuations and spectral analysis 
or correlation of the data. 

2.3.1 Su b-Ca tegory 3-A Literature - Without Spectral Analysis of 

Correlation 

Twenty-two different studies provided basic measurements 
of turbulent fluctuations [2,18,20,25,74,75,83,92,108,110,111, 
112,115,118-126]. The majority of these studies were concerned 
with vertical fluctuation [2,18,20,25,74,83,92,108,110,115, 

120, 122-126 ] . The data were generally presented in the form 
of a plot of altitude, temperature or pressure versus time. 

The period of oscillation were generally relatively short, being 
measured in minutes. 

Ten studies provided measurement of horizontal fluctua- 
tions [75,110-112,118-121,124,1251. As with vertical fluctuation 
the data were generally presented in the form of a plot of 
displacement versus time. Two separate groups could be identi- 
fied. The first of these was concerned with large-scale fluc- 
tuation with periods of oscillation measured in hours or even 
days [111,112] while the second group dealt with small-scale 
fluctuations with periods of oscillation measured in minutes 
[110,118-121]. 

Sub-Category 3-A literature contains a wealth of . informat ion 
but time and manpower limitations precluded any attempt to 
carry out any detailed analysis of the measured data. 

2.3.2 Sub-Category 3-B Literature - With Spectral Analysis or 

Correlat ion 

The most significant group of literature bearing on the 
current study consisted of twenty-one articles which provide 
both turbulence measurements and spectral analysis or correla- 
tions of measured data [68 , 72 , 78 , 84 , 127-143 ] . It was appropriate 
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to identify two groups within Sub-Category 3-B, In the 
first of these, the authors generally provided turbulence measure 
ments, based on tracking the CVB, and some analysis of the 
resulting data [68,72,78,84,127,129,130,135,136,138-141,143]. 
These types of analyses produced Lagrangian or quasi-Lagrangian 
representations of turbulence. It is important to note that 
a wide range of frequencies and wave numbers were covered in 
these studies. A number of the studies were primarily concerned 
with low frequency oscillations (.001 to 0.1 cycles/hr) of 
the horizontal winds [68,130,139,140,143], The remaining dealt 
with intermediate to high frequency oscillations (0.1 to 100 
cycles/hour) of the vertical winds [72,78,84,127,129,135,136, 

138]. In several of these cases the natural frequency of 
oscillation was detected [72,84,129,136]. One study utilizing 
neutrally buoyant floats for measuring vertical velocity 

-2 

fluctuations in the ocean (in the spectral range from 10 to 1 
cycle/hour), also detected the natural frequency of oscillation 
[141] . 

The seven remaining Sub-Category 3-B papers [128,131-134, 
137,142] were especially pertinent to the current investigation 
because in addition to turbulence measurements and analysis 
based on the motion of CVB's they also provide corresponding 
data and analysis based on measurements taken at a stationary 
point. Thus both Lagrangian (or quasi-Lagrangian) and Eulerian 
descriptions of turbulence were available for comparison. 

The first of these studies is especially noteworthy 
because of its completeness. Gifford [128] calculated the 
vertical velocity energy spectra over the range from 3 to 200 
cycles/hour at a height of 300 feet based on measurement, by 
fixed anemometers, CVB’s and gust equipment mounted aboard an 
airplane. He demonstrated that in terms of the frequency, 
(corresponding to the spectral maxima), the data could be 
correlated by the relation 

w = + k Uj (2-1) 


8 



where 


oo = Eulerian frequency (measured by the fixed anemometer) 

I 

= (quasi) Lagrangian frequency (measured by the CVB) 
k = wave number (measured by the airplane) 
u^ = mean wind velocity^ 4 

Angell and Pack [131] , in a study of low-level CVB 
flights from Wallops Island, obtained measurements of vertical, 
longitudinal, and transverse fluctuations. Based on such measure- 
ments they calculated values of transverse velocity variance 
and transverse turbulence intensity at altitudes from 2500 - 
3000 feet. They provided a comparison of such values with 
Eulerian values (obtained at heights from 6 to 300 feet) and 
with other Lagrangian values (obtained at altitudes from 1000 
to 3000 feet). In general the values of Lagrangian variance 
and intensity obtained by Angell and Pack were less than the 
corresponding Eulerian values and were also less than the other 
Lagrangian values. Because the mean wind velocities were not 
provided and because of differences in height at which the 
various measurements were obtained, it is difficult to draw 
quantitative conclusions from these comparisons. 

Angell [132] in another study carried out an analysis 
of the trajectories of CVB's launched from Cardington at 
altitudes ranging from 1200 to 4200 feet. Based on measurement 
of vertical fluctuations, vertical velocity spectra were 
derived extending from .01 to 1.0 cycles/minute. The spectral 
peaks of these quasi-Lagrangian spectra were then compared 
with corresponding peaks from Eulerian spectra, derived from 
measurements by wind vanes attached to a barrage-balloon cable 
at heights ranging from 600 to 3500 feet. The ratio, 6, of the 
frequency for the Eulerian spectral peak to the frequency for 
the Lagrangian spectral peak was computed and tabulated. Values 
of 3 ranged from 1.1 to 8.5. An approximate correlation of 3 
with turbulence intensity was developed. 
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Kao and Bullock I; 133] performed a comparison of 
Lagrangian and Eulerian correlations and energy spectra of 
geostrophic velocities. The frequency range for the spectra 
extended from .001 to .07 cycles/hour. The. curves for the 
Eulerian auto correlation coefficients of both horizontal 
velocity components resembled their Lagrangian counterparts 
tijut displayed larger integral time scales. The Eulerian velo- 
city spectra also resembled the Lagrangian but were shifted 
toward lower frequency. The value of $ based on the ratio of 
the integral time scale was 0.53. 

In a separate paper Kao [134] computed and analyzed the 
Eulerian and Lagrangian autocorrelations and energy spectra of 
large-scale turbulent motion at the 300 mb level. The results 
were similar to those presented in the preceding study. The 
frequency range was the same and the same resemblance of the 
Eulerian autocorrelation curves and spectra t;o their Lagrangian 
counterparts was again observed. The value of 3 was found to 
be 0.33. It is important to note that the Lagrangian values 
were based on data collected by Angell [15] from CVB flights 
from Japan to the United States while the Eulerian values were 
based on data collected by Kao over Salt Lake City. The mean 
wind speeds were not equal in the two experiments. 

Angell, Pack, Hoecker, and Delver [137] performed a com- 
parison of Lagrangian and Eulerian time-scales based on CVB flights 
past a 460 meter tower in Nevada equipped with a bidirection wind 
vane. The range of frequencies extended from ^.0003 to .01 cycles/ 
sec. Based on the frequencies corresponding to the spectral peaks, 
values of 3 from 1.5 to 4.7 were obtained. A limited correla- 
tion of 3 as a function of turbulence intensity was developed. 

In a closely related study, Angell [142] calculated the 
Lagrangian and Eulerian time scales based on CVB flights past 
a 460-meter television tower near Oklahoma City. The time scales 
were based on the spectral peaks of composite Lagrangian and 
Eulerian spectra with a frequency range from .0001 to 1.0 
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cycles/second. Values of 8 ranged from 2.4 to 13. Angell 
suggested that the presence of an urban area might reduce the 
Lagrangian time scale, thus causing larger values of 8. Corre- 
lation of 8 with turbulence intensity was demonstrated to a 
limited degree. 

2 . 4 Catego ry 4 Literature - Theoretical Studies of CVB's 

and Other Bodies Immersed in a Turbulent Flow Field 

Many of articles, papers, and reports already described 
also provided some theoretical treatment regarding the behavior 
of CVB’s [2,20,22,46,58,65,72,72,97-99,102,103,119,122,124, 
126,127,136]. Some of these dealt with predicting the expansion 
of the balloon due to pressure differences and predicting the 
equilibrium altitude [58,65,72 ]. Others dealt with predicting 
the natural frequency of vertical oscillation for the CVB and/or 
an air parcel [20,58,72,127,136]. Still others were concerned 
•with the response of the CVB to oscillations in the flow 
field [22,46,97,126]. Two papers dealt with the numerical 
simulation of the dispersion of CVB [98,99], while two others 
were concerned with relating balloon dispersion to atmospheric 
diffusion [102,119]. The behavior of CVB’s in the vicinity 
of mountains was the subject of two other studies [122,124]. 

Much of the remaining Category 4 literature was concern- 
ed with the behavior of a body immersed in a fluid [144-167] . 

One of the earliest treatments of this problem was that by 
Bassett [144], who took into account transient viscous effects. 
For the case of low Reynolds numbers, corresponding to small 
particles, a sizeable number of studies have been performed 
[145-154] . These studies generally involved particle densities 
which were much greater than the fluid density and thus differed 
significantly from the CVB problem. The rigor with which the 
governing equations were derived, however, proved useful in 
establishing the governing equations for CVB motion. Two 
other papers [155,156] dealt with the behavior of bubbles in 
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fluids, when the density of the bubble is much less than the 
density of the fluid. The basic question of the natural vertical 
oscillation in a stratfied fluid has been addressed in a number 
of studies [20,58,72,127,136] as previously noted. In addition 
to these, the original work of Brunt [157], for vertical 
oscillations of an air parcel in the atmosphere, and the work 
of Larsen [158], for a neutrally buoyant sphere oscillating in 
a stratified fluid, are worthy of note. 

Clearly the most relevant studies involving immersed 
bodies were those dealing with balloon motion (159-167]. Two 
papers [161,163] were especially pertinent. The paper by Hirsch 
and Booker [161] dealt with the response of superpressure 
balloons to vertical air motions. In developing the governing 
equation for the balloon response, however, the authors 
apparently neglected apparent mass and pressure gradient 
effects as well as the Bassett terms. In addition the equations 
describing the motions of the air and the balloon, as presented 
in the paper, appear erroneous. For these reasons the resulting 
balloon trajectories are of questionable value. 

The paper by Hanna and Hoecker [163] was concerned with 
the response of constant-density balloons to sinusoidal varia- 
tions of vertical wind speeds. The equation governing the 
balloon motion was derived with more rigor than noted above, but 
several important assumptions and/or simplifications were made 
which were not clearly stated. First, the f ollow-the-f luid- 
particle total derivative was assumed identical to the follow- 
the-balloon total derivative. Second the Bassett term was 
neglected. The first assumption is valid because the authors 
considered only a periodic velocity field, which was spatially 
homogeneous. The second simplification appears acceptable for 
large Reynolds numbers. The authors presented dimensionless 
plots for calculating phase lag and amplitude response as a 
function of the properties of the atmosphere and the balloon. 
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It is important to note that the authors did not predict any 
frequency m shift of the balloon motion with respect to the air 
motion . 

In addition to literature dealing with the behavior of 
immersed bodies, certain other Category 4 literature was identi- 
fied. This included certain studies dealing with CVB trajectories 
[168-170], Also, a number of important references dealing with 
atmospheric turbulence [171-1803 . Especially useful in this 
regard were the works of Lumley and Panofsky [1753 and Slade [1763 . 
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3 . ANALYTICAL DEVELOPMENT 


In order to gain a better understanding of the behavior 
of a constant-volume balloon immersed in a flow field, a mathe- 
matical model has been developed based on the appropriate 
governing equations. The equations describing the balloon 
motion and the fluid velocity field are presented in subsection 
3.1. These equations are presented in both dimensional and 
dimensionless form. By means of dimensional analysis, as des- 
cribed in subsection 3.2 certain important dimensionless groups 
can be identified, which are useful in presenting results with 
more general application. 

As an alternative method of treating the problem, a 
first-order perturbation analysis of a perfectly responding CVB 
has been performed as discussed in subsection 3.3. 

3 . 1 Governing Equations for Balloon Motion 

The governing equations for the motion of a balloon 
submerged in a turbulent flow field consist of the equations 
for the conservation of momentum of the balloon coupled with 
the equations for conservation of momentum of the fluid. 

3.1.1 Conservation of Momentum Equations for the Balloon 

The conservation of momentum equations were originally 
derived for small Reynolds numbers (>0.1) as discussed in 
Appendix B. Based on the derivation provided in that a.ppendix, 
but without the assumption of low Reynolds numbers, the conser- 
vation of momentum equations for a balloon can be written as 
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A combination of Equations (3-1) and (3-2) yields 


dv i d ( v i ~ u ± ) 
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Y/ith these identities Equation (3-5) can be written 
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Now the drag coefficient C D can be approximated by the relation 

S = is + Sn (Ro < v ’” > < 3 - 18 > 

where 

C D4 = .5 (3-19) 
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(3-20) 

v 

A .combination of Equations (3-17), (3-18) and (3-20) yields 



Further development of the dimensionless form of the equation 
is provided in Appendix C. 
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3.1.2 


Mathematical Model for the Flow Field 


The equations governing the motion of fluid are commonly 
referred to as the Navier-Stokes equations. Although numerical 
solutions of these equations are possible, such a procedure is 
beyond the scope of the current study. Accordingly a simple 
mathematical model was developed representing a periodic three - 
dimensional flow field which conserves mass, as described in 
Appendix D. The final dimensional form of the model is: 

u i * “i 5 ii + A sin (C “ ait + e ± ) (3-22) 

where 

£ = k (x 1 + x 2 + x 3 ) (3-23) 


0 i = 0 

02 = 2tt / 3 

03 = ~ 2 n 1 3 


(3-24) 


It is important to note that the flow field model is three-dimensional 
and does contain a mean translational velocity.* To increase 
the generality of the model it can be cast in dimensionless 
form as described in Appendix E. In such form the model can 
be written: 


Ui = u 1 6^ + sin (C - (iot + 0^) 

where 

Ui = n ± /& 

5i = u x /A 


(3-25) 


(3-26)’ 

(3-27) 


* These two features of the flow field model distinguish it from 
the model of Hanna and Hoecker [163 J , who, as already noted 
in Section 2.4, employed a model which contained no mean trans- 
lational motion and was spatially homogeneous. 
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C = k ( *1 + *2 + (3-28) 
k = k D (3-29) 
X i = x/D (3-30) 
to = toD/A (3-31) 
t = At/D (3-32) 


3 . 2 Dimensional Analysis 


Inspection of Equation (3-21) reveals certain dimen- 
sionless groups. As shown in Appendix C these groups can be 
written as: 


L 
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(3-35) 

(3-36) 

(3-37) 

(3-38) 


A comparison of Equations (3-35) and (3-38) reveals 


L 6 = V 3L 3 /,r 


(3-39) 


For a CVB the balloon density is essentially equal to the 
atmospheric density. Thus, 
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C* -10) 
(3-41) 

For simplicity the balloon is assumed to have a spherical 
shape. For this case 

C m 2 (3-42) 

and 

L 4 ~l/ 4 (3-43) 

Thus the six dimensionless groups are reduced to two: Lg and 

Lg . Furthermore, the dimensionless group, Lg , can be approxi- 
mated as: 

12v 

T ~ 0 

3 ~ AD (3-44) 


L i : v 3 

and 

h 2 ~ ’ *^3 


where v q = kinematic viscosity at the equilibrium altitude. 

This parameter is clearly the reciprocal of a type of Reynolds 
number as shown in Appendix F. It is important to note that the 

o 2 2 

two dimensionless groups, 12 v /(AD) and - g (g/R - y)D /(A T ) 

3 ° 

are essentially constant for any given balloon problem and appear 
sufficient (to first-order accuracy) for characterizing the balloon 
motion. At the same time, the six original dimensionless groups 
are not truly constants and the numerical solution of Equations 
(3-21) described in Section 4, has allowed for the actual variation 
of these groups to achieve higher-order accuracy. 

Examination of Equation (3-25) indicates that the 
flow field model contains five dimensionless parameters: u- 

■*- J 

K, to , and f. Two of these, x^ and f, represent the balloon 

coordinates in space and time as obtained from the solution to the 
equation governing the balloon motion, and thus these two parameters 
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The true significance of Equation (3-45) can best be understood 
in terms of certain examples as presented in Appendix H. The 
example based on Taylor's Hypothesis is especially noteworthy. . 
For this case the relation between the two spectra becomes 


^L(TAYLOR) (fi ’ 


X 3C>) = 



(TAYLOR) x 30 )dK 


(3-48) 


According to Equation (3-48) the quasi-Lagrangian power spectrum 
is simply a spike at ft=0, and is thus independent of the shape 
of the Euler ian space-time power spectrum. It therefore appears 
according to the first-order perturbation analysis, that the 
quasi-Lagrangian power spectrum, obtained from observing the 
balloon motion, could not be used to describe the Eulerian 
spectrum, if Taylor's Hypothesis holds. 


22 



4. COMPUTER PROGRAM BALLOON 

The equations presented in Section 2 were incorporated 
into a digital computer program entitled "BALLOON". Essentially, 
the program utilizes a fourth-order Runge-Kutta technique to 
integrate the differential equations governing the balloon 
motion. The program consists of a driver routine (MAIN); 9 
subroutines (GEOMET, INTEG, ACCEL, PROPTY, BUOYNT, APARNT, DRAG, 
BASSET, and POTFLU) ; and 2 functions (DVAL and MREF) . The 
basic orginization of BALLOON is indicated in Figure 4-1. 

A copy of the source program is contained in Appendix K. 

A description of each subroutine or function is contained in 
Table 4-1. In subsection 4.1 an explanation of all input para- 
meters is provided including the input format. A description 
of the output is contained in subsection 4.2. 

4 . 1 Inputs 

All inputs to the program are read in through MAIN. 

These inputs can be divided into two segments or blocks. To 
facilitate explaining the sequence and format for each input 
item, the actual FORTRAN statements (both READ and FORMAT) 
associated with the inputs are provided as part of the discussioi 
which follows. 

4.1.1 Block #1 Inputs 

The first set of inputs consists of data setting up 
certain initial parameters. The form of the input is as follows: 

READ (5,7002) NNN 

READ (5,7001) BCDX , BCDY1 , BCDY2 , BODY 3 , FLDY1 , FLDY2 , FLDY3 

READ (5,100) G,X,K, THETA, UB,V 

READ (5,100) A ,$MEGA , SIGMA ,D , MUREF , TEMREF , T,TLIM,DELTAT , 
♦TEMPO , PO , TEMLAP , R , X30 
100 F0RMAT ( 8E10 , 3 ) 

7001 FORMAT (12A6.8X) 

7002 FORMAT (12) 
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Figure 4-1. Organization of Balloon Program 
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Table 4-1 


Summary of Subroutines and Functions in BALLOON 


Name 

Purpose 

Called by 

Calls Up 

ACCEL 

Calculates acceleration of 
balloon 

INTEG 

APPRNT 

BASSET 

BUOYNT 

DRAG 

PROPTY 

APPRNT 

Calculates forces due to 
apparent mass and pressure 
gradients in the fluid 

ACCEL 

DVAL 

POTFLU 

BASSET 

Calculates Bassett force 

ACCEL 

DVAL 

BUOYNT 

Calculates force due to 
buoyancy 

ACCEL 


DVAL 

Calculates correction to 
convection term due to 
difference in balloon and 
fluid velocities 

APPRNT 

BASSET 


DRAG 

Calculates drag 

ACCEL 


GEOMET 

Calculates time-indepen- 
dent geometric factors 

MAIN 


INTEG 

Performs numerical 
integration 

MAIN 

ACCEL 

MAIN 

Serves as driver and 
provides for input and 
output 


GEOMET 

INTEG 

MREF 

Calculates viscosity as 
a function of temperature 

PROPTY 


POTFLU 

Calculates flow field 
derivatives 

APPRNT 


PROPTY 

Calculates temperature 
and density as a function 
of altitude 

ACCEL 

MREF 
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Table 4-2 provides a definition of each input item. Notice should 
be taken than the table is arranged in the same order as the items 
are read. 

4.1.2 , Block # 2 Inputs 

The second set of data occur in the following form: 

HEAD ( 5 , 101 )XD( 2 ) , XD( 5 ) , BCDX( 11 ) , BCDX( 12 ) , T , TLIM , DELTAT , IDT , IDS 
101 FORMAT (2E10 . 3 , 10X, 2A6 , 8X, 3F5 .0,215) 


The data in Block #2 are somewhat redundant of Block #1 data. 

This redundancy results from certain modifications to the program 
to permit dimensionless results and to provide for multiple cases 
to be run in sequence. The meaning of each Block #2 input is 
provided in Table 4-3. As before, the items in the table are list- 
ed in the same order as they are read into the program. One 
set of these inputs must be provided for each case to be run 
in sequence . 

4 . 2 Outputs 

All outputs to the program occur in MAIN. These outputs 
can be divided into 4 blocks. As with the description of in- 
puts, to aid in explaining the sequence and format of each 
output item, the actual FORTRAN statements (both WRITE and 
FORMAT), associated with the outputs are provided as part of 
the discussion which follows. 

4.2.1 Block #1 Outputs 

The first output set occurs in the following form: 

WRITE ( 6 , 9004 ) A,B,C 
9004 FORMAT ( / //10X, 1P3E12 . 4 ) 

The variables A, B, C, refer to the amplitudes of the u^ , Ug 
and Ug velocity components. The variable A is read in, but 
B and C are computed. Sample values of this output is included 
in Appendix L. 
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Table 4-2 


Variable 

NNN 

BCDX 

BCDY1 

BCDY2 

BCDY3 

FLDY1 

FLDY2 

FLDY3 

G 

X 

K 

THETA 

UB 

V 

A 

OMEGA 

SIGMA 

D 

MUREF 

TEMREF 

T 

TLIM 

DELTAT 

TEMPO 

PO 

TEMLAP 

R 

X30 


Definitions of Block #1 Inputs 
Meaning 

The number of different cases to be run in sequence 

alphanumeric label for abscissa of all plots 

alphanumeric label for ordinate of first plot 
(normally x^coraponent of velocity) 

alphanumeric label for ordinate of second plot 
(normally Xg-component of velocity) 

alphanumeric label for ordinate of third plot 
(normally x^-component of velocity) 

alphanumeric label for ordinate of fourth plot 
(normally x^-coordinate of balloon) 

alphanumeric label for ordinate of fifth plot 
(normally X 2 ~coordinate of balloon) 

alphanumeric label for ordinate of sixth plot 
(normally x^-coordinate of balloon) 

gravitational acceleration vector 

position vector for initial position 

wave number vector for the fluid 

phase angles for the components of the fluid 
velocity model 

initial fluid velocity vector 
initial balloon velocity vector 
fluid wave amplitude for x-^ direction 
Eulerian frequency for the fluid 
balloon density 
balloon diameter 

reference dynamic viscosity coefficient 

reference temperature 

initial time of balloon fluid 

final time of balloon flight 

time step for numerical integration 

temperature at altitude x^g 

pressure at altitude x 3Q 

temperature lapse rate 

gas constant 

initial balloon altitude 
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Table 4-3 


Definitions of Block #2 Inputs 


Variable 

'Meaning 


XD( 2 ) 

dimensionless wave number, k 


XD(5) 

dimensionless frequency, w 


BCDX(ll) I 
BCDX( 12 ) j 

alphanumeric label to permit identifying 
run by a "test” caption 

each 

T 

initial time of balloon flight 


TLIM 

final time of balloon flight 


DELTAT 

time step for numerical integration 


IDT 

number of integration time steps between 
data generations 

output 

IDS 

factor for increasing TLIM (final time = 

IDS+TLIM) 
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4.2.2 Block #2 Outputs 


The second set of output provides a label for each case 
and provides a list of most inputs in the following form: 

WRITE ( 6 . 0915 ) BCDX( 11 ) , BCDX( 12 ) , (XD(L2 ) , L2=l . 5 ) 

WRITE(6, INPUT) 

9015 FORMAT ( 1H1 , 25X , ' TEST CASE’2A6//’ DIMENSIONLESS GROUPS’//’ N7= ’ 1PE1 
.12.4/’ N8= ' E12 . 4/ ’ N9=’E12.4/ T N10=’E12.4/’ N11='E12.4) 

NAME LI ST/ IMPUT/G , X , K , THETA , UB , V A A , OMEGA , SIGMA , D , MUREF , TEMREF , T , 
*TLIM , DELTAT , TEMPO , PO , TEMLAP , R , X30 


The definition of all of these items are provided in Tables 4.2 
and 4.3. A sample set of this output is included in Appendix L. 

4.2.3 Block #3 Outputs 


The form of this output is as follows: 


WRITE( 6 , 9002 ) 

WRITE( 6 , 9003 ) ( ( WRA( I , J) , J=1 , 11 ) , 1=1 , 50) 
WRITE (6, 9008) 

WRITE ( 6 , 9009 ) ( ( WRB( I , J ) , J=1 , 6 ) , 1=1,50) 
WRITE (6 , 9006 ) * 

WRITE( 6 , 9007 ) ( (WRC(I, J) , J=1 , 17 , 1=1 , 50 ) 


9002 FORMAT ( ' 1 TIME SCALED TIME 
1 SCALED BALL00N VELOCITIES 

9003 FORMAT ( IX, OPF8 . 2 , 1P10E11 . 3 ) 

9008 FORMAT ( ’ 1 TIME SCALED TIME 
!'/■/) 

9009 F0RMAT( IX. 0PF8 . 2 , 1PS5E11 . 3 ) 
9006 FORMAT ( ' 1 TIME SCALED TIME 

1 FAP1/FBASS 


SCALED ATMOSPHERIC WIND 
VELOCITY RATIOS,//) 


DZ 


RE 


FDRAG/FBUOY 
FAP2 ’ / / ) 


9007 FORMAT ( IX , OPF8 . 2 , 1P10E11 . 3/20X, 1P6E11 . 3 ) 


Basic data describing the computed wind and balloon motion are 
stored in the fields WRA , WRB, and WRC. The WRA field provides 
a tabulation of wind and balloon velocity as a function of time. 
The WRB field contains the balloon coordinates and the Reynolds 
number as a function of time. The WRC field contains a tabu- 
lation of the various forces acting on the balloon as a function 
of time. Sample values of block #3 data are included in 
Appendix L. 

4.2.4 Block #4 Outputs 

This output set consists of calls to the NASA/MSFC 
UNIVAC 1108 plotting routines as follows: 
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YB = 8.8 
YT = 11.2 

CALL QUIK3L ( -1 , XL , XR , YB , YT , 1H+ , BCDX , BCDY1 , - IDX , PT , PY1 ) 
CALL QUIK3L(C, XL , XR , YB , YT , IHO , BCDX , BCDY1 , -IDX , PT , PZ1 ) 
YB = -1.5 
YT = 1.5 

CALL QUIK3L ( -1 , XL , XR , YB , YT , 1H+ , BCDX , BCDY2 , - IDX , PT , PY2 ) 
CALL QUIK3L ( 0 , XL , XR , YB , YT , IHO , BCDX , BCDY2 , - IDX , PT , PZ2 ) 
CALL QUIK3L( -1 , XL , XR , YB , YT , 1H+ , BCDX , BCDY3 , -IDX, PT , PY3 ) 
CALL QUI K3L ( 0 , XL , XR , YB , YT , IHO , BCDX , BCDY3 , - IDX , PT , PZ 3 ) 
CALL QUIK3V( -1, IHO, BCDX, FLDY2,- IDX, PT,PY) 

CALL QUIK3V ( -1 , IHO , BCDX , FLDY3 , - IDX , PT , PZ ) 

The fluid velocity components are stored in the PY1, PY2 and 
PY3 fields while the balloon velocity components are stored in 
the PZ1, PZ2 , and PZ3 fields. The lateral and vertical co- 
ordinates of the balloon position are stored in the PY and PZ 
fields respectively. The preceding parameters are plotted as 
a function of time which is stored in the PT field. Examples 
of the plots generated by this block of output data are found 
in Section 5. 
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5. NUMERICAL RESULTS 


II 

l 

t 

The computer program, BALLOON, has been used to numeri- 
cally solve the equations governing the motion of the balloon in 
a periodic, three-dimensional flow field. The numerical results 
can be grouped in two categories. The first category contains the 
results of test runs to validate the numerical model incorporated 
into the program. This category is presented in subsection 5.1. 

The second category represents the results of a series of runs, which 
are designed to provide general insight into the balloon response to 
the flow field. This category is described in subsection 5.2 

Based on these numerical results certain observations con- 
cerning the behavior of constant-volume balloons have been made. 

These observations are presented in subsection 5.3 

5 . 1 Test Case Results 

Several types of test runs were carried out. The first 
were designed to provide a means of initially validating the 
model for certain simple conditions for which first-order theo- 
retical solutions exist. These runs are described in subsection 
5.1.1. The second type of test run was designed to verify the 
generality of the dimensionless form of the numerical solution. 

The results of these runs are presented in subsection 5.1.2. 

A number of other test cases were conducted to determine the 
effect of changing or omitting certain features in the program. 

The results of these tests are discussed in subsection 5.1.3. 

5.1.1 Initial Test Cases 

A series of four test cases were carried out which were 
simple enough to permit comparison with first-order theory. 

The simplifications employed in these four cases are summarized 
in Table 5-1. For each of these cases, as discussed in Appendix M, 
the period of oscillation based on first-order theory agreed 
closely with the results of the numerical solution. One 
additional test case (#5) was conducted to ensure that the 
simplification involving an isothermal atmosphere was not a 
special case. The period of oscillation for this test run, 
according to first-order theory, and the period obtained from 
the numerical solution were nearly identical. 
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Table 5-1 



Summary of 

Simplifications for 

Initial Test Runs 





Simplifications 


Test 

Run 

# 

Drag 

Force 

Apparent 

Mass 

Bassett 

Force 

Fluid 

Acceleration 

Force 

Atmosphere 

1 

zero 

zero 

zero 

zero 

isothermal 

2 

linear 

zero 

zero 

zero 

isothermal 

3 

zero 

i balloon 
mass 

zero 

zero 

isothermal 

4 

zero 

zero 

linear 

zero 

isothermal 

5 

zero 

zero 

zero 

zero 

nonisothermal 
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5.1.2, Dimensionless Test Cases 


In order to verify the dimensionless formulation of the 
governing equations and the resulting dimensionless solution, 
two test cases (#6 and #7) were conducted with the program 
BALLOON. The dimensional input values for the two cases are 
presented in Table 5-2. While the dimensional inputs are 
clearly different for the two cases, the values of the five 
dimensionless parameters developed in Section 3.3 are the same 
for both' cases. Values of the dimensionless groups are pre- 
sented in Table 5-3. 

In terms of dimensionless time, space, and velocity the 
differences in the numerical solutions for the two cases were 
less than 1%. The dimensionless plots of the vertical compon- 
ents of the balloon and wind velocities versus time for these 
cases are presented in Figure N-l of Appendix N. Thus the 
validity of the dimensionless formulation was confirmed. 

5.1.3 Other Test Cases 

Because of the fact that evaluation of the Bassett term 
represents 90% of the computation time in BALLOON, a test run 
(test case #8) was conducted, in which the Bassett term was 
set to zero with the same inputs as given for test case #6. 

The results are presented in Figure N-2 of Appendix N. As 
indicated in the figure, the absence of the Bassett term caused 
an 11% reduction in the period of the balloon motion. Although 
such a reduction is considered significant, the general shape 
of the plots of both balloon and wind vertical velocity compon- 
ents closely resemble those shown in Figure N-l. Furthermore, 
the Reynolds number associated with this test case was -^480 
which is considerably above the Stokes flow regime in which the 
Bassett term is rigorously applicable. For these reasons, 
computation of the Bassett term was omitted in subsequent runs 
unless the Reynolds number was a.1 . 
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Table * 5-2 


Dimensional Input Values for Test Cases #6 and #7 


Variable 

Case #6 

Case #7 

Units 

T 

o 

218 

300 

°K 

R 

2 . 870224xl0 6 

2 . 870224xl0 6 

cm 2 /sec 2 

y 

0 . 

0 

°K/cm 

g 

980.6 

980.6 

cm/sec 2 

% 

1 . 40646x1 0~ 4 

1.84540xl0" 4 

g/cm sec 

0) 

3.20036xl0~ 2 

2.72746xi0" 2 

rad/sec 

A 

15.44 

19.15 

cm/sec 

D 

12.0616 

17.5562. 

cm 

k 

2 . 072693xl0” 3 

1.42399xl0' 3 

cm“*^ 

5 1 

45.7555 

56.75966 

cm/sec 

a 

3.61494xl0" 4 

2 . 626856xl0~ 4 

gm/cm 3 

p o 

3.61494xl0“ 4 

2.626856xl0“ 4 

gm/cm 3 



t; 


• Table 5-3 

Values of Dimensionless Parameters for 



Test Cases #6 and #7 


Variable 

Definition 

Value 

L 3 

12^ 

AD 

2.5 x 10 -2 

L 5 

2 g (g/E-v)D 2 

3 a 2 t. 

6.25 x 10" 4 


o 


u 1 u^/A 2.96 



kD 


2.5 x 10 2 


o)D 

A 


2.5 x 10 2 
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In the setup of test cases #6 and #7, the initial balloon 
velocity was set equal to the mean. wind velocity and, in 
the flow field model, the phase angle was set to 0, with ©£ 
set to 2 t r/g, and 0g set to -2 tt / I n examining the results 
obtained, the question arose as to whether or not a change in 
the value of the phase angle would have a significant effect on 
performance of the balloon. This question is equivalent to 
the question of the influence of the initial phase relation 
between balloon and wind velocity on the periodic solution. 

To resolve this issue, three additional test cases 
(#9, 10, and 11) were conducted with phase angles as indicated 
in Table 5-4. All other input variables in test cases #9, 10 
and 11 were the same as case #6. The results obtained are 
presented in Figures N-3 through N-5 of Appendix N. Examina- 
tion of these figures reveals that for each different phase 
angle immediately following time zero there is a different 
transient buildup for the balloon motion. Subsequent to the 
transient buildup, a periodic pattern occurs. Comparison of 
Figures N-l , N-3, N-4, and N-5 indicated that the periodic 
motion of the balloon relative to the wind, following the 
transient buildup, is essentially independent of the phase 
angle . 

5 . 2 Results of Numerical Experiments 

Based on the results of test cases #1 - #11, the computer 
program BALLOON was considered acceptable for conducting a 
series of numerical experiments. To provide maximum results 
with a minimum number of computer runs the entire experiment 
was set up in dimensionless form. The experiment was carried 
out in two phases as described in subsections 5.2.1 and 5.2.2 
which follow. 

5.2.1 Phase I Numerical Experiments 

Considerable effort was devoted to devising a series 
of dimensionless runs involving variation of the primary 
dimensionless groups discussed in Section 3.2. With five inde- 
pendent parameters a complete investigation would entail a 
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Table 5-4 


Phase Angles 

Used in 

Test Cases #6 

- #11 

Test Case # 

9 1 

®2 

6 3 

6 

0 

2it /3 

-2u/ 3 

7 

0 

2w/ 3 

- 2 "/ 3 

8 

0 

2*/ 3 

- 2 t ,/ 3 

9 

7r / 2 

tt/ 2 +2tt/ 3 

tt/ 2 -2tt/ 3 

10 

TT 

TT +2tt/ 3 

TT -2tt/ 3 

11 

3tt / 2 

37r/ 2 +27r/ 3 

3tt/ 2 -27 r/ 
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very large number of runs. Accordingly in Phase I the decision 
was made to hold the variables and u^ constant with the 

values given in Table 5-5 while systematically varying k and oj 
over the ranges indicated in the same table. The individual 
runs were visualized as elements in a matrix, as shown in 
Figure 5-1. The values of k and (o for each run were assigned 
according to the position of the element in the matrix. The 
matrix was designed such that along the main diagonal 

w = k 

which satisfies Taylor’s Hypothesis. 

The dimensionless form of the input variables corres- 
pond to a large number of different physical situations. In 
this form, however, it is difficult to visualize the type of 
physical situation represented. Actually in Phase I the 
dimensionless parameters were computed from "typical" dimensional 
values of the physical parameters. The dimensional values used 
are presented in Table 5-6. It is important to note that in 
certain portions of the report dimensional results are presented 
in lieu of dimensionless values. In these cases the correspon- 
ding dimensional inputs are those presented in Table 5-6. 

5. 2. 1.1 Phase I - Original 16 Runs 

For the sixteen runs shown in Figure 5-1 a quantita- 
tive analysis of the numerical results is presented in Appendix 0. 
A qualitative view of the results of these experiments is pro- 
vided by means of Figures 5-2 through 5-6. These figures depict 
the transverse and vertical position of the balloon and the 
velocities in all three directions as functions of time for 
each test case within the matrix. In each figure, the sixteen 
plots are arranged in order of test case number so as to corres- 
pond to the elements in the test case matrix. Note that the 
scales differ from plot to plot within these figures. The 
scales were selected to provide a maximum of 10 complete 
oscillations of the forcing function. Also note that the equi- 
librium point is represented by a darkened line and that the 
balloon velocity may exceed the upper or lower limits of the 
graphs over certain intervals. In such situations portions of 
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Table 5-5 


Values of Dimensionless Parameters for Phase I 


Parameter 


Definition 


Value 


12 V o 
AD 


3.5914x10 


-4 


— g (g/ R _ 1 . 7317X10 -3 


3 AT 


u. 


u-/A 


10 


kD 


8.168x10 4 to 8 . 168x10*" 1 


03 


o)D 

A 


8 . 168xl0~ 3 to 8.168 


N, k 

8.168x10 4 

8.168X10" 3 

8.168xl0“ 2 

8 . 168xl0 _1 

8.168xl0" 3 

A 11 

A 12 

A 13 

A 14 

8.168xl0“ 2 

A 21 

A 22 

A 23 

A 24 

8.168X10 -1 

A 31 

A 32 

A 33 

A 34 

8.168 

A 41 

A 42 

A 43 

A 44 


Figure 5-1. 


Basic Matrix 


A. . for Cases in Phase 
1 3 


I 
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Table 5-6 


Typical Dimensional Values for Phase I 


Physical 

Parameter 


Value 


218° K 


R 

Y 


2.870224 x 10 6 cm 2 /sec 2 °K 


0.0 °K/cm 


S 


980.6 cm/sec 


o 


1.40646 x 10 ^ g/cm sec 


100.0 cm/sec 


D 


130.0 cm 


u. 


1000.0 cm/sec 


3.61494 x 10“ 4 g/cm 3 


3.61494 x 10“ 4 g/cm 3 


6.28 x 10 6 to 6.28 x 10 3 cm 1 * 


03 


6.28 x 10 3 to 6.28 x 10~^ sec - ^ * 


♦Such ranges of k and u> correspond to spatial wave lengths 
from 10 m to 10 km and time periods of 1 sec to 1000 sec. 
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TRANSVERSE POSITION 












Figure 5-2. x 2 -Coordinates of Balloon Position 

(Phase I) 
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VERTICAL POSITION 




TIME 




Figure 5-3. Xq-Coordinates of Balloon Position 
* (Phase I) 
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LONGITUDINAL VELOCITY 
COMPONENT 



b3 





VERTICAL VELOCITY 
COMPONENT 



1-5 




velocity plots are "chopped off". Examination of Figures 5-2 
through 5-6 reveals that all runs, which are off of the main 
diagonal, display a periodic motion which is in good agreement 
with first-order theory. As shown in Appendix 0, the period 
of oscillation according to first-order theory is given by 

~ - 2tt 

k | u x -c| (5-1) 

where 

c = (5-2) 

k 

Along the main diagonal where Taylor's Hypothesis is satisfied, 
no long-term periodic motion is observed. The balloon experiences 
one or two oscillations after ^ihich the velocity components 
appear to reach steady-state values. It is important to note, however, 
that under such steady-state conditions the balloon drifts 
laterally (perpendicular to the mean wind direction) and general- 
ly remains below its equilibrium altitude. 

5.2. 1.2 Phase I-Intermediate Runs in Cruciform Arrangement 

The marked difference between on-diagonal and off- 
diagonal runs suggested the need for runs at intermediate values 
of wave number and frequency in the vicinity of the main dia- 
gonal. Accordingly two sets of runs in cruciform arrangements , 

C 22 and , about elements A^ 2 and A 33 > as shown in Figure 5-7, 
were carried out. A quantitative description of the results of 
these runs is presented in Appendix P. The basic trend observed 
in both sets of runs was one of increasing disagreement between 
first-order theory and observation, as conditions satisfying Taylor’s 
hypothesis were approached along any branch of either cruciform. 

This trend indicates that nonlinear effects become significant 
under conditions which (nearly) satisfy Taylor's Hypothesis. 

It was also observed that the bounds of the region in which such 
nonlinear effects occurred were related to the magnitude of 
the product of 1 the dimensionless period of oscillation, as computed 
by Equation (5-1) according to first-order theory, and the dimen- 
sionless wave number. 7 > 
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Figure 5-7. Relative Position of Runs in 

Cruciform Arrangements C 22 and C 33 
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5 . 2 . 1 . 3 Phase I~Special Runs Concerning Lat eral D rift 


In all cases of balloon motion studied during Phase I, 
a lateral(x 2 -direction) drift was noted. That is, the balloon 
appeared to be subject to a mean lateral velocity. This drift 
was observed for the 16 cases comprising ^the original experiment 
matrix as shown in Figure 5-2. An example of this drift is shown in 
Figures 5-8 and 5-9. Figure 5-8 represents a case off the prin- 
cipal diagonal of the experiment matrix and shows the super- 
position of the oscillation due to the fluid oscillations and 
mean drift. Figure 5-9 represents a case on the principal 
diagonal and shows an initial or transient displacement of the 
balloon (a feature found in many of the cases examined both on 
and off the principal diagonal) followed by a final mean drift 
in the opposite direction. 

In examining long periods of balloon motion, it was 
found that the final direction of the lateral balloon motion 

•j 

depended on which side of the principle diagonal of the experi- 
ment matrix the experiment occurred. Thus, the final direction 
of drift appeared to depend in general on the sign of the 
quantity (u - c), where u is the mean wind velocity and c is the 
phase velocity. When c >u, a positive lateral drift was observed 
for the flow field under examination and when c <u, a negative 
lateral drift was observed. Examination of the flow field model 
confirmed in general that with such inequalities, the balloon 
would be subject to lateral forces consistent with the drift 
observed. However, since a negative lateral final drift was 
observed when u = c (as shown in Figure 5-9), this rule was not 
considered sufficient and further tests were necessary. 

Several runs were prepared to examine the effects of 
modifications in the flow field on the balloon drift. These 
cases consisted of (1) varying the initialization of the balloon 
motion, (2) modifying the initial phase relationships among 
the velocity perturbations, and (3) setting u = 0. These tests 
were carried out for the two matrix points, and A.^, as 

indicated in Figure 5-7. 
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TRANSVERSE POSITION, X 2 lent) 



Figure 5-8. Transverse Balloon Drift (u > c) 







For both cases the first two procedures produced results 
which generally resembled those of the original experiment . 

Neither altering the initialization by setting 

v ="5 (5-3) 


rather than 

v = u i 


(5-4) 


nor changing the phase angles of the y and z perturbations. 


0 2 < — >© 3 ( 5 - 5 ) 

quantitatively affected the direction of the drift. Setting u . — 0 
alters the original experiment such that c>u at all times. With 
this point in mind, the results, which showed the balloon drift- 
ing off in the positive x 2 ~direct ion , were consistent with the 
■earlier experiments. 


The preceding results indicated that there was no 
special bias- in the model which caused the lateral drift along 
the main diagonal where Taylor's Hypothesis for the fluid was satis- 
fied. Further study revealed that the mean balloon velocity, v, 
(instead of the mean wind velocity, u,) was the appropriate 
variable to be compared with the phase velocity, c. Thus the 
dividing line for lateral drift corresponds to 

v = c (5-6) 


representing Taylor's Hypothesis for a finite parcel. This 
relationship is satisfied slightly to the left (c>u) of the 
principal diagonal of the original matrix. Three test cases 
were performed to confirm this relationship. The first, desig- 
nated T^. , corresponded to v <c; the second, Tg , represented v '= c; 
and the third, T 3 , v > c. Figures 5-10 through 5-12 illustrate 
how the lateral velocity asymptotically reaches a constant 
value in accordance with 


>0 

(v < 

c) 


=0 

(v = 

c) 


<0 

(v > 

c) 

(5-7) 
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DIMENSIONLESS TRANSVERSE COMPONENT OF VELOCITY, V 




DIMENSIONLESS TIME, t SAl-OttH 


Figure 5-12. Transverse Balloon Velocity for Case (v > c) 
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Similarly, Figures 5-13 through 5-15 show the lateral drift 
changing direction.* 

Thus the question of lateral drift was resolved and 
the distinction between Taylor's Hypothesis for the fluid 
(u * c) and Taylor's Hypothesis for the parcel (v = c) was 
clearly established. 

5.2.2 Phase I I -Numerical Experiments 

As noted in subsection 5. 2. 1.3 the results of Phase I 
revealed the distinction between Taylor's Hypothesis for the fluid 
(THF) and Taylor's Hypothesis for a finite parcel (THP). In 
addition, there was some indication that nonlinear effects, 
became significant when the product of the dimensionless wave 
number and the dimensionless period of oscillation, according 
to first-order theory, exceeded some limiting value. For these 
reasons a second experiment matrix B^ k was designed. The rela- 
tion between the original matrix A. . and the second matrix is 

J 

depicted in Figure 5-16. In the matrix B^ k along each row the 
product of k • £ was held constant while along each column the 
ratio of oi/k remained constant. The arrangement of the matrix, 
along with the value of k and w for each element is shown in 
Figure 5-17. The values of all other parameters are the same 
as for the first matrix. 

Also included in the figure is the product of the 
dimensionless wave number and the dimensionless period of oscil- 
lation (according to first-order theory) for each element of 
the matrix. It is important to note that along each column this 
product remains constant , with THF being represented by the 
fourth column. Because the mean balloon velocity is a dependent 
variable, it is not possible to satisfy THP exactly along any 
column of the matrix. Based on analysis of the special runs of 
Phase I, however, the third column was empirically designed to 


* Notice should be taken that the direction of lateral drift when 
v 4 c depends upon the flow field model. Thus the inequalities 
presented in Equation (5-7) might be reversed in some cases. 
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DIMENSIONLESS TIME, t 


Figure 5-13. Transverse Position of Balloon for Case T 
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■mu limn mill in mu mi ■ ■ 


1 (v < c) 







TRANSVERSE DISTANCE, Xj (cm) 


400,000*. 



DIMENSIONLESS TIME, t 


SAI-OKM 


Figure 5-15. Transverse Position of Balloon for Case (v > c) 
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Figure 5-16. Relation Between Experiment Matrix A. . 

of Phase I and B k ^ of Phase II J 



2.0182xl0' 4 


2 . 0182xl(T 2 


3.7124xl0" 3 

5.4361xl0“ 2 


B 12 

B 13 

B 14 

B 15 

4 . 2529xl0 -3 

4 . 3674xl0 -3 

4 . 4924xl0” 3 

5 . 5651xl0” 3 

4 . 7466xl0“ 2 

4 . 6209xl0” 2 

4 . 4924xl0” 2 

3.6265xl0” 2 


B 21 

B 22 

B 23 

B 24 

B 25 

3.7124xl0 -2 

4 . 2529xl0 -2 

4.3674xl0" 2 

4 . 4924xl0 -2 

5 . 5651xl0 -2 

5.4361xl0 -1 

4 . 7466xl0 _1 

4 . 6209xl0 -1 

4.4924xl0 -1 

3 . 6265xl0 -1 


2.0182 


B 31 

B 32 

B 33 

B 34 

B 35 

3.7124X10 -1 

4 . 2529x10” 1 

4.3674X10" 1 

4.4924X10” 1 

5. 5651x10” - 1 

5.4361 

4 . 7466 

4.6209 

4 

4.4924 

3.6265 


Note: Beneath each matrix element there are two numbers. The upper is 

and the lower is &. 


Figure 5-17. Matrix B^ for Phase II 
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approximately satisfy this hypothesis based on the value of the 
product, xk. 

A complete set of detailed plots from Phase II is 
provided in a separate document [1] . Qualitative results are 
given graphically in Figures 5-18, 5-19 and 5-20. The period 
of observation is the same for each column but none of the other 
coordinates necessarly coincide. Figure 5-18 presents the ver- 
tical velocity of the balloon and the wind at the balloon's posi- 
tion as a function of time. These results are similar to those 
found in previous cases. The vertical position of the balloon, 
shown in Figure 5-19, is seen to correspond to the vertical 
velocities. Away from THP corresponding to column 3, the balloon 
is seen to oscillate regularly. Close to THP the balloon is 
seen to reach equilibrium at progressively lower altitudes 
to the right (along the rows), or toward the bottom (along the 
columns), of the figure. Along the bottom row, it can be seen that 
the balloon does not reach equilibrium within the time period shown, 
although the balloon velocities are approaching zero. Also in the 
last case shown along the third row the oscillation is not found. 

Figure 5-20 shows the lateral position of the balloon 
as a function of time. The results are very similar to those 
discussed in subsection 5.2.1. When the phase velocity exceeds 
the mean balloon velocity (c>v) } corresponding to column 4 and 5, 
the balloon drifts toward the positive y-direction. When c<v 
corresponding to column #1 and #2, the balloon dirfts toward the 
negative y-direction. It is important to note that along column 
#3 Taylor’s Hypothesis for the parcel is not exactly satisfied. 

In row #1 of column #3 the balloon drifts toward the positive y- 
direction in near equilibrium, while in subsequent cases the 
balloon drifts toward the negative y-direction. Thus, the 
assumption that THP depends only on the product fk is found to 
be incorrect as it probably also depends on the scale of motion 
(which may also be responsible for the differences seen in row 
#3 of Figures 5-18 and 5-19) . 
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5.3 


General Observations Concerning Balloon Response 

Examination of the numerical results from Phases I and 
II of the numerical experiments revealed a number of significant 
features of the balloon behavior. First of all, the mean motion 
of the balloon in the horizontal plane did not in general corres- 
pond to mean horizontal wind. As already noted the balloon exhib- 
ited a lateral drift except when Taylor's Hypothesis (for the par- 
cel) (THP) was satisfied. In addition, when the x^- component of 
the mean wind velocity greatly exceeded the phase velocity (u 1 >>c), 
the balloon mean x^- velocity component generally was slightly less 
than the wind velocity component (v^< u^ ) as shown in Figure 5-21. 
However, as c increased relative to u^, v-^ also increased until, 
as Taylor’s Hypothesis (for the fluid) (THF) was approached, v^ 
first equaled and then exceeded u-^ . Subsequently when THF was 
satisfied, exceeded both u^ and c. For a value of c slightly 
greater (-5%) than u^ , THP was satisfied with c becoming equal to 
v r Further increases in c resulted in a reduction in v-^ until, 
for cases where c was very large compared with u^ , v^ became 
essentially equal to u^. 


Comparison of the quasi-Lagrangian frequency of the 
balloon, with the Eulerian frequency of the wind, ui, also 

proved of interest. The quasi-Lagrangian frequency was less than 
the Eulerian frequency for all values of c greater than u 1 /2. For 
smaller values of c the quasi-Lagrangian frequency exceeded the 
Eulerian frequency. This variation is shown in Figure 5-22. As 
indicated in the figure the quasi-Lagrangian frequency is zero 
when THP is satisfied (c = v 1 ). Also the quasi-Lagrangian fre- 
quency equals the Eulerian frequency for c = u^/2. 


As indicated by the figure the numerical results 
generally agreed with linear theory except in the vicinity of 
Taylor's Hypothesis where non-linear effects became significant. 

The relation between the quasi-Lagrangian frequency and the Eulerian 
frequency was observed to be 


°L 


k c - k v.^ 


(5-8) 
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Figure 5-22. General Variation of Quasi-Lagrangian 
Frequency with Phase Velocity 
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where 


a) = |kc | (5-9) 

This relation is quite similar to that suggested by Gifford [128] 
which, as. shown in Section 2.3, was of the form 

oj = + k Uf (2-1) 

In Gifford’s case the Eulerian frequencies were all larger than 
the Lagrangian frequencies and thus the absence of the absolute 
value signs was of no significance. With this fact in mind, 
combined with the observation that 

v 2 ~ u 1 (5-10) 

it can be seen than the two relations are nearly equivalent. 

Notice should be taken that the first-order perturbation analysis 
presented in Appendix H also suggested this same type of relation. 

Gifford demonstrated [128] that a relation of the form 
of Equation (2-1) was sufficient to correlate the measured Lagran- 
gian spectral peak frequencies to the corresponding Eulerian 
spectral peak frequencies. Three investigations by Angell et.al. 
[132,137,142] also developed and compared Eulerian and quasi- 
Lagrangian spectra as discussed in Section 2.3. Unfortunately, in 
none of these studies were the values of the spatial wave number k 
determined.* In every case, however, the Lagrangian frequencies 
were less than the Eulerian and there was a tendency for the ratio^ 

8 ( = a>/fi£) to increase as the wind velocity increased. These 
observations are consistent with the form of Equation (5-8). 

In two other studies discussed in Section 2.3, Kao et.al. 
[133,134] also presented Eulerian and Lagrangian spectra for 
comparison. In these cases, which involved much higher altitudes 
and lower frequencies than the studies of Gifford and Angell, the 
Lagrangian peak spectral frequencies were larger than the Eulerian. 


* It is not sufficient to calculate k according to the relation, 
k = o)/c, because this relation is based on the assumption that 
THF holds. 


66 



, Based on Equation (5-8) and Figure 5-18, it would appear 

that the phase velocities encountered by Gifford and Angell exceeded 
u^/2, while the phase velocities involved in Kao's investigation 
were less than this value. This result is generally similar to the 
observations of Mizuno and Panofsky [180] in the atmospheric 
surface layer. However, because of the fact that both Gifford's 
and Angell’ s data involved measurement of small-scale vertical 
fluctuations at low altitudes (300-4200 feet) while Kao's data 
involved measurement of large-scale horizontal fluctuations at 
higher altitudes (18000 to 30000 feet), a certain degree of caution 
must be exercised before drawing any final conclusions concerning 
the relative magnitude of phase velocities and mean wind velocities. 

Many of the investigations discussed in Section 2.3 
involved measurement of the ratio of the Eulerian frequency to the 
quasi-Lagrangian frequency 3. It is important to note that 

3 = 

.= |kc|/0£ (5-11) 

A combination of Equations (5-8) and (5-11) yields 



c 


(5-12) 

Except in the immediate vicinity of Taylor's Hypothesis, Equation 
(5-12) can be approximated as 

3 - I 1 + 
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6. CONCLUSIONS AND RECOMMENDATIONS 

Clearly the response of a constant-volume balloon to 
atmospheric turbulence is a complex problem. The current investi- 
gation has utilized a more rigorous mathematical model than pre- 
viously employed, but even so, certain simplifications were made 
to keep the problem tractable. Representing the flow field as a 
periodic function, as opposed to a random variable, is of special 
significance. Some of the responses of the balloon to the periodic 
function may not occur in the presence of the real random process. 
Thus, caution must be exercised in reaching general conclusions. 

At the same time, the numerical results obtained are in general 
agreement with observations, and a number of points are worthy of 
note. These matters are presented in subsection 6.1. In the 
course of the study a number of questions arose for which no 
answers could be obtained due to time and funding limitations. 
Certain recommendations aimed at answering such questions are 
provided in subsection 6.2. 

6 . 1 Conclusions 

First, the distinction between Taylor's Hypothesis (for 
the fluid) and Taylor's Hypothesis (for the parcel) is of consider- 
able interest because this distinction sheds light on a number of 
characteristics of the balloon motion. The balloon does not move 
with the same mean velocity as the wind and tjius when THF is satis- 
fied the balloon does not "lock in" with the wind velocity. When 
THP is satisfied, however, the balloon does match the mean wind 
direction (but not its speed) . Under such conditions the balloon 
displays no oscillation or lateral drift. 

Except in the vicinity of Taylor's Hypothesis the quasi- 
Lagrangian frequency of the balloon could be related to the Eulerian 
frequency of the wind by means of a Doppler shift law, in accor- 
dance with first-order theory. In the vicinity of Taylor's Hypothe- 
sis first-order theory was not sufficient to relate the Lagrangian 
and Eulerian frequencies. When THP was exactly satisfied, the 
Lagrangian frequency was, of course, zero. 
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Based on available measured data, it would appear that 
const ant -volume balloons generally display some oscillatory mo- 
tion and thus THP is not being satisfied. Because THP is gene- 
rally close to THF, this suggests that the latter is also not 
being satisfied. 

An examination of the available Euler ian and Lagrangian 
turbulence spectra obtained from earlier studies indicates that 
a simple Doppler shift relation such as that suggested by first- 
order theory is possibly sufficient for correlation. Unfortun- 
ately such an approach requires values of the spatial wave number 
(based on direct measurement as opposed to calculation based on 
Taylor’s Hypothesis) and these values are not generally available. 

Although the balloon did not respond exactly to either 
the mean wind velocity or the periodic velocity fluctuations, 
the wind velocity served as a forcing function. The natural 
oscillatory frequency of the balloon was only observed during the 
initial phase following the start of a test case. 

Casting the problem in dimensionless form and providing 
a dimensionless numerical solution are also significant. Unfortun- 
ately a complete evaluation of the relative importance of the per- 
tinent dimensionless groups was not accomplished and it is not 
clear whether or not a more "universal" solution can be obtained. 

6 . 2 Recommendations 

The most important recommendat ion would be to replace 
the current periodic model for the flow field with a random model. 
The resulting balloon motion could then be subjected to Fourier 
analysis and realistic Eulerian and Lagrangian spectra could.be 
generated and compared. 

Further analysis of the existing numerical data should 
be made to provide better understanding of the relationship be- 
tween the various parameters. Such analysis might lead to a more 
precise relationship between Taylor's Hypothesis (for the fluid) 
and Taylor's Hypothesis (for the parcel). 
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The dimensionless groups associated with the dimensionless 
form of the equation should be varied over a wider range of values 
to establish their relative importance. A more general solution 
should be obtained if possible. 

Further study should be made of existing Eulerian-Lagrangian 
turbulence spectra to provide more conclusive proof of the Doppler 
shift relation between these two types of systems. Methods for 
calculating the spatial wave number as a function of altitude, 
without use of Taylor's Hypothesis, should be investigated. 

Further investigation should be made into the question of 
the difference between the mean direction and speed of the balloon 
motion and that of the wind. Conditions under which such differ- 
ences are significant in the atmosphere should be defined if 
possible . 



7 . REFERENCES CITED 

1. Tatom, F. B. and King, R. L. : "Computer Plots of the 

Motion of Constant-Volume Balloons", SAI-77- -HU, 

Science Applications, Inc., May 1976. 

2. Ley, Captain C. H.; "Report on. Balloon Experiments at 
Blackpool," Quarterly Journal of the Royal Meteorological 
Society , Vol. 37, 1911, pp. 33-58. 

3. Ley, Captain C. H. : "An Automatic Valve for Pilot Balloons' 

Quarterly Journal of the Royal Meteorological Society, 

Vol. 37, 1911, pp. 247-250. 

4. Akerman , J. D. , Piccard, J. F. : "Upper Air Study by Means 

of Balloons and the Radio Meteorograph", Journal of the 
Aeronautical Sciences , Vol. 4, pp. 332-337, 1937. 

5. Spilhaus, A. F. , Schneider, C. S., Moore, C. B.: 

"Controlled Altitude Free Balloons", Journal of Meteorology 
Vol. 5, pp . 130-137, Aug. 1948. 

6. Murray, W. D. , Schneider, C. S. and Smith, J. R.: 
"Development and Utilization of Constant-Level Balloons" 
Transactions, American Geophysical Union , Vol. 31, No. 6 
pp . 843-849 December 1950. 

7. Hopper, V. D. , Wilson, A. R. W.: "A Method for Obtaining 

Long Exposures at High Altitudes Using Rubber Balloons", 
Journal of Scientific Instruments , Vol. 30, pp . 211-212, 
June, 1953. 

8. Herz, A. J., Tennent, R. M, : "Ballooning for Research" 

Science News , Vol. 27, 1953, pp . 46-58. 

9. Masterbrook, H. D. , Anderson, A. D. : "The Transosonde" 

Weatherwise , Vol. 7, pp. 79-81, Aug. 1954. 

10. Haig, Major T. 0. : "Plastic Balloons for High Altitude 

Research", Journal of Aviation Medicine , Vol. 25, pp. 351- 
353, 1954. 

11. Anderson, A. D. , Mastenbrook, A. J. : "A New Upper Air 

Data System - The Transosonde" , Bulletin American Meteoro - 
logical Society , Vol. 37, pp. 342-350, 1956, 

12. Durst, C. S.: "Constant Level Ballooning", Journal of 

the Institute of Navigation , Vol. 10, pp. 70-77,. 1957. 

13. Laby , J. E., Lim, Y. K. , Hopper, V. D. : "Level Flights 

with Expansible Balloons" II Nerovo Aimento Series V; 

Vol. 5, pp. 249-259, 1957. 


71 


14. Haig, T. O. , Lally, V. E. : "Meteorological Sounding 

Systems" Bulletin of the American Meteorological Society 
Vol . 39, pp. 401-409, 1958. 

15. Angell, J. K. : "The Use of Transosonde Data as an Aid to 

Analysis and Forecasting During the Winter of 1958-1959" 
Journal of Geophysical Research, Vol. 64, pp. 1845-1853 
1959. 

16. Lally, V. E. : "Satellite Satellites: A Conjecture on 

Future Atmospheric-Sounding Systems" Bulletin American 
Meteorological Society , Vol. 41. pp . 429-432, 1960. 

17. Hopper, V. D. , Laby, J. E,: "Importance of Global Wind 

Study in the 30 km Region" Symposium on Antarctic 
Meteorology, 1959 , pp . 466-480, published 1960. 

18. Angell, J. K.: "Use of Constant Level Balloons in 

Meteorology", Advances in Geophysics, Vol. 8, 1961 pp . 
137-219. 

19. Young, Murray J. : "Comparison of Methods for Determining 

Probable. Impact Areas in Planning Short Range Instrumented 
Balloon Flights", Journal of Applied, Meteorology , Vol. 1 
December 1962, pp . 531-535. 

20. Angell, J. K. : "Meteorological Applications of Constant- 

Volume Balloons", The Meteorological Magazine , Vol. 92, 

No. 1087, pp. 37-48 February 1963. 

21. Giles, K. C., Angell, J. K. : "A Southern Hemisphere 

Horizontal Sounding System: A Preliminary Study", 

Bulletin of the American Meteorological Society, Vol. 44. 
pp. 687-696, 1963. 

22. Booker, D. R. , Cooper, L. W. : "Superpressure Balloons for 

Weather Research", Journal of Applied Meteorology , Vol. 4 
Feb. 1965, pp. 122-129. 

23. Calder, Nigel: "Ten Thousand Balloons", New Scientist , 

Vol. 26, pp. 494-496, 1965. 

24. Gildenberg, B. D. : "Meteorological Aspects of Constant- 

Level Balloon Operations in the Southwestern United 
States" AFCRL-66-706 , October 1966. 

25. Reynolds, A. D. , Lambertb, R. L. : "Ambient Temperature 

Measurements from Radiosondes Flown on Constant-Level 
Balloons" Journal of Applied Meteorology, Vol. 5, June 
1966, pp. 304-307. 


72 



26. Lichfield, G. W. , Frykman , R. W.: "Ghost Balloons Riding 

The Skies Will Report the World's Weather", Electronics 
Vol . 39, pp. 98-106, Nov. 1966. 

27. Lally , V. E. , Lichfield, E. W. , Solot, S. B. : "The 

Southern Hemisphere Ghost Experiment" World Meteoro - 
logical Organization Bulletin , Vol. 15, pp . 124-128, 1966 

28. Blamont , J. E. , "The EOLE Experiment", Spaceflight , 

Vol. 8, pp. 139-140, 1966. ” 

29. Lally, Vincent: "Trial Balloons in the Southern Hemis- 

phere", Science , Vol. 155, pp . 456-459, Jan. 1967. 

30. Lichfield, E. W. , Frykman, R. W. ; "Electronics for 
Around-the-World GHOST Balloon Flights", Proceedings , 
Fourth AFCRL Scientific Balloon Symposium , Jan. 1967, 
pp . 59-67. 

31. Lally, V. E.: "The Use of Satellites to Track Ghost 

Balloons", American Astronautical Society 13th Annual 
Meeting, Dallas, Texas, May 1967. 

32. Booker, D. R. , Cooper, L. W. , Hart, H. E.: "Updraft and 

Air Flux Studies Utilizing Instrumented Aircraft and 
Super-Pressure Balloon-Transponder Systems" , Weather 
Science, Inc. Norman, Okla., Contract E22-79-67 (N) 

Final Report, Nov. 1967. 

33. Franch, K. E.: "Ascent Trajectory of a Superpressure 

Balloon", Journal of Spacecraft and Rockets, Vol. 4 
(11), pp. 1557-1559, 1967. 

34. Morel, P.: "Preliminary Flights of EOLE Balloons", 

Space Research VIII - North Holland Publishing Company, 
Amsterdam, pp . 1005, 1011. 1967. 

35. Morel, P., Fourrier, J. , Freif, M. , Sitbon, P.: "Eole- 
Pacific IV Field Exercises", U. S. National Aeronautics 
and Space Administration, Technical Translation of 
"Campagne Eole-Pacif ique IV, Premiere Partie" , March, 

1968 , 34 pages . 

36. Lally, V. E.: "Comments on "The Occurrence of Icing 

on Constant Level Balloons", J ournal of Applied Meteoro - 
logy , Vol. 7, April, 1968, p. 958. 


73 



37. Mason, Earl S. : "Building and Tracking of Mylar Constant 

Density Balloons at Low Cost", Journal of Applied Meteor- 
ology , Vol . 7, June, 1968, pp . 512-513. 

38. Morel,. P. , Fourrier, J., and Sitbon, P. : "The Occurrence 

of Icing on Constant Level Balloons", Journal of Applied 
Meteorology , Vol. 7, pp. 626-634 August, 1968. 

39. *, Ginsburg, Th. : "DAS GHOST-Projekt" , Wal traumf ehrt , Sept./ 

Oct. 1968, pp. 133-140. 

40. Tefft, J. D. : "Joint NCAR-GSFC Meteorological Experiment 
Employing NIMBUS-D/ IRLS , U.S. Air Force, Cambridge Research 
Labs., Bedford, Mass., Special Reports. No. 85, Dec. 1968. 

41. Jessup, E. A.: "Application of Constant-Volume Balloons 

for Determining Air Divergence and Trajectories", U.S. 
National Severe Storms Lab., Norman, Okla. , Technical 
Circular No. 7, Dec. 1968. 

42. Alexander, H. : "Flight Analysis of a Constant Level 

Expendable-Type Balloon", Proceedings, Fifth AFCRL 
Scientific Balloon Symposium , Dec. 1968, p ; p . 19-31. 

43. Lally, V. E.: "Southern Hemisphere Ghost Flight Program; 

the First Year", Space Research, Amsterdam, Vol. 8, 1968, 
pp. 1002-1004. 

44. Morel, P. : "Preliminary Flights of EOLE Balloons", Space 

Research, Amsterdam , Vol. 8, 1968, pp . 1005-1011. 

45. Nolan, George F. : "Meteorological Considerat ions for 

Tethered and Hovering Free Balloons", Earth Observations 
from Balloons Sponsored by the American Society of Photo- 
grammetry, Feb. 1969. 

46. Booker, D. R. : "Evidence of Severe Storm Rotational 

Characteristics Obtained from Superpressure Balloons 
Trajectories", Conference on Severe Local Storms, 6th, 
Chicago, April 1-10, 1969 (Preprints of papers presented 
at the Conference, 1969, p. 32-37. 

47. Hard, H. D. , Hall, D. C., Sutherland, J. L.: "Super- 
pressure Balloon-Transponder Systems: A summary of three 

years use in thunderstorm research". Weather Service, Inc., 
Norman , Okla . , WR 903, Contract E22-65-69 (N) , June 15,1969. 

48. Allen, P.W., White, R. E.: "Method of obtaining long-range 

air trajectories at low altitudes using constant-volume • 
balloons", Upper-Air Instruments and Observations: Pro- 

ceeding of the World Meteorological Organization Technical 
Conference, Sept i 8-12 , 1969 , pp^ 285-298^ 


74 



49. Fourrier, J. : "Possibility of Measuring Humidity using 

Constant Level Balloons," Upper-Air Instruments and 
Observations: Proceedings of the World Meteorological 
Organization Technical Conference, Sept. 8-12, 1969, 
pp. 313-327. 

50. Lally, V. E., Lichfield, E. W. : "Summary of Status and 

Plans for the GHOST Balloon Project", Upper-Air Instruments 
and Observations: Proceedings of the World Meteorological 

Organization Technical Conference, Sept. 8-12, 1969, 

pp. 297-312., 

51. Anonymous: "Nimbus Balloons from Ascension Island", Facili- 

ties for Atmospheric Research, No. 10, Sept. 1969. 

52. Anonymous: "Balloons, Buoys, and Satellites", NCAR Quarter - 
ly , Boulder, Colo. , No. 25, Nov. 1969, pp . 8-11. 

53. Lally, V. E., Lichfield, E. W. : "Summary of Status and 

Plans for the GHOST Balloon Project", American Metero - 
logical Society, Bulletin, Nov. 1969, pp . 867-874. 

•54. Morel, Pierre: "Eole Constant Level Balloon Flights in the 

Troposphere", Space Research , Vol. 9, 1969, pp . 586-589. 

55. Fourier, J., Morel, P. and Sitbon, P.: "Ambient Tempera- 

ture Measurements from Constant-Level Balloons", Journal 
of Applied Meteorology , Vol. 9, pp . 154-157, February 1970. 

56. Ohta, S., Ito, T.: "Superpressure Balloon and Precise 

Radio-Sonde System for the Research of Vertical Air-Motion", 
Papers in Meteorology and Geophysics , April, 1970, pp. 45-72. 

57. Spier, J. L.: "Balloons above the Southern Hemisphere", 

Hemel en Dampkring, The Hague, Vol. 68, April, 1970, 

pp . 84-86. 

58. Lally, V. E.: "Constant-Level Balloons for Sounding 

Systems", Meteorological Monographs , Vol. 11, No. 33, 
pp. 392-396, Oct. 1970. 

59. Ballard, H. N. , Beyers, N. J., Izquierdo, M.: "A Constant- 

Altitude Balloon Experiment at 48 Kilometers", Sixth AFCRL 
Scientific Balloon Symposium , June 1970. 

60. Biamont , J. E., Heinsheimer, T. F. , Porameraen , J. P. : 
"Position Determination Methods Used* to Track Superpressure 
Balloons", Monthly Weather Review , Vol. 98, pp . 756-764, 
1970. 

61. Kulinchenko, T. M.: "Primenenie Aerostatov V Issledovanii 

Svobodnoi Atmosfery", Tsentral 1 Naia Aerologicheskaia 
Observatoriia , Trudy , 100; 1970, pp . 147-155. 


75 


62. Muller, J. : "Le Programme Eole", La Recherche Spatiale, 

Paris, Vol. 10, Aug., 1971. 

63. Pinglier, A.: "Le Programme Franco-Americain Tiros N.", 

La Recherche Spatiale , Paris, 10, August 1971. 

64,. Cherry, N. J.: "Characteristics and Performance of Three 

Low-Cost Superpressure Balloon (Tetroon) Systems " , Journal 
of Applied Meteorology , Vol. 10, Oct., 1971, pp . 982-990. 

65. Sofiyev, Ye. I.: "Methods of Calculating the Ballest 

Weight for Compensated, Nonexpandable Balloons at a Given 
Altitude", Tashkeont , Sredneaziatskiy Nauchno-Issledovatel ' 
skiy Gidrometeorologicheskiy Institut , Trudy. 53, 1971. 

66. World Meteorological Organization: "Present Status of 

Constant-Level Balloon Programmes", World Weather Watch 
Planning Report No. 33, 1971. 

f t 

67. France, Meteorologie Nationale, Paris: " EOLE, 1971-1972 ", 

Met-Mar Bulletin No. 76, July 1972, pp . 25-30. 

68. MOrel, P.: "The Eole Experiment", Proceedings of the 

World Meteorological Organization Technical Conference. 
Toyko , Japan, Oct. 2, 1972. 

69. Garger , Ye. K., Pis'menskiy, V. F.: "Kotsenke Kharak- 

teristik Ballonov Postoyannogo Ob" Yema (Tetronov)", 
Leningrad. Institut Eksperimental 1 noy Meteorologii . Trudy , 



Vol. 2, 

pp, 

. 141-151, 

1972. 


70. 

Mezeix , 

J. 

F. , Isaka, 

H, , Soulage, R. G. : 

"Etude De La 


Retention D'Eau Par Des Corps Plastiques: Application 

Amx Ballons Meteorologiques Plafonnants" , Journal de 
Recherches Atmospheriques , Vol. 8, pp.. 91-98, April/June, 
1973. 

71. Hoecker , W. H. : "Tetroon Drag Coefficients from Experi- 

mental Free-Flight Data", Journal of Applied Meteorology , 
Volume 12, pp 1062-1065, September 1973. 

72. Levanon, N. et.al.: "On the Behavior of Superpressure 

Balloons at 150 mb", Journal of Applied Meteorology , 

Volume 13, pp., 494-504, June 1974. 

73. Carten, A. S., Jr.: "An Investigation of the Applicability 

of High Altitude, Lighter-Than-Air (LTA) Vehicles to the 
Tactical Communications Relay Problem", Report #AFCRL-TR- 
74-0399 , August, 1974. 


76 


74. Cadet, D. , Ovarlez, H., Ovarlez, J.: "Superpressure 

Balloon. Flights in the Tropical Boundary Layer", Journal 

of Applied Meteorology , Vol. 14, Dec. 1975, pp . 1478, 1484. 

75. Durst, C. S. and Gilbert, G. H.: "Constant-Height Balloons- 

Calculation of Geostrophic Departures", Quarterly Journal 

of Royal Meteorological Society , Vol. 76, pp . 75-86, 195. 

76. Murray, W. D. , Schneider, C. S., and Smith, J. R. : 
"Development and Utilization of Constant-Level Balloons", 
Transactions, American Geophysical Union , Vol. 31. No. 6 
pp. 843-849, December 1950. 

77. Berry, F. A.: "Constant-Level-Balloon and Meteorological 

Trajectory Studies at Project AR0WA" , Transactions American 
Geophysical Union, Vol. 35, Section of Meteorology, Paper 
#4, May 1954, p.368. 

78. Ohnsorg, F. , Mantis, H. T. : "The Structure of Upper Winds 

as Determined by Constant-Altitude Balloons", Transactions 
American Geophysical Union , Vol. 35, Section of Meteorology, 
Paper #5, May, 1954, p. 368. 

79. Moore, C. B., Smith, J. R., and Gaalswyk, A.: "On the 

Use of Constant-Level Balloons to Measure Horizontal 
Motions in the Atmosphere" Journal of Meteorology , Vol. 11, 
No. 3, pp. 167-172, June 1954. 

80. Arakawa, H.: "Basic Principles of the Balloon Bomb", 

Papers in Meteorology and Geophysics, Vol. 6, pp . 239-243, 
1956. 

81. Neiburger, M. , Angell, J. K. : "Meteorological Applications 

of Constant-Pressure Balloon Trajectories", Journal of 
Meteorology , Vol. 13, pp . 166-194, April, 1956. 

82. Holzworth, G. C., Kauper , E. K. and Smith, T. B.: "Some 

Observed Low-Level Air Trajectories over Los Angeles, 
California", Monthly Weather Review , pp. 387-392, August 
1963. 

83. Pack, D. H., Angell, J. K. : "A Preliminary Study of Air 

Trajectories in the Los Angeles Basin as Derived from 
Tetroon Flights", Monthly Weather > Review , Vol. 9, pp . 583- 
604, Dec. 1963. 

84. Angell, J. K., Pack, D. H. , Holzworth, G. C., Dickson, C.R.: 
"Tetroon Trajectories in an Urban Atmosphere", Journal of 
Applied Meteorology , Vol. 5, pp . 565-572, 1966 . 


77 



85. Angell, J, K. , Hass, W. A.: "Effect of Clustering Upon 

a 500-mb Horizontal Sounding System", Monthly Weather 
Review , Vol. 94, pp . 151-165, 1966. 

86. Lund, I. A.: "Fifty Millibar Balloon Trajectories" 

Proceedings, Fourth AFCRL Scientific Balloon Symposium , 

Jan. 1967, pp . 215-232. 

87. Solot, S. B. : "Some Meteorological Results of the Southern 

Hemisphere GHOST Experiment", Proceedings, Fourth AFCRL 
Scientific Balloon Symposium , Jan. 1967, p. 297. 

88. Allen, P. W. , Jessup, E. A., White, R. E.: "Long Range 

Trajectories”, Proceedings of the USAEC Meteorological 
Information meeting, Chalk River Nuclear Labs, Sept 11-14, 
1967, pp. 176-190. 

89. Solot, S. B. : Ghost Balloon Data, Vol. II: Balloon 28X, 

Jan. 1968, NCAR-TN-34. 

90. Solot , S. B. : Ghost Balloon Data, Vol. Ill: Balloon 34J, 

Jan. 1968, NCAR-TN-34. 

91. Solot, S.- B. : Ghost Balloon Data, Vol X: Balloon Data 79R, 

Jan. 1968, NCAR-TN-34. 

92. Ahrens, D. , Miller, A.: Air Trajectories in the San 

Francisco Bay Area , Report No. 5, San Jose State College 
Meteorology Department, Sept. 1968.- 

93. Hall, D. R. , Hall, D. C., Hart, H. E.: "Balloon Trans- 

ponsder Studies of Severe Storm Systems", Weather Science, 
Inc., Norman, Okla. , Final Report, Oct. 1968, 34 pages. 

94. Solot, S. B. : "Ghost Balloons are Solving the Riddle of 

The Southern Hemisphere Circulation", Proceedings, Fifth 
AFCRL Scientific Balloon Symposium , Dec. 1968, p. 223. 

95. Lally, V. E. , Carlson, N. , Frykman , R. W., Lichfield, E. W. , 

Rickel, A. B. , Verstraete, M. : "Superpressure Balloon 

Flights from Christchurch, New Zealand, August 1967 - June 
1968", National Center for Atmospheric Research, Boulder , 
Colo., NCAR Technical Notes , Feb. 1969, 74 pages. 

96. Lally, V., Schumann, A., Reed, R.: "Superpressure Balloon 

Flights in the Tropical Stratosphere", Science, Vol. 166 
(3906), pp. 738-739, Nov. 1969. 


78 



97. Hoecker, W. H., Angell, J. K.: "Effect of a Sudden Change 

in Terrain Height on the Three-Dimensional Low-Level Air 
Flow, as Estimated from Tetroon Flights”, Monthly Weather 
Review , Dec. 1969, pp. 845-849. 

98. Mesinger, F. , Mintz, Y. : "Numerical Simulation of the 

1970-71 Eole Experiment", L.A. Dept, of Meteorology , 
Numerical Simulation of Weather and Climate Technical 
Report No. 4 , Jan, 1970, 51 pages. 

99. Mesinger, F. , Mintz, Y.: "Numerical Simulation of the 

Clustering of Constant-Volume Balloons in the Global 
Domain" , L.A . De p t, of Meteorology, Numerical Simulation 
of Weather and Climate Technical Report No. 5 , Feb. 1970 
35 pages . 

100. Lally, V. E., Rickel, A. B. , Verstraete, M. : Superpressure 

Balloon Flights from Christchurch, New Zealand, July 1968 - 
December 1969", National Center for Atmospheric Research , 
Boulder, Colorado, NCAR Technical Notes TN-48 , May 1970, 

104 pages. ~~ ~~~ 

101. Pack, D. H. : "Tetroon Flights in Los Angeles, Calif., 1969" 

U. S. Air Resources Lab, Silver Spring, Md. , Technical 
Memorandum C-RLTM-ARL 19, June 1970. ~ 

102. Morel, Pierre: "Large Scale Dispersion of Constant-Level 

Balloons in the Southern General Circulation", Annales de 
Geophysique , Vol . 26 (4), pp . 815-828, 1970. 

103. Angell, J. K.: "Helical Circulations in the Planetary 

Boundary Layer", Journal of the Atmospheric Sciences , 

Vol. 28, Jan. 1971 pp . 135-138. 

104. Angell, J. K., Allen, P. W. , and Jessup, E. A.: 

"Mesoscale Relative Diffusion Estimates from Tetroon 
Flights", Journal of Applied Meteorology , Vol. 10, pp . 43-46 
February 1971 . 

105. Angell, J. K. , Pack, D. H. , Dickson, C. R. , Hoecker, W. H. : 
"Urban Influence on Nighttime Airflow Estimated from 
Tetroon Flights", Journal of Applied Meteorology, April, 
1971, pp. 194, 204, Vol. 10. 

106. Angell, J. K. , Pack, D. H. , Machta, L. , Proulx, R., 

Dickson, C. R. : "Mesoscale Trajectories, Helicopter 

Oxidant Measurements and Surface Wind Fields", Conference 
oh Air Pollution Meteorology, April. 5-19, 1971, Preprints 
of Papers , pp . 119-120. 


79 


107. Blamont, J. E. , Heinsheimer, T. F. , Pommereau, J. P. : 

"An Experimental Demonstration of the Isentropic Nature 
of the Motion of Stratospheric Air Masses", Journal of 
the Atmospheric Sciences, Boston, Vol. 28, No. 6, Sept. 
1971. pp. 1015-1020. 

108. Ballard, H. N. , Beyers, N. J., Miers, B. T . , Izquierdo , M. , 

Whitacre, J. : "Atmospheric Tidal Measurements at 50km 

from a Constant-Altitude Balloon", U. S. Army Electronics 
Command Research and Development Technical Report EC0M- 
5417 , December 1971. 

109. Angell, J. K. , Pack, D. H. , Machta, L. , Dickson, C. R. , 

Hoecker, W. H. : "Three-Dimensional Air Trajectories 

Determined from Tetroon Flights in the Planetary Boundary 
Layer of the Los Angeles Basin", Journal of Applied 
Meteorology , April 1972, Vol. 11, pp . 451-471. 

110. Angell, J. K. : "Comparison of Circulations in Transverse 

and Longitudinal Planes in an Unstable Planetary Boundary 
Layer", Journal of the Atmospheric Sciences , Vol. 29, 

pp. 1252-1261, Oct. 1972. 

111. Webster, P. J., Curtin, D. G. : "Interpretat ions of the 

E0LE Experiment I. Temperal Variation of Eulerian Quanti- 
ties", Journal of the Atmospheric Sciences, Vol. 31., 

pp. 1860-1875, Oct. 1974. 

112. Morel, P., Bandeen , W. : "The EOLE Experiment: Early 

Results and Current Objectives", American Meteorological 
Society Bulletin , Vol. 54, April 1973, pp. 298-306. 

113. Levine, J., Garstang, M. LaSear, N. E.: "A Measurement 

of the Velocity Field of a Cumulus Cloud" , Journal of 
Applied Meteorology , Volume 12, pp. 841-846, August 1973. 

114. Angell, J. K. , Hoecker, W. H. , Dickson, C. R. , and Pack, 

D. H. : "Urban Influence on a Strong Daytime Air Flow as 

Determined from Tetroon Flights" Journal of Applied 
Meteorology , Volume 12, pp . 924-936, September 1973. 

115. Leahey, D. M., Hicklin, H. S.: "Tetroon Studies of 

Diffusion Potential in the Airshed Surrounding the 
Crowsnest Pass Area", Atmosphere, (Toronto, Canada), 

Vol. 11, pp. 77-87, 1973. 

116. Morel, P., Desbois , M. : "Mean 200-mb Circulation in the 

Southern Hemisphere Deduced from EOLE Balloon Flights". 
Journal of the Atmospheric Sciences , Vol. 31, pp . 394-407 
March, 1974. 


80 



117. Angell, J. K.: "Mean Meridional Winds Estimated from 

Constant Level Balloon Flights in Southern Hemisphere 
Temperate Latitudes”, Royal Meteorological Society , 
Bracknell, Eng., Quarterly Journal , Vol. 100, pp. 212-220, 
April 1974. 

118. Giles, K. C. , Peterson, R. E.: "A Study of Wind Velocities 

Measured from Constant-Pressure Balloons”, Journal of 
Meteorology , Vol. 14, pp. 569-572, Dec. 1957. 

119. Wilkins, Eugene, M. : "observations on the Separation of 

Pairs of Neutral Balloons and Applications to Atmospheric 
Diffusion Theory", Journal of Meteorology , Vol. 15, 

pp. 324-327, June 1958. 

120. Laby, J. E. , and Sparrow, J. G. : "Wind Studies Using Level 

Balloons" Journal of Applied Meteorology , Vol. 4, 

pp. 585-589, October 1965. 

121. Angell, J. K. , Pack, D. H.: "Mesoscale Diffusion Derived 

from Tetroon Flights", Proceedings of the USAEC Meteor- 
logical Information Meeting. Chalk River Nuclear Labs, 

Sept. 11-14, 1967, pp . 241-251. 

122. Reynolds, R. D. , Lamberth, R. L.: "Investigation of a 

Complex Mountain Wave Situation", Journal of Applied 
Meteorology , Vol. 7, June 1968, pp . 353-358. 

123. Dickson, C. R. and Angell, J. K.: "Eddy Velocities in the 

Planetary Boundary Layer as Obtained from Tetroon Flights 
at Idaho Falls", Journal of Applied Meteorology , Vol. 7 
Dec. 1968, pp . 986-993. 

124. Vergenier, J., Lilly, D. K. : "The Dynamic Structure of 

Lee Wave Flow as Obtained from Balloon and Airplane 
Observations" Monthly Weather Review , Vol. 98, March 1970, 
pp . 220-232. 

125. Lile, R. C.: "Tetroon-Derived Eddy Viscosity in the Region 

of the Elevated West Coast Temperature Inversion", 

Scientific Report No. 1 , July, 1970. 

126. Reynolds, R. D. "Superpressure Balloons as Isentropic/ 
Isopynic Tracers" Journal of Applied Meteorology , Volume 12, 
pp. 369-373, March 1973. 

127. Emmons, G. , Haurwitz, B. , Spilhaus, A. F. : "Oscillations 

in the Stratosphere & High Troposphere" , Bulletin American 
Meteorological Society, Vol. 31, No. 4, April 1950, 

pp. 135-138. 


81 


128. Gifford, Frank, Jr. ; ’"A Simultaneous Lagrangian-Eulerian 

Turbulence Experiment" Monthly Weather Review , Vol. 83, 

No. 12, December 1955, pp. 293-301. 

129. Angell, J.K., Pack, D. H. : "Analysis of some Preliminary 

Low-Level Constant Level Balloon (Tetroon) Flights", 
Monthly Weather Review , Vol. 88, No. 7, pp. 235-248, 

July, 1960. 

130. Angell, J. K. "An Analysis of Operational 300 MB 
Transonde Flights from Japan in 1957-58", Journal of 
Meteorology , Vol. 17, pp. 20-35, 1960 

131. Angell, J. K. , Pack, D. H.: "Analysis of Low-Level Cons- 

tant Volume Balloon (Tetroon) Flights from Wallops Island" 
Journal of the Atmospheric Sciences, Vol. 19, pp. 87-98 
1962. 

132. Angell, J. K. : "Measurements of Lagrangian and Eulerian 

Properties of Turbulence at a Height of 2500 Feet" 
Quarterly Journal of the Royal Meteorological Society, 

90 (383) pp. 57-71, 1964. 

133. Kao, S. K. and Bullock, W. S.: "Lagrangian and Eulerian 

Correlations and Energy Spectra of Geostrophic Velocities" 
Quarterly Journal of the Royal Meteorological Society, 

Vol. 90, No. 384, pp . 166-174, 1964. 

134. Kao, S. K. : "Some Aspects of the Large-Scale Turbulence 

and Diffusion in the Atmosphere" Quarterly Journal of the 
Royal Meteorological Society , Vol. 91, No. 387 pp. 10-17 
1965. 

135. Kao, S. K. : "Some Characterist ics of Relative Particle 

Dispersion in the Atmosphere's Boundary Layer", Pro ceedings 
of the USAEC Meteorological Information Meeting, Chalk 
River Nuclear Labs, Sept . 11-14 , 1967, pp. 437-452. 

136. Angell, J. K. , Pack D. H. , Delver, N. : "Brant-Vaisala 

Oscillations in the Planetary Boundary Layer", Journal of 
the Atmospheric Sciences , Vol. 26, Nov. 1969, pp . 1245-1252. 

137. Angell, J. K. , Pack, D. H. , Hoecker , W. H. , Delver, N. : 
"Lagrangian-Eulerian Time-Scale Ratios Estimated from 
Constant-Volume Balloon Flights Past a Tall Tower", 

Royal Meteorological Society, Quarterly Journal , Jan. 1971 
pp. 87-92. 

138. Garger, Ye. K. et.al.:, "Razraketka i Naturnyye Ispytaniya 
Ballonov Postoyannogo Ob" Yema (Tetronov)", Leningrad , 
Institut Eksperimental ' noy Meteorologii , Trndy, Vol. 2, 
pp. 125-140. 


82 



139. ' Morel, P., Necco, G. : "Scale Dependence of the 200-mb 

Divergence Inferred from EOLE Data", Journal of the 
Atmospheric Sciences , Vol. 30, July 1973, pp. 909-921. 

» ; t 

140. Desbois, M. : "Tranferts Macroturbulents a 200 mb d apris 

les donnees de 1' experience EOLE dand 1' hemisphere 

Sud" , Annales de Geophysique , Vol. 29, pp . 252-537, 

Oct /Dec, 1973. 

141. Gascard, J. C. "Vertical Motions in a Region of Deep-Water 
Formation", Deep-Sea Research . Vol. 20, pp . 1011-1027, 

1973. 

142. Angell, J. K.: "Lagrangian-^Eulerian Time-Scale Relation- 

ship Estimated from Constant Volume Balloon Flights Past 
a Tall Tower" , Turbulent Diffusion in Environmental 
Pollution, Advances in Geophysics , Vol. 18B, pp . 419-431 
Academic Press, 1974. 

143. Desbois, M.: "Ondes St at ionnaires et progressives de la 

Circulation Atmospherique a 200 mb dans 1' Hemisphere 
Sud e'tudiees a' 1 ’ aide de 1 'experience EOLE", Annuals 
of Geophysical Transactions , Vol. 30, 1974, pp . 459-471. 

144. Basset, A. B.: A Treatise on Hydrodynamics , Vol. 2, 

Chapter 5, Deighton, Bell and Co., Cambridge, 1888 

145. Tchen, C. M. : "Mean Value and Correlation Problems Con- 

nected with the Motion of Small Particles Suspended in a 
Turbulent Fluid", Ph.D. Dissertation Delft, Martinus 
Nijhoff, the Hague, 1947. 

146. Corrsin, S., and Tumler, J.: "On the Equation of Motion 

for a Particle in Turbulent Fluid", Applied Scientific 
Research , Section A, Vol. 6, pp . 114-116, 1953. 

147. L'umley, J. L. : "Some Problems Connected with the Motion 

of Small Particles in Turbulent Fluid, Ph.D. Dissertation, 
John Hopkins University, 1957. 

148. Hinze, J. 0.: Turbulence , McGraw-Hill Book Co., N.Y. 1959. 

149. Choo, B. T. : "Turbulent Transport Behavior of Small 

Particles in Dilute Suspension" Osterichisches Ingenieur - 
Archi v . Sonderabdruck aus Bd., XVIII, Heft 1 and 2, 1964. 

0 

150. Hjelmfelt, A. T. Jr. and Mockros , L. F. : "Motion of Dis- 

crete Particles in a Turbulent Fluid", Applied Scientific 
Research , Vol. 16, pp . 149-161, 1966. 


83 


151. Shirazi, M. A..: "On the Motion of Small Particles in a 

Turbulent Fluid Field", Ph.D. Dissertation, University of 
Illinois, 1967. 

152. Gorbis, Z. R. and Spokinyi, F. E.: "Pulsating Motions of 

a Solid Particle in Turbulent Flow", Mekhanika Zhidkosti i 
Gaza , Vol. 4, No. 1, pp. 137-141, 1969. 

153. Hill, M. K. and Zukoski , E. E., "Behavior of Spherical 
Particles at Low Reynolds Numbers in a Fluctuating Trans- 
lational Flow - Preliminary Experiments", Aerospace 
Research Laboratories 72-0017 , 1972. 

154. Morsi, S. A. and Alexander, A. J. : "An Investigation of 

Particle Trajectories in Two-Phase Flow Systems", Journal 
of Fluid Mechanics , Vol. 55, Part 2, pp . 193-208, 1972 

155. Contaxes, N. A. and Tatom, F. B. : "The Motion of a Void 

in a Highly Viscous Rotating, Accelerating Flow", 
Nortro n ics-Huntsville Internal Note , No. 76, 1968 

156. Calton, I. and Schwartz, .. H. : "Motion of Bubbles in a 

Rotating Container", Journal of Spacecraft, Vol. 9, No. 6, 
pp. 468-471, 1972. 

157. Brunt, D. : "The Period of Simple Vertical Oscillations 

in the Atmosphere" , Quarterly Journal Meteorological 
Society , Vol. 53, pp . 30-31, 1927. 

158. Larsen, L. H. : "Oscillations of a Neutrally Buoyant 

Sphere in a Stratified Fluid", Deep-Sea Research , Vol. 16, 
pp. 587-603, April 1969. 

159. Kao, S. K. : "An Analysis of the Trajectory of a Constant- 

Pressure Balloon wn a Long Atmospheric Wave in Comparison 
with the Nonlinear Theory", Journal of .Meteorology , Vol. 18, 
pp. 251-256, 1961. 

160. Reed, W. H. Ill,: "Dynamic Response of Rising and Falling 

Balloon Wind Sensors with Application to Estimates of 
Wind Loads on Launch Vehicles", NASA Technical Note 
D-1821 , October 1963. 

161. Hirsch, J. H. , Booker, D. R. : "Response of Superpressure 

Balloons to Vertical Air Motions", Journal of Applied 
Meteorology / Vol . 5, pp . 226-229, April, 1966 

162. Scoggins, J. R. : "Spherical Behavior and the Measurement 

of Wind Profiles" NASA TN D-3994 , June, 1967. 


84 



163. Hanna, S. R. and Hoecker, W. H.: "The Response of Constant- 

Density Balloons to Sinusoidal Variations of Vertical Wind 
Speeds”, Journal of Applied Meteorology , Vol. 10, 

pp . 601-604 , June 1971 . 

164. Grant, F. C.: "Proposed Technique for Launching Instru- 

mented Balloons into Tornadoes", NASA- TN D-6503 , Langley 
Research Center, Nov. 1971.. 

165. Fichtl, G. H. "Aerodynamic Properties of Spherical Balloon 
Wind Sensors", Journal of Applied Meteorology, Vol. 11, 

No. 3, pp. 472-481, April 1972. 

166. Fichtl, Q. H. : "Spherical Balloon Response to Three- 

Dimensional Time-Dependent Flows", NASA TN D-6829, July, 
1972. ’ ” ‘ 

167. Fichtl, G. H.: "Behavior of Spherical Balloons in Wind 

Shear Layers", Journal of Geophysical Research, Vol. 77, 

No. 21, pp. 3931-3935, July 20, 1972. 

168. Ohnsorg, F.: "Error Lstimates for Constant-Level Balloon 

Tracking" Journal of Meteorology, Vol. 14, pp . 81-83 
February 1937 . 

169. Mesinger, Fedor: "Behavior of a Very Large Number of 

Constant-Volume Trajectories", Journal of the Atmospheric 
Sciences , Vol. 22, pp . 479-492, 1965. 

170. Nolan, G. F. (Editor A. O. Korn): "The Constancy of Winds 

in the Lower Stratosphere and Constant-Level Balloon 
Planning", Proceedings, AFCRL Scientific Balloon Workshop 
1965 , (AFCRL-66-309), May, 1966 p.9. 

171. Etkin, B.: Dynamics of Flight , John Wiley & Sons, Inc. 1959. 

172. Wandel, C. F. and Kofoed, Hansen, 0: "On the Eulerian- 

Lagrantian Transform in the Statistical Theory of 
Turbulence" Journal of Geophysical Research, Vol. 67, No. 8 
pp. 3089-3092, July 1962. 

173. Panofsky, H. A.: "Scale Analysis of Atmospheric Turbulence 

at 2m" , Quarterly Journal of the Royal Meteorological 
Society , Vol. 88, pp . 57-69, 1962. 

174. Corrsin, S.: "Estimate of the Relation Between Eulerian 

and Lagrangian Scales in Large Reynolds Number Turbulence" 
Journal of Atmospheric Sciences , Vol. 20 pp . 115-119, 1963. 

175. Lumley, J. L. and Panofsky, H. A.: The Structure of 

Atmospheric Turbulence , Interscience Publishers, 1964. 


85 


176. Slade, D. H. (Ed.): Meteorology and Atomic Energy 1968 , 

U. S. Atomic Energy Commission, Office of Information 
Services, July 1968. 

177. Priestley, M. B. : "Time Dependent Spectral Analysis and 

its Application in Prediction and Control", Journal of 
Sound and Vibration , May 1971, Vol. 17 (4), pp. 517-534. 

178. Panchev, S.: Random Functions and Turbulence , Pergamon 

Press, 1971. 

179. Etkin, B. : "Flight in a Turbulent Atmosphere", Dynamics 

of Atmospheric Flight , John Wiley & Sons, Inc., 1972, 
pp. 529-563. 

180. Mizuno, T. and Panofsky, H. A.: "The Validity of Taylor’s 

Hypothesis in the Atmospheric Surface Layer" Boundary- 
Layer Meteorology, pp . 375-380, D. Reidel Publishing Co. 
1975. 

181. Hess, S. L. : Introduction to Theoretical Meteorology 

Holt, Rinehart and Winston, New York, 1959. 


86 



APPENDIX A 


Balloon Terminology 


During the past 65 years a number of different terms have 
been used to describe that class of balloons which are designed for 
horizontal flight with a constant volume at essentially a constant 
altitude. The definitions which follow are designed to aid the 
reader in understanding the similarities and differences between 
the types of balloons to which the terms apply. For convenience 
the definitions are arranged in alphabetical order. 

CONSTANT ALTITUDE BALLOON - A helium-filled, zero-pressure balloon 
made of polyethelene for high altitude research. 

CONSTANT DENSITY BALLOON - A pillow-shaped balloon made of Mylar. 

CONSTANT-LEVEL BALLOON - A balloon which by means of ballast or 
superpressure is designed to operate at essentially a 
constant altitude. 

CONSTANT-PRESSURE BALLOON - A balloon of the type used in the 
transosonde system. 

CONSTANT-VOLUME BALLOON - A balloon whose volume remains essen- 
tially constant during flight. This condition is usually 
achieved by inflating the balloon to a pressure which is 
considerably greater than the ambient air pressure at 
the equilibrium altitude. 

CONTROLLED ALTITUDE FREE BALLOON - A balloon made of polyethelene 
with a control device to maintain constant altitude. 

HORIZONTAL SOUNDING BALLOON - A superpressure balloon made of Mylar 
designed for level flight for periods in excess of 60 days. 

LEVEL BALLOON - An expansible balloon made of rubber or Neoprene 
equipped with a valve in the neck to maintain constant 
levels at high altitudes. 

NEUTRAL BALLOON - A balloon with zero lift. 

PILOT BALLOON - A balloon made of rubber or Neoprene used for 
observing winds aloft. 

SUPERPRESSURE BALLOON - A balloon made of nonstretchable material 
(normally Mylar) so that its volume is essentially 
constant with excess internal pressure. 
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TETROON - A one-cubic-meter const ant -volume balloon of tetra- 
hedronal shape constructed of Mylar. 

TRANSOSONDE - A constant-level balloon equipped with meteoro- 
logical sensing instruments, a radio transmitter, and a 
power supply. 
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APPENDIX B 


Equation of Motion for Immersed Bodies 

The governing equation for the motion of a body submerged 
in a turbulent fluid consists of the equation for the conservation 
of momentum of the body coupled with the equation for conservation 
of momentum of the fluid. The governing equation has been developed 
in various forms 141,144-156,161,163,165,166 with various simpli- 
fying assumptions for a number of different applications, but there 
is not total agreement as to the validity or equivalence of all 
such forms. The development which follows is relatively general and 
the result is an equation which is representative of the type en- 
countered in the analysis of the motion of bodies immersed in a 
turbulent flow. 

Consider first the conservation of momentum for a spheri- 
cal body of density <*, diameter D, and mass m, slowly moving 
(Re<0.1)* with a velocity v^ in an otherwise stationary fluid in 
the presence of a gravitational field, as originally developed by 
Basset 144 . This equation can be written in the form 


dv . dv. 

m . . ■ + m — (E_). + 

dt a dt v D'i 

Vi + Vi + Vi 

(B-l ) 

where 

m a = | m 

(apparent mass) 

(B-2) 

(Pj,)^ -3U DWj 

(drag force) 

(B-3) 

(Fp)^ -m/a9p/3x i 

(pressure force) 

(B-4) 

<Vi = “ g i 

(gravitational force) 

(B-5) 

(Vi = 3/ 2 pl)2 fiT f 

v o 

dv i (t , )/dt' dt , 

— p= (Basset term) 

J t - t' 

(B-6) 


* This condition is not necessarily satisfied with a CVB 
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Before proceeding further it is important to note that 
the apparent mass term, m dv./dt', can be derived from inviscid 
flow theory and has the same value for both viscous and inviscid 
flow. The apparent mass term represents the time-rat e-of -change 
of fluid momentum due to inviscid effects. The Basset term is 
produced by viscous effects and thus can be interpreted as the 

r 

time-rate-of-change of momentum of the fluid due to viscous 
effects. Together the apparent mass and Basset terms produce an 
additional drag on the body which is a function of the time-rate- 
of-change of the relative velocity of the fluid with respect to 
the body. 

The next step in the development involves considering the 
case of a body slowly moving in a fluid with non-uniform velocity 
u i (t,x. ). The non-uniform velocity (t.x^) represents a turbulent 
flow process with the following assumptions [148]. 


1. The turbulence is steady and homogeneous. 

2. The domain of turbulence is infinite in extent. 

3. The body is spherical with a motion relative to the 
fluid which is characterized by a Reynolds number 
less than 0.1.* 

4. The body is small relative to the smallest turbulence 
wavelength present. 

5. While the body is in motion the fluid immediately 
surrounding it will be composed of the same fluid 
particles . 

6. The only external force acting on the body is produced 
by a gravitational field (or other potential field). 


The conservation of momentum equation for this case can be written 
dv. d(v.-u.) 


m dt 

+ m — ^ — = (F-O. + (F ). + (F ). + (F d ). 

a dt D'l p'i v g l B'l 

(B-7) 

where 



<Vi = 

-StrDu^ - u i ) 

(B-3) 

<Vi = 

-m/a 3p/3x^ 

(B-9) 


* The condition of small Re ( 0.1) is not necessarily satisfied 
with a CVB . 
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(B-10) 


( V± = mg i 

t 

<Vi ■ - I m 2 '"'/ 


h 


/ 


(t’> - u.(t') 

t -t r 


\/at' 


dt 


(B-ll) 


Notice should be taken that in apparent mass term, the drag term, 
and the Basset term the body velocity has been replaced by the 
relative velocity (v^ - u^), because in each case the resistance 
is a function of the relative velocity between the body and the 
fluid. 

The next step in the development of the body momentum 
equation consists of expressing the pressure force term as a 
function of the fluid velocity. For the case of an incompressible 
fluid with a velocity u^(t) which is a function of time but is 
spatially uniform, conservation of fluid momentum can be expressed 
as 

Du i 

P = -3p73x jL + pg ± (B-12) 


For the case of the non-uniform fluid velocity u^(t,x^) the instan- 
taneous equation for the conservation of fluid momentum can be 
written for the incompressible case as: 

Du, 3 2 u. 

0 or = " 3p/3x i + 0 3x - ' ax . + pg i 

0 0 


Notice should be taken that in equation (B-13) the instantaneous 
fluid properties are represented and not the time-averaged. Thus, 
although the flow may be turbulent, no turbulent or Reynolds 
stresses appear in the equation for momentum conservation. 


In his original development Tchen [145] neglected the 
viscous stress terms and effectively equated the follow-the-fluid 
particle derivative, Du^/Dt, to the follow-the-solid body 
derivative, du^/dt. Such a process corresponds to the simpli- 
fying assumption: 


3p_ - ^i 
a$T dt 


pg i 


(B-14) 
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Tchen used this relation to replace the pressure force term in 
equation (B-l). Corrsin and Lumley [146] first noted the inexact- 
ness of this approximation while Hinze [148] indicated that such 
an approach was valid if: 


3u 
v 3x 


<<1 


(B-15) 


and 


D 2 (3 2 u/3x 2 ) 


>>> 1 


(B-16) 


Hinze ' s development of Equations (B-15) and (B-16) is not 
altogether rigorous as indicated by the absence of subscripts. 
More precisely the two conditions are: 


D‘ 


18v 


(u k - v(Sj} (U i ' V i ) 


<< 1 


and 

v k (3u i /3x k ) 


>> 1 


2/3 v ( 3 ^ x j ) 

If the flow is near isotropic it is reasonable to assume: 


U 1 - V 1 = U 2 ~ V 2 * U 3 - V 3 


V 1 * v 2 - v 3 


s 

3x^ 3 x 2 3x 3 

and 


3 2 u i 3^u^ 3 2 u <i 


2 „ 2 
3x 1 3 x 2 


3x, 


(B-17 ) 


(B-18) 


(B-19 ) 
(B-20) 

( B-21 ) 


(B-22) 


With these assumptions equations (B-17) and (B-18) reduce to 
2 , , << ^ (B-23) 


D 
6v 


3u i /3x J 
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and 



'llj/cjx k 


>>1 


(B-24) 


Equations (B-23) and (B-24) are seen to agree with equations (B-15) 
and (B-16) except for the subscripts, the numerical constants, 
and the use of absolute values, all of which Hinze omitted. 


Because the spatial derivatives of the fluid velocity 
represent instantaneous values as opposed to time-averaged values, 
for the case of homogeneous turbulent flow, the order of magnitude 
of such derivatives can be expressed in terms of the scale and 
intensity of the turbulence as follows: 

I 3u ± /3Xj IrV^ri /x j ■ (B— 25) 


| 3 2 Ui /3Xj 2 | ^2 /Xj 2 (B— 26 ) 

where 

X. - micro-scale of turbulence. 

3 

A combination of equations (B-23) and (B-25) yields: 





<< 1 


while combining equations (B-24) and (B-26) produces: 



v 


. | A. : 

i 1 .1 


V u! 2 


>>>1 


Now 


(B-27 ) 


(B-28) 



(B-29) 
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Thus equation (B-28) can be reduced to 


< 5 i/ V 


u 


,2 


+ 1) A 


>» 1 


D' 


(B-30) 


For a turbulent flow field with known values of kinematic viscosity, 
turbulence intensity, turbulence micro-scale, and mean fluid 
velocity, equations (B-27) and B-30) provide a means of quanti- 
tatively determining the maximum size bodies for which Tchen's 
approximation of the pressure force is valid. 


For spherical bodies which satisfy the preceding restric- 
tions the particle momentum equation can be written 


dv. d(v -u. ) 

m -rr^ + m — — — 
dt a dt 


du. 

-3 7r D y (v i - u i ) + mp/a ^ - mpg ± /o 



d fv (t^-u-Ct’^/dt’ 

, dt’ (B-31) 

4 t - V 


Now 

m = -L D 3 a (B-32 ) 


and 


m 


a 



(B-33) 


Introduction of the two preceding relations into equation (B-31) 
produces 


dv . du^ ^ 

+ a v. = a u. + b -jr + c / 
dt i i dt J ^ 


where 


a = 


36y 


(2o + p) D' 


b = 3 p 


2a + p 


d ru.ct’j-v.ct'q/dt' 

- - dt’ + d g. (B-34) 

4 t - t* 1 


(B-35) 
(B-36 ) 
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c 


18 

(2a + p ) D 


(B-37) 




= 2 ~ P 1 

2 a + p 


(B-38 ) 


Equation (B-34) is the most familiar form of the immersed body 
momentum equation. It should be realized, however, that in 
this form the application of the equation is generally limited 
by the original assumptions, as well as the restrictions imposed 
by equations (B-27) and (B-30). 
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APPENDIX C 


Derivation of Dimensionless Equations 
for Conservation of Balloon Momentum 

As developed in Section 3.1.1 the equations for conser 
vation of balloon momentum in dimensionless form can be written 



dt 


" ^(r^-p ) 0 (fl i - V * I (fi i - VI ® - 


( p ^ 4 

p 

h: - v 

3u 

2 (1 - p) 


2 + p 

3x j " 

2 + P 

j ( p )nr 

t 

f 

d(v i - Ui) 
dt’ 

- dt’ 


3 V 2 + p r * 

J 

o 

yjt - t' 






(C-l) 


Now the density p is not a constant but varies with altitude, Xg. 
For an ideal gas 

P = (C-2) 

For the case of a constant temperature lapse rate, if, in a 
hydrostatic atmosphere. 


P 


d) 


g_ 

ry 


(C-3) 


A combination of Equations (C-2) and (C-3) produces 



(C-4) 
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where x 3Q is the equilibrium altitude at which o =p Q . If 
Equation (C-4) is divided by o the result is 



p o r, 

Y , v ,1 RY " X 


p = 

— L 1 ■ 

' (X 3 *30 J j 


= 

£l - a A 

- 1 b-1 
X 3j 

(C-5) 

where 

a = 

YD 

T 

o 


(C-6) 

AX 3 

" x 3/° ” 

x Sq/ d 

(C-7) 

b = 

g/(RY) 


(C-8) 


Now Equation (C-5) can be substituted into Equation (C-l) but the 
result is quite cumbersome. A term by term examination of the 
R. H. S. of Equation (C-l) reveals that only in the numerator of 
the fifth term, representing buoyancy, does the variation of p 
appear significant. If Equation (C-5) is substituted into only 
this term in Equation (C-l), the result is sfill awkward because 
of the nonlinear form of Equation (C-5). To avoid this difficulty, 
by means of series expansions, 

(1 - a Ax 3 ) b_1 = 1 - (b - 1 ) a Axg + 0 (A 2 x 3 ) (C-9) 

Neglecting second order terms, 

P - 1 - (b - 1) a Ax 3 (C-10) 

If Equation (C-10) is substituted into the numerator of the fifth 
term of the R. H. S. of Equation (C-l), the result is 
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Equation (C-ll) represents a dimensionless form of the balloon 
equations from which certain dimensionless parameters can be 
obtained. A term by term inspection of the R. H. S. of Equation 
(C-ll) results in the following dimensionless groups: 



(C-12 ) 
(C-13 ) 

(C-14) 

(C-15) 

(C-16) 

(C-17) 


In terms of the six preceding dimensionless groups. Equation (C-ll) 
can be written 
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dv 


-r - 3L 1 L 3 <W + 3L 1 L 4 
dt 


Du. 


3 U . 


3L 1 -T + L 1 ( W ~ " 3L 2 L 5 Ax 3 S 13 

Dt J J 


3Xj 


-3 Ll L 6 f d( W dt 

n j* » 


0 dt 

VI i" 


(C-18 ) 
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APPENDIX D 

Dimensional Inviscid Flow Field Development 

The velocity components are assumed to be of the form 
vl^ - + A sin + kg + k g Xg - cut + Q^) (D-l) 


Ug = B sin (k^ x 1 + k 2 x 2 +k 3 x 3 -wt + 0^) 


u 3 = C sin (k^ x^ + kgXg + kgXg - uit + Qg) 


For brevity let 

C = k 1 x + k 2 y + k 3 z 

Then 

3u- 

-^= Ak^ cos (£ - wt + 0^) 

= E &2 cos (£ - wt + 0 2 ) 
3iu 

-g— = Ck 3 cos (£ - wt + 0 3 ) 

u 

For mass conservation, 

^ ^3 __ 

3x- ax- ax- u 


**2 9X 3 


(D-2 ) 
(D-3) 

(D-4 ) 

(D-5 ) 
(D-6 ) 
(D-7) 

(D-8) 


Thus 


A^ cos (£ - wt + 0^) + Bkg cos (C - wt + 0 3 ) + cos (£ - wt + 0 3 ) 
- 0 


(D-9 ) 


Ak^ (cos ( £ - wt) cos 0^ - sin (£ - wt) sin 0^> + Bkg (cos (£ - wt) cos 

- sin (C - wt) sin 0 2 ) + Ck 3 {cos (£ - wt) cos 0^ 


- sin (£ - wt) sin © 3 } = 0 


To satisfy Equation (D-10), it is necessary for 
(Ak^ cos 0^ + Bkg cos 0g + Ckg cos © 3 ) cos (£ - wt) = 0 


(D-10) 


(D-ll) 
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and 


(Afc, sin + Bkg sin 0 2 + Ck^ sin 0^) sin (5 - wt) = 0 (D-12) 

This in turn requires 

Ak^ cos ©^ + Bkg cos 0g + Ck^ cos ©3 = 0 (D-13) 

and 

Ak^ sin 0^ + Bkg sin ®2 + sin 0 3 ~ 0 (D-14) 

Equation (D-13) can be multiplied by tan 0^ to obtain 

Ak^ sin 0 ^ + 0k 2 cos ©2 tan 0 ^ Ck^ cos © 3 tan 0 ^ 0 ( 0 ™ 15 ) 

Subtracting (D-15) from (D-14) yields 

Bk 2 (sin © 2 - cos © 2 tan 0^) + Ckg (sin © 3 - cos 0 3 tan 0^) = 0 (D-16) 


Solving for B, 


B = 



(sin © 3 - cos © 3 tan 0 ^) 
(sin 0 2 cos 0 2 tan 0 ^) 


(D-17) 


Next Equation (D-13) can be multiplied by tan 0 g to obtain 

Ak^ cos Q x tan © 2 + Bkg sin © 2 + Ckg cos © 3 tan © 2 = 0 (D-18) 

Subtracting (D-18) from (D-14) yields 

Ak^ (sinQ^ - cos©^ tan © 2 + Ckg (sin0 3 - cosQg tan® 2 ^ = 0 (D-19) 


Solving for A, 

Ckg (sin0 3 - cos0 3 tan0 2 ) 

k^ (sin - cos 0 ^ tan © 


A = - 


(D-20 ) 


Thus if Equation (D-17) and (D-20) are satisfied, the inviscid 
flow field model described in Equations (D-l) through (D-3) 
satisfies the continuity equation. 

Certain other partial derivatives of velocity with respect 
to space and time can be derived from the model as follows: 

3u- 


3x, 


Ak, 


cos ( £- wt + 0 ^) 


(D-21 ) 
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= Ak 3 cos (£ - wt + ■ 0 1 ) 


(D-22) 


3u. 

3xv 


3u, 

■ d 

3xi 


Bk^ cos ( C - wt + Og j 


3u, 

4 

3x3 


= Bk„ cos (£ - mt + 0 Q ) 


3u, 

3xi 


Ck 1 cos (£ - (^t + 0g) 


3u 

3x 

3u 

at 


— = Ck 2 cos (£ - wt + 0g) 
2 


— = _A to cos (£ - cot + 0^) 


(D-23 ) 
(D-24 ) 
(D-25) 
(D-26 ) 
(D-27 ) 


3 u o 

= -Bco cos - cot + 0 2 ) 

3Uq 

Oo cos ( ^ cot 0g ) 


(D-28) 

(D-29) 


As formulated the model contains a total of eleven 
different parameters (u^, A, B, C, k^, k 2 , k 3 , to , 0^, 0 2 , and 0g). 
As given by Equations (D-17) and (D-20) for conservation of mass, 
two of the three amplitudes (A, B, and C) must be functions of 
the third. Thus there are nine independent parameters in Equations 
(D-l) through (D-3). In the current study, it was not feasible 
to consider all possible combinations of these nine parameters. 

The following simplifications were made: First, all wave numbers 

(k^, k 2 and k 3 > were taken to be equal, 

k x = k 2 - k 3 - k (D-30 ) 


Second, by proper choice of the phase angles (0^, 0 2 , and 0 3 ) the 
amplitudes were taken to be equal. 

A = B = C (D-31 ) 

It can be readily shown by means of Equations (D-17) and (D-20) 
that the phase angles, 

0 1 = 0 

0 2 = 2tt/3 
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(D-32 ) 
(D-33 ) 


Og * -2ir/3 


(D-34) 


or any similar combination of phase angles set 120° apart, satisfy 
Equation (D-31). 

Based on the simplifications noted, the fluid velocity 
components can be written as 

u i * 6 il + A sin (C - <ot + 0 ± ) (D-35) 

where 

5 * k (x ± + x 2 + x 3 ) (D-36) 

* 0 (D-37) 

©2 * 2ir / 3 (D-38) 

© 3 -2tt/3 (D-39) 
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APPENDIX E 


Dimensionless Inviscid Flow Field Development 

As developed in Appendix D, the wind velocity components 
can be expressed as: 

u^ * u^ <5^ + A sin ( C- cot + © i ) (E-l) 

where 


C = k (Xj + x 2 + x 3^ 
0 ! " 0 


© 2 = 27T/3 
0 3 =-2tt/3 


(E-2) 

(E-3) 

(E-4) 

(E-5) 


Based on the nondimensional procedure followed in 
Appendix C, Equation (E-l) can be written in the following 


dimensionless form: 

Ui - u 1 6^ + sin ( C - St + ©^) (E-6) 

where 

u A = u./A (E-7) 

u = u^/A (E-8 ) 

C = k + x 2 + x 3 ) (E-9) 

S = wD/A (E-10) 

t = At/D (E-ll) 

k = kD (E-12) 

x ± = x i /D (E-13 ) 


Equation (E-6) represents the dimensionless form of the inviscid 
flow field. Notice should be taken that in dimensionless form 
the flow field can be described in terms of three dimensionless 
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groups, u^, k, and £3. The dimensionless variables 3L and t are 
not necessary for characterizing the flow field because they 
represent the coordinates of the balloon in space and time (to 
be obtained from the solution of the differential equation 
governing the balloon motion). 
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APPENDIX F 


as : 


Reynolds Number Relation to Dimensionless Groups 

The Reynolds number for flow past the balloon Is defined 


Re » p l tt ~ v l D (F-l) 

Now, the maximum difference between the wind velocity and the 
balloon velocity must be roughly equal to the amplitude A in 
the inviscid flow field model. Thus, 




A 


(F-2) 


Therefore , 

Re * £^2- (F-3) 

An examination of the first dimensionless groups given in Appendix 
C reveals that to first-order accuracy. 

Re * 12/Lg (F-4) 


where 


L 


3 



(F-5) 
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APPENDIX G 
Dimensional Analysis 

Dimensional analysis is a standard technique for reducing 
the number of variables in a problem and for identifying the 
important parameters associated with the problem. This is 
particularly useful when one has a physical problem the variables 
of which are well known but for which no analytic relationship 
is known. But, it is also useful for determining' the necessary 
nondimensional groups for a given set of analytic relationships 
describing a physical process. 

For the analysis of a balloon in a realistic atmosphere, 
this analysis is presented as an extension of the development of 

the nondimensional equations discussed in Appendices C and E. 

The pertinent variables for the analysis are: 

v^ balloon velocity 

x^ balloon position 

u i ambient wind velocity 

t elapsed time 

p ambient density 

y dynamic viscosity 

g gravitational acceleration 

o density of the balloon 

D diameter of the balloon 

C D £ drag coefficient (10^ <Re< 10 5 ) 

The first two variables, v^ and x^ represent the unknown properties 
of the balloon. These two variables are not independent of each 
other since v^ is simply the time deviative of x^ . The variable 
u i , representing the ambient wind velocity can be expressed as 

u i = u i (u-^Ajk.co) (G-l) 
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(G-2 ) 


The ambient density, p, can likewise be expressed as 
P = P(P Q , Y,T o ,g,R,x 3 ) 

The viscosity, y, can also be expressed as 

V = y(y Q , y,T q ,x 3 ) (G-3) 

Thus there are 16 variables (v i ,x i ,u 1 ,A,k,w,t ,p Q , Y,T o ,R,y o ,g,a,D,C m ) 
involved in the dimensional analysis. The primary dimensions 
involved are mass (M), length (it), time (t) and temperature (T). 

The dimensional matrix for this case is as follows: 

v. x. A k a) t p Q y T q R y Q g c D C pjl 

M 0 0000001 00 0 ^ 0100 

l 1 1 1 1 -1 0 0 -3 -1 0 2 -1 1 -3 1 0 

t -1 0 -1 -1 0 -1 1 0 0 0 -2 -1 -2 0 0 0 

T O 0000000 11 -1 0 0000 


According to the Buckingham IT theorem, with 16 variables and 4 
primary dimensions, there will be a maximum of 12 dimensionless 
groups. The 12 dimensionless groups take on a variety of forms 
depending on the choice of primary variables. If A, p Q , D and 
T q are chosen as the primary variables, then the following dimen- 
sionless groups result: 


n i = v l/A 

n 2 = Xi /D 

II 3 = U-^/A 
n 4 = wD/A 
JI 5 = At/D 

H 6 = yd/t o 
n 7 = RT o /A 2 

n 8 = V (Ap o D) 

n 9 = gD/A 2 

n i0 = °/ p o 
n ll = kD 
n i2 = C DA 


( G-4 ) 
(G-5) 
(G-6) 
(G-7) 

( G-8) 
(G-9) 
(G-10) 
( G-ll) 
( G-12 ) 
(G-13) 
( G-14) 
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(G-15 ) 



Some of the preceding dimensionless groups are identical 
to or closely resemble the dimensionless groups derived in 
Appendices C and E of this report. Others, however, are not so 
easily identified. The relationships between the two sets of 
dimensionless groups can be developed as follows: 


V. = 
1 

n i 

(balloon 

velocity) 

(G-16) 

II 

•H 

n 2 

(balloon 

position) 

(G-17) 

t = 

n 5 

(elapsed 

time) 

(G-18) 

L 1 s 

2l 

1 

10 + 1 


(G-19) 



n io 



L 2 a 

2Tf 

O 

+ 


(G-20) 

L 3 £ 

12 

n 8 


(G-21 ) 

L 4 S 

1 

2 

n i2 


(G-22) 

L 5 S 

2 

3 

n 9 ( " 

n 6 ) 

( G-23 ) 

L 6 ^ 

6 

^S/ir 


(G-24) 

k = : 

"ll 



( G-25 ) 

£> = : 

n 4 



(G-26 ) 

5 i = 

n 3 



( G-27 ) 


The use of the approximate equalities in Equations (G-19) through 
(G-22) and in Equation (G-24) reflects the fact that through 
L 4 and Lg contain the variables p and y while n g , n iQ and 
contain their equilibrium counterparts p Q and y Q . An examination 
of Equations (G-16) through (G-27) reveals that all 12 dimension- 
less groups (n x through H.^) appear in the definitions of the 
dimensionless terms which arise in the dimensionless differential 
equation for balloon motion and the dimensionless inviscid flow 
field model. Furthermore, as previously noted, and 3L (H ^ and 
n ^ ) are not independent of each other but represent the solution 
of the differential equation for the balloon equation. The 
variable t (II 5 ) is obviously dimensionless time, with respect to 
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which both 5L and v i can be expressed. The remaining eight dimen- 
sionless groups can be reduced to five if the ratio p/a is taken 
as unity and if Cjj -is taken as constant.* Thus, with this simpli- 
fication the entire problem can be Characterized in terms of 
five dimensionless groups: . 


12v 

AD 


2 g(g/R-V)D 2 
3A 2 T 

o 



, and 



The solution in terms of dimensionless balloon position (x^yD) 
and velocity (v i /A) as functions dimensionless time (At/D) 
should remain unchanged (to first-order accuracy) provided these 
five groups remain constant. 


* This reduction from eight to five dimensionless groups is possible 
in part, because of the manner in which Ilg, n^, Ilg, n g , and JI^ 
occur in the governing equations. 1 J 
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APPENDIX H 



... Balloon Motion in Turbulent Plow 


The problem under consideration is a constant-volume 
balloon floating in a two-dimensional turbulent flow field with 
stationary turbulence which is homogeneous in the x-direction. 
The coordinate system is shown in Figure H-l. 

v BALLOON INITIAL BALLOON 

3 POBITION PATH 



Figure H-l. Coordinate System for First-Order Perturbation Analysis 

The longitudinal component of the wind velocity can be expressed as 


u ± (x.,,t,x 3 ) * (x 3 ) + (x 1 ,t,x 3 ) (H-l) 

The fluctuating portion of the longitudinal component can be 
expressed in terms of a Fourier integral. 


u’ ( X- , t , x« ) 


= f f B (k,w;x 3 )e 1 ( kx i~wt) dk dco 

%/-00 */-co 


The coordinates of the CVB are functions of both time and the 
initial balloon position. Because the turbulence is homogeneous 
in the x-direction and the mean wind velocity is assumed to be 
constant with respect to x^ and t, only the x 3 ~coordinate of the 
initial position is significant. Thus 


(Balloon Coordinates) = ^ X 1 (t ,X 3Q ) ,X 3 (t ,X, 


(H-3 ) 


The horizontal component of the balloon velocity can be 


- - “i [xjCt.XgoJ.t.Xg^.Xao)] 
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Equating the horizontal component of the balloon velocity to the 
corresponding wind velocity component at the balloon location is 
based on the assumption that the balloon responds perfectly to 
the wind velocity in the horizontal direction. 

The vertical component of balloon velocity can be 

written: 

v 3 (t ' x 30 ) - TT (H - 5) 


Notice should be taken that v g cannot be set equal to vertical 
component of the wind velocity, u^, because a constant-volume 
balloon cannot be perfectly responsive in the vertical direction. 


The balloon coordinates can be written as: 


X 1 ( t,X 30^ 


X 30> dt 


“/ V 1 (t - 

•'o 

= (X 3Q ) t + / (X^t.Xg) dt 


- X 1 (X 30 ,t) + X^(t,X 30 ) 


and 


■/ 


X 3 (t ’ X 30> =7 v 3 (t > X 30> dt + X 
o 


30 


(H-6) 


- X 3 (t > X 30 ) + X 30 (H ~ 7 > 

The fluctuating portion of the longitudinal component of the 
balloon velocity, v£, can be expressed in terms of a Fourier 
integral analogous to that for the wind velocity component, u^. 
Thus : 

v i ( t,X 3(P = V 1 ^ t,X 30^ ” ^ X 3^ 


f f i fk(X 1 + X’)-ut] 

= / / B (k,u;X 3 )e x dkdw (H-8) 


00 00 
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Expansion of B (k,w;X 3 ) in a Taylor series about X 3Q yields 

B (k,*;^) = B Oc.MiXgo) + || <V*»> * § 0 (V*30> 2 + — 

= B (k,(D;3Cj 0 ) + § X^ + | 


3Z 


ikX. 


o%r + ... 


(H-S) 


Likewise the exponential series for e 1 ca n be written 


e ikX i = 1 + ik | (k X|) 2 +... 


(H-10) 


A combination of Equations (H-8) through (H-10) yields 

H <*’**>> -££[* + 1 * 3 + I 5 < X 3) 2+ ---][ 1 


+ ikX[ - | (kX^) 2 + . . . J e l(kX l ^ dk 


dw 




To first-order accuracy* 


oo oo 

v' (t.Xg.) = f f B (k.wjXg ) e i(kX l' ut) dk 

•Coo •Coo 


do) 


(H-12 ) 


In similar fashion the velocity component (t +t , X 3Q ) 
can be expressed to first-order accuracy as 

(t + t »X j Q ) = y*y BCk.iojX^) e 1 [ fe(X i +;L L T) " “ (t + T) J dk do) (H-13) 

Because v ^ is real, the Fourier integral representation for 
(t + t »X 3 q) can be written as 


v i (t + T - X 30> 


00 

-// 

— 00—00 


B (k,w;X 30 )e 


-X 0(5^ + ^ T) - u(t + t)^ du 

(H-14) 


* The first term in the series expansion of B(k,^;Z) is first-order 
because B(k,w;Z) itself represents a perturbation and is first- 
order. Thus the only first-order term in the product,-.. ikXl 

is the product of the first terms of the two series. 
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The product v^(t,X 30 ) v£(t + t,X 3Q ) can be expressed as 

Vi(t,Xg 0 ) V {(t + T.Xgo 5 

= / f f f B( k ',<o';X 30 ) B*(k,u;X 30 ) 

—oo —co an — co 

e -i [k(X x + u ± T) - w(t + T)] dk . da) r dk dw (H- 15 ) 


Equation (H-15) holds for each realization of the ensemble. Thus 
the ensemble average representing the quasi-Lagrangian time auto- 
correlation function, can be written as 


<y{ (t,x 30 ) vi (t + t,x 30 )> 

OO 00 00 CO 

= /”/■// <B(k',o';X 30 

-Zco — OO OO OO 

e -i [x 1 (k , -k)-t( 03 ’-a))] e ~i(k u 


) B* (k,o3 ;X 3 q)> 
i T ~ “ T > dk’du'dk du 


(H-16) 


A more useful form Of Equation (H-16) can be obtained by the 
following change of variables: 


k' ' = k - k' 

0) ’ ’ = 03 - 0) ’ 


k ' = k - k' ' 

03 ' = 03 - 03 ’ ’ 


(H-17) 


Introduction of k' 1 and oj 


v-j (t ,X 3Q ) v{ (t + t,X 30 ) 


30 


00 00 00 OO 

■ / / / / 


—00 «— OO —00 


-i(k' - w"t) g -i(ku 1 


in Equation (H-16) produces 

,uj-o>' ' ;X 3Q ) B*(k,ai;X 30 )> 
T “ WT) dk’ 'dm' ’dk do) 


(H-18) 
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The general Eulerlan space-time power spectrum, 

<f>g(k, tu , , t ; Xg ) , of the horizontal wind velocity u^,(x 1 ,t,x 3 ) 

by definition is: 


r°° r°° * 

4> E (k,^,Xi ,t ;Xg) ~ J J <B(k-k ? ' ' ;Xg) B (k,u);Xg) > 

•LDO ^-iX) 


Kk'V't) ^ 


'dw’ 


(H-19) 


At the balloon coordinates (X 1 , X 3Q ) the general power spectrum 
of the wind velocity component ^ (X lf t,X 30 ) would be 


♦e (k 


= J°° f <B (k-k' - ,u-ui' ' ;X3 0 ) B* (k,ui; Xg 0 )> 


e Kk"^ - «' -t) 


dk' ’do)’ 1 


(H-20) 


For the case under consideration the turbulence is homogeneous in 
the x^-direction and stationary, and thus the corresponding 
Eulerian space-time power spectrum is not a function of X^ or t. 
For such a case, 

(k,u>,X^, t;Xg Q ) •> (k,co; Xgg) (H-21) 


A combination of Equations (H-20) and (H-21) produces 

7 ^00 ^*oo 

J <B(k-k” ,0)-<0' ';X- ) B Ck > o ) ;}^ 0 )> 

__oo •'.oo 

e i(k-X 1 -0)-t) dkfldulf 


( H-22 ) 


In order for the RHS of Equation (H-22) to be consistent with the 
LHS 

. < B(k-k f 1 ,03 — uj* 1 ;Xg 0 ) B*(k,ai;X 30 )> = F^^X^) 6(k”) 6 (oj' 1 ) (H-23) 

where 

F (k^jX^) - < B(k,oi;X3 0 ) B^k.oi;^) > (H-24) 
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A combination of Equations (H-22) through (H-24) yields 


♦e (k^X^) “ f F (k.miXg 0 ) 6(k") 6(o>" ) 


e l(k’ ’X^ - w f ’t) ^rr^r. 
= F (k,oj; X3 Q ) 


and 


^ (k^jX^) <r < B (k,di;X 30 ) B (k.wjXgg) > 


(H-25 ) 


(H-26 ) 


The quasi-Lagrangian time auto correlation for stationary 
turbulence with homogeneity in the x-direction can be developed 
by a combination of Equations (H-18 ) , (H-23 ) and (H-25) as follows: 


< v£ (ti'Xgo^ v i ^ + t;X 30^ > Jf* ^ k,w;X 30^ 

e "i (kii L - “> T dk dm 


The LHS of Equation (H-27) must be independent of time because, 
due to stationarity the RHS is not a function of time. Thus, 


< v i <*5*30) v i <* + *;*»> > — = ■ Hl <’•*»> 


(H-28) 


Thus the quasi-Lagrangian time auto correlation function for 
stationary turbulence, with x^-direction homogeneity, can be 
expressed as : 

\ (t;^ 0 ) m f f *e ( k >«; X 30 ) 6 _i(k51 ~ “ )T (H-29) 

• 4-00 * 4-00 ^ 

Equation (H-29) can be cast in a more convenient form 
by means of a coordinate transformation where 


K - k | 

S2 = ku x - w j (H-30) 
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The quantity (-ft) is the doppler frequency observed in a coordinate 
system moving with a velocity u^. In K, ft - space Equation (H-29) 
becomes 


^ t;X 30^ J[ < K > KU 1 n ;X 30 * 


-iftx 


9k 

9k 

9K 

9ft 

9oj 

9o) 

9K 

9ft 


dK dft 


- rr *e < k - k5 i - e 

—CO ~ +oo 


-iftx 


1 0 
"i - 1 


dK dft 


= f f <(u (K.KSj -ft;Xg 0 ) e _lSlT dK dft 

*^-00 ''-CO 


(II- 31) 


At the same time, based on generalized Fourier integral theory, 

- iftx 


*L (t;X 30> = /*»L(a;^o) 


dft 


(H-32 ) 


where V (^; X 3 q) quasi-Lagrangian time power spectrum. Thus 

by comparison of Equations (H-31) and (H-32), 


\ (ft,Xg ) = f ^ (K, - ft;Xg ) d K 

i * 9 / j0O ° 


(H-33 ) 


Notice should be taken that Equation (H-31) represents a relation- 
ship between the quasi-Lagrangian time auto correlation function 
and the Eulerian space-time power spectrum for stationary turbu- f. 
lence with x-direction homogeneity, while Equation (H-33) relates 
the quasi-Lagrangian time power spectrum with the same Eulerian 
space-time power spectrum. Both relations are based on first- 
order perturbation theory. 

The true significance of Equations (H-31) and (H-33) can 
be ascertained by applying them to certain situations involving 
both nondispersive and dispersive media. First, in a nondispersive 
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media Taylor's hypothesis holds and thus, 

w = k 
or 

£2 = 0 

For this case Equation (H-33) becomes 

(TAYLOR) = ^E (TAYLOR) K **l " Q;X 30^ dK 

= £ * E (TAYLOR) *< fl > * 

6(Q) L $E ( TAYIj0R ) (K;X 30> dK 


(H-34) 


(H-35) 


(H-36 ) 


where $>£ (K,X^q) is the Eulerian space power spectrum Equation 
(H-36) indicates that the quasi-Lagrangian time spectrum based 
on Taylor’s hypothesis is simply a spike at £2 = 0. A combination 
of Equations (H-31) and (H-36) yields 


H L (TAYLOR) 



•/! $ E (TAYLOR) (K; W ^ 
^ E (TAYLOR) (K i ^ 



(H-37) 


Notice should be taken that the RHS of Equation (H-37) is inde- 
pendent of t . Thus the quasi-Lagrangian time auto correlation 
function is independent of t and 


^(TAYLOR) ^ T ’ *»> H L (TAYLOR) ^*30^ 


(H-38 ) 


Equation (H-38) combined with Equation (H-28) indicates that 


H, 


L (TAYLOR) (X 30 ) < V 1 **30* V 1 (X 30 ) > 


«° Vl <*30> 


(H-39) 
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Based on Equation (H-39) the quasi-Lagrangian time auto correlation 
function (for the case of Taylor’s Hypothesis), for the balloon 
velocity component v£(t,XgQ), equals the variance of the same 
velocity component . The latter in turn equals the variance of the 
^component of the wind velocity. These relations are based on 
the first-order perturbation analysis with the assumption of a 
perfectly responsive balloon (in the x-direction) with stationary 
turbulence which is homogeneous in the x-direction. 

For the case of a dispersive media, Taylor's hypothesis 
does not hold. Instead 

to ** u 1 k + F.(k) (H-40) 

or 


a = - F (K) (H-41 ) 

The inverse of Equation (H-41) proves more useful and can be 
written 

K = f (ft) (H-42 ) 


A combination of Equations (H-33) and (H-42) yields 
V L (DISPERSIVE) ^ ;X 30^ = f ^ E (DISPERSIVE) ^ K ’ K **1 ** 


jf * E (DISPERSIVE) (K;X 30 ) K] dK 

$ E (DISPERSIVE) W;X 30^ 


Thus in a dispersive media the quasi-Lagrangian time power spectrum 
can be related to the corresponding Eulerian spatial power spectrum 
if the function f ( £2) is known. Furthermore, the quasi-Lagrangian 
time auto correlation function can be related to the same Eulerian 
spatial power spectrum by a conbination of Equations (H-32) and 
(H-43) as follows: 

*L < T ’ X 30 ) = X ^(DISPERSIVE) [ f(0 > ;X 3o] e (H-44) 


119 



For dispersive Rossby waves the phase velocity c is 
given by the relation 

c — “ B/k 2 (H-45) 

By definition 

c = w/k (H-46) 

and thus 

to = u 1 k - g/k (H-47 ) 

or 

C2 = g/k (H-48 ) 

Thus for a Rossby wave, based on Equation (H-42) and (H-48), 

f (ROSSBY) ^ = (H-49) 


Then 

y L (ROSSBY) $ E (ROSSBY) ^ n;X 3(P (H-50) 

The relation between $ E(ROSSBY) ^ X i X 30^ an< * X is known to be °* 
the form 

* E( ROSSBY) ( K = X 30) « K ^ <*-51) 

where 

p i <0 

Therefore 

^ L(ROSSBY) ^’ X 30^ a ^ 1 (H-52) 

For dispersive gravity waves, 

w= k ± /gk (H-53 ) 

or 

ft = + /gTc (H-54) 

Thus, based on Equation (H-42) and (H-54), 

"^(GRAVITY) ^ Q / g (H-55) 
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In this case 


¥ L (GRAVITY) ( ft;X 30^ ** **E (GRAVITY) / g,X 3(P 


(H-56) 


The relation between (gravity) (K,X 3Q ) and K is of the form 


*‘E (GRAVITY) (K ’ X 30> * 


(H-57) 


where 


P 2 < 0 

Therefore 

2P 2 

y L (GRAVITY) <«' X 30> " (H-58) 

Rossby waves possess long wave lengths and thus small 
wave numbers while gravity waves have relatively short wave 
lengths and hence larger wave numbers. A plot of (K; X 3Q ) 
for these two types of waves is given in Figure H-2. The 
corresponding plot of ¥ L (ft;X 3Q ) is shown in Figure H-3 . 
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Figure H-3 . Lagrangian Time Power Spectra 
for Rossby and Gravity Waves 



APPENDIX I 


NOTES ON SPACE-TIME SPECTRUM <f>_, (k, co : X™) 

^E 30 

The function <J) E (k, co ; X 3() ) represents the Eulerian 
Space-time spectral density at altitude X 3Q of turbulent eddy 
with wave number k and frequency w . If no relation between k 
and a) exist, the general variation of <J> with k and oj might be 
pictured as shown in Figure 1-1. 

0 E «c f w;X 3O ) 



Figure 1-1. General Eulerian Space-Time Power Spectrum 

If, however, k and m are related according to Taylor's hypothesis, 

w = ku 1 (1-1) 

^E(TAYLOR) as a ^ unct:i - on of k and w , appears as shown in 
Figure 1-2. 
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Figure 1-2. Eulerian Space-Time Power Spectrum According to 
Taylor's Hypothesis 


Notice should be taken that, with increasing u^ , the plane con- 
taining 4> e (TAYLOR) k ^l ; X 30^ rotates CCW . 

For the case of dispersive media, 

w = ku x + F(k) (1-2) 

For this case, as a function of k and to appears as shown in 

Ji 

Figure 1-3. 
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(k, co; Xjq) 



0, (k, kit. + F(k); X~») 

^E( DISPERSIVE) 1 30 


SAI-0941 


Figure 1-3. Eulerian Space-Time Power Spectrum for 
Dispersive Media 


The preceding figure represents the general concept of the spectral 
density function for a dispersive media. Two specific examples 
are useful. For a Rossby wave, 

“ = k5 l ~ k - (1-3) 

Although Rossby waves are limited to relatively low frequencies, 
and thus low wave numbers, it is beneficial to picture the 
variation of 0 ^ (roSSBY) k and to over a wide range of 

values of k and oj as shown in Figure 1-4. 
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0 (k, o> ; X 3 q) 



Figure i-4. Eulerian Space-Time Power Spectrum 
for Rossby Waves 


Notice should be taken in the preceding figure that the heavily 

shaded portion of the plots of g (ROSSBY) X 3(P 

represents the regions where Rossby waves would most likely be 
present . 

For gravity waves 


w = u 1 k + /gk (1-4) 

Gravity waves normally occur at relatively high values of k, but 
as before it is beneficial to plot the curve over a range of values 
of both k and m. The result is as shown in Figure 1-5. 
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k 



co* kuj —\J gk 


SAI-0943 


Figure 1-5. Euler ian Space-Time Power Spectrum 
for Gravity Waves 


As before the more heavily shaded portions of the plot indicates 
the part of the spectrum in which gravity waves most likely will 
occur . 
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APPENDIX J 


Notes on Space-Time Spectrum (K, ft; X 3 q) 

/ 

It is desirable to introduce a two-dimensional space- 
time spectrum ^ (K, ft; X 3Q ) which satisfies the relation 

y (ft; X 30 ) = f (K, 8: X 3Q ) dK (J-l) 

too 

The function ^ (K,-ft; X 3 q) is the counterpart of (k, u>; X^^). 
Each of the five figures given in Appendix I could be redrawn to 
represent (K, ft; X 3 q) as a function of K, ft. For the sake 
of brevity only the three specific cases, illustrated by Figures 
1-2, 1-4 and 1-5, will be developed. 

For the case involving Taylor’s hypothesis, 

ft = 0 ( J-2) 

as a function of K and ft is shown in Figure J-l. 



Figure J-l. Lagrangian Space-Time Power Spectrum 
According to Taylor's Hypothesis 
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A comparison of Figures 1-2 and J -1 reveals that as 
approaches zero, 4* - g (TAYLOR ) > ^30^ approaches 4*- g (TAYLOR ) 

( k , (o , Xg q) . 

For the case involving the Rossby wave. 

Si = 3/K ( J-3) 

or 

K = S/Si ( J-4 ) 


The variation of ^ g (ROSSBY ) ^ K ' ^ * X 30^ with K and ^ is shown i n 
Figure J-2. As before the more heavily shaded regions indicate 



Figure J-2. Lagrangian Space-Time Power Spectrum 
for Rossby Waves 
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that portion of the spectrum where Rossby waves are more likely / 
to occur. A comparison of Figures 1-4 and J-2 indicates that as / 
5 1 approaches zero, * L(H0SSBY) ' _(*. X 30 ) approaches 4> E(R0SSBY) 

(k, oj; X 30 ) / 

In addition, inspection of the shaded portion of Figure J-2 
repeals that in that portion of the spectrum where Rossby waves 
occur, ^ l (ROSSBY ) increases with increasing ft while, as. shown 
in Figure 1-4, <f> E (ROSSBY ) decreases with increasing k. This 
provides a slightly better understanding of the relationship 
between $ E ( R0SSBY ) and (ROSSBY ) which was originally presented 

in Appendix H. 

For the case involving gravity waves, 


ft — + 



( J-5) 


and 



(J-6) 


The variation of ^‘ L (GRAVITY ) ^ k ’ 
in Figure J-3. 


ft; with K and ft is shown 


The more heavily shaded portion of the spectrum represents the 
region in which gravity waves normally occur. Comparison of 
Figures 1-5 and J-3 reveals, that as u.^ approaches zero’, ip GRAVITY ) 

(K — ft; X 3Q ) approaches (p GRAVITY) ^30^’ Furthermore, 

although difficult to determine quantitatively from Figure 1-5 and 
J-3, it is possible to observe qualitatively that ^l(GRAVITY) 
decreases more rapidly with increasing ft than ^ E ( GR avitY) decreases 
with increasing k. 
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Figure J-3. Lagrangian Space-Time Power Spectrum 
for Gravity Waves 
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APPENDIX K 


/ 


Source Listing for BALLOON Program 

..... / 

Balloon is written in FORTRAN IV for the Univac il08 
Computer. The listing which follows represents the current 
form of the program. 


COMMON/CONS/Pl. D.SIC.M A, AP.M. VOL » 

♦MUREF.TEMREF, TEMPO, TEMLAP.X3n.R,GSRTL»Gf3).PO,X33 

C0MMQN/tEKMS/F0RA<,<3)«FAPi<3).FAP2(3)»FB-UoY(3).FbASS(3), 
*MA.NU,MU,RH0,TfeMP t x(3).FACT0R,KE 
C0MM0N/VEL0C/U13),V<3).UMAG.VMAG.MAGDIF 
C0MM0N/INTL/K(3 ).ThETA(3),0MEGA,Uh<3).A,B,C 
dimension fu(3),Fvo)«vdu(3) 

DIMENSION X D ( 5 J « X n ( 7 ) 

DIMENSION aka ( SO • 1 I ) . WRfl ( SO . 6 ) i WKC 150*17) 

DIMENSION BCD X < l 2 ) ,RCDYI (12) »BCoy2< 1 2 ) • B C D Y 3 ( \ 7 ) . F.L&Y I ( 1 2 ) * FLOY 2 < 

M2) • F L D Y 3 ( 1 2 ) 

DIMENSION PT(Sr?G),PX<5nn),PY<50r),PZ<SUO)iPYl(e,00) , PY2(5UG) , PY3(t>0 
10)»PZM50O),PZ2(5Cn).Pz3(S00) 

EQUI vaLEnCE (K ( 1 ) iK | ) . ( K ( 2) ,K2 ) . f K ( 3 ) »K3 ) . ( THETa M ) .THETAA) * ( TH tTA( 
12) » THE TAB) , (THETA ( 3 ) .ThETAC) 

REAL M » L , K 

REAL MA#nUiMU»MUP£F 

REAL MAGOIF 

namelist/ input/g»x,k» theta. ub.v, a. omega. sigma. n. mur ef»temref,t, 

*TlIM,DELTaT.TEmP 0,PO»T£MLAP.R.X3O 
P I *3 . 1 4 J b V 2 6 5 4 
READ ( 5 , 7uU2 ) NnN 

READ(b»7&0l) BCDX.rcDyi »BCDy2,BCUy3,FLDYI*FL0y?*FLl>Y 3 
READlS, 100) G.x .K ,th£Ta ,UB.V 

REAOIB.lCu) A, OMEGA, SIGMA. D. MURE F.TEMReF.T.TLim.DelTAT. 
*TEMP0*P0,TEmLAP,R.X3u 
igc Format ( 8Eio. 3) 

700 1 format ( 1 2A6 , $ x > 

7002 FORMAT ( 12) 

DO 20 1*1.3 
FB ASS ( I ) - o. 

20 V ( I ) * UB ( I ) 

DGTOKDep 1/180, 

TmETAI sTHET A A*I)GT(jWD 
ThETA2*ThETAB*dGT0RD 
ThETA3*THETAC*0GT0R0 

I F I ThET A l ,EQ ,c. ANU . THET a 2 .£0 ,r. ) go TO 1 

C*-A*KI*ISIN(ThETai ) -Cos { THET A 1 )^TaN(THE7A2))/k3/<SIN1Th£Ta3) 

* -C0SITHETA3)«TAN(thETa2) ) 

B»-A*Kl*( S IN(THETAj)-CoS(THETAl)*TAN(ThETA3))/K2/(SIN(THETA2) 

• -COS l THETA2 ) *TAN ( THETA3 ) ) 
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GO TO 2 

1 6*-A«K t*.B/K2 
C*-A«K 1 • • S/K3 

2 CONTINUE 
*RlTt(6,9OO0) a • 0 » C 

9000 FORMAT ( /// IfcX * IP3E I 2 * 9 > 

UCON * US ( 1 ) /A 

00 1001 L l * 1 .NmN 

READ l 5 , 101 ) XD ( 2 ) *X0 (5 > .BCDX < 1 1 1 *RCnX ( 1 2 ) ,T *T|, tM tOELTAT, IDT * IDS 
101 F0RMAT12E10.3. lOX.2A6,8X.3FS*n,21&> 

TRF * T 
TRLIM « TtlM 
K ( 1 ) « XD ( 2 ) /D 
K ( 2 ) ■ K C 1 ) 

K ( 3 ) * K ( 1 ) 

OMEGA * XO < 5 ) *i /O 
ETA « 1 • Ufc“3G 
DO 10 1*1 .3 

Ml) * AMaXI t K ( 1 ) .ETA) 

10 CONTINUE 

call geomet 
call propty < x i 
6R ■ -G<3) 

GAM * GR/R/9 • 1 B9E+7 
XO ( 1 ) « < 1 2 • • N U ) / ( A*D > 

XD(3) * •66AA6a7*GR*D*0*(GAM-TEmLAP)/<A*A*TEMPo) 

XD<9) * UCON 

WRITE(6,90lS) rtCDx ( 1 1 1 c BCOX( 12) . ( *0 ( L2 ) *L2 a 1 • b ) 

WR I TE ( 6 , I nPUT ) 

90 1 B FORMAT ( 1H1 |25X .« TEST CaSE* 2A6//* DIMENSIONLESS GROUPS*//* N7«* lPtl 
12.9/* N8* * E 1 2 • S/ • N9«*El2*9/* NlO**£l2.0/* Nil**El2.9) 

CALL 1DENT < 8 * 1 ) 

DO 1002 LS * 1.IDS 
1 D ■ 0 

1 D X * 0 

b CALL 1NTEG IT »V .X .DFLTAT . IDT ) 

DO 200 1*1*3 
Full) * U ( I > / A 
F V ( I ) * V ( I ) / A 
V DU ( I ) * V 1 1 ) /Of I ) 

200 CONTINUE ..r 

AK T ■ A*K < l ) • T 
DZ * X ( 3 ) “ X 3n 
1 0 * ID ♦ 1 
Ii)X » IDX ♦ I 
PX(IDX) « X(l) 

PY ( IDX ) * X ( 2 ) 

PZ 1 I OX ) * DZ 
PYl(lOX) m FU(1) 

PY2 ( IDX )■ FU ( 7 ) 

PY3( IDX) « FU(3) 

PZ 1 < IDX ) . FV I l ) 

P Z 2 ( IDX) - FV(Z) 

P 2 3 < IDX) * FV(3) 

ARA C ID » k ) * T 

AR6 (10*1) * T 
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WRC ! ID 

, 1 ) * 

T 

WR A ! ID 

* 2 ) » 

AKT 

W R 6 ! I 0 

,2 ) «= 

AKT 

WRC ! ID 

,2 i * 

AKT 

PTt IDX 

) ■ AKT 

WR B ! ID 

• 3 ) » 

X ! 1 ) 

WRB ! ID 

, H ) a 

X !2) 

WRB ! I D 

*5) « 

DZ 

WRB l ID 

» 6 ) « 

RE 

DO 300 

I • 1 » 3 


WRA ( ID 

• I *2 ) 

* FU< I ) 

WRA! ID 

• I+S) 

■ fv ( n 

WRA { ID 

. 1+8 ) 

■ vDU | I ) 

WRC ! ID 

• 1+2) 

« FDRAG ! I ) 

W R C t I D 

• 1+5) 

* F A P 1 ( I ) 

WRC ( I D 

• 1+8) 

* F A P 2 ( I ) 

WRC t ID 

.1 + 11) 

• FBUO Y ! I ) 

WRC ! ID 

• I + 1 H) 

* FBASS t I ) 

CONT I NUE 


IF ! ID 

- so) 

3?S ,30? .302 


302 CONTINUE 

*R 1 TC (6 ,9002 ) 

Wk 1 T£ ( 6 , 9303 ) ( ( 'Vr a ( I . J > . J* I . 1 1 ) ♦ I * t . 50 ) 

WR 1 TE ! 6 ,9QQ6 ) 

WR I TE ! 6 , 9009 ) ( ( WRB ( I , J) , J«1 ,6) , 1*| 50) 

AH I TE ( 6 , 9006 ) 

WR I TE ! 6 ,9007 ) C ( WR C ( I . J ) . J- 1 . I 7 > , I « i , 50 1 

9002 FORMAT ! * I TIME ScaLED TIME SCALED aTmOsPhfRIC WIND 

1SCALE0 BALLOON VELOCITIES VELOCITY RATIOS*//) 

9003 FORMAT! IX.QPF8.2. iPlOEl 1 .3) 

9006 FORMAT ( * 1 TIME SCALED TIME X Y Dz RE 

1 *//) 

900V FORMAT! 1X.0PF3.?, ipSEl 1 ,3) 

9006 FORMAT!*! TIME SCALED TIME FDRAG/FrUOy 

I FAP1/FBASS FAp?»//) 

9 00 7 format! ix»opfb.2»ipioei i • 3 / 2 o x * 1 p we 1 1 • 3 ) 

10 « 0 

30s continue 

IF ( T • L T ,Tt 1 m ) GO TO 5 
XL ■ A*K ! 1 ) «TRF 
XR » A*K ! 1 ) *TL I M 
YB = b • 8 
Y T * 1 I . 2 

call ouik3L < - 1 ,xl ,xr . yb , yt • ih + .bcox .bcdy i ,-l dx ,pt ,py 1 > 

CALL WUIK3L!C,xL.XR.YB,YT,1H0.BCDx,qC0Y) .-IDX,PT*PZ1) 

YB • -1 .5 
YT ■ 1 .5 

Call wUlK3L!-l.XL,xR«Ya.YT,lH+,ftCDX t BCDY2,-lDx.PT # PY2) 

CALL 0UIK3L!0,xL»XR,YB,vT, 1 HO, DC OX, BCDY2, -IDX.pT, PZ2) 

CALL 0UIx3L!-1.XL,xR,YB,YT.IH+,BCdX,BCDY3,-I0x.PT,PY3) 

CALL 0UIK3L(O.xL.XR.YB.YT,1HO,6COx,BCDY3,-IDX,pT»PZ3) 
call ouik3V<-i . 1 ho. bcdx.fldy 2 . -iox.pt ,py> 

CALL 0UIK3V!-1,IH0.RCDX,FL0Y3.-IUX.pT,PZ) 

TRF * T 

XL I M » TL I M + TRL I M 
1 C0 2 CONTINUE 
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CALL E ND JOB 
T « O* 

X( 1 ) « o* 

X ( 2 ) - 0. 

X ( 3 ) -X30 
v ( i ) * ubi n 
V I 2 ) - ua<2) 
V (3 i * UB < 3 ) 
1001 CONTINUE 
STOP 
End 


SUBROUTINE ACCEL < T .V I #X I tOVDT .DEUT aT *S > 

COMMON/ CONS /PI .D.SIGMA.aP.M, VOL* 

•MuREF.TEMkEF, TEMPOtTEMLAP ,X30.R .GSRTL *GI 31 *P0*X33 

COMMON/ TEKMS/FoRAGI 3) iFAPK3) .FAP21 3) .FBUOT 13) . FB A SS(3) * 

• MA»NU.MU*RH0»Tt'MP*X(3) .FACTOR 
COMMON /VEL0C/U(3>«V(3>*UHAG.VNAG*MAG0IF 

real M.L.K 

REAL M A ,Nu • MU , MUREF 
REAL M AGO I F 
REAL MEFF 

DIMENSION VI (3) .Xl ( 3) .OV0T13) 

CALL PROP T Y ( X I ) 

CALL BUOYNT 

CALL APARNT ( T ,X I * V I ) 

VMAG— SORT (VI ( 1 I • * 2 ♦ V l (2)**2+Vl (31**?) 

UMAG-SORTiUl 1 ) **2«-u ( 2 ) **2*U ( 3 ) • *2 ) 

SUM«0. 

00 S I « 1 . 3 

5 SUM-SUM*. C V I < I ) -U < I > ) **2 

MAGO 1F-SQRT ( SUM ) 

CALL OR AG < V I I 

B-F ACTOR #SQR T (0ELTaT*S » 

MEFF-M+MA 
Do 1 U 1-1,3 

DV0T1 1 )-(FORAG( I >*FAP! ( I )*.FAP2i I )+FnUOTf I I^FBacS I 11 )/(MeFF*-B) 

io continue 

RETURN 

EnQ 
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SUdROuTlNE APARNT C T ,X I . V I) 

DIMENSION X I ( 3 ) • V I ( 3 ) 

c This routine calculate forces rfsulting from apparent mass and fluid 
COMMON/CONS/P I .DiSIGMA. AP,M, VOL, 

♦MUREF * TEMREF » TEMPO • TEML AP » X3f» , R ,gSRTL » G < 3 ) ,PO » X33 

COMMON/TERMS/FDR A G I 3) • F AP I C 3 ) . F A P2 t 3 ) .FBUOYC3) . F&ASS < 3 ) , 

• MA.NU, MU, RHO, TEMP, x(3). FACTOR 
COMMON / VEL0C/UI3) *V<3) ,UM A G,V MAG, MaGOIF 

COMMTON/DERV/PDuOX I 3,3) ,P0U0T(31 ^ 

COMMON/T JME/N . . - 

DIMENSION DUDTFI3) 

REAL M , L , < 

REAL MA ,NU ,mu ,muRef 
REAL MAPMRS 
realmagoif 

MAPMNS*MA«>M*RHu/S 1 GMA 

X I 3 * X M3)- X 30 

call potfluit.xi i i > ,xic 2 ) ,xi3) 

DO 10 1*1,3 
U D UD A *D V A L t U ,0. • I ) 

DuOTF ( I )*PDUDT I I ) ♦UDUDX 
FAPM I )*mapmrs*dudtfi I ) 

F A P2 ( 1 ) ■ MA*DVAL(vI , U , I ) 

ID CONTINUE 
RETURN 
End 


subroutine buoynt 

COMMON/CONS/PI, D, sigma. AP,M. VOL, 

♦ MUREF.TEMREF, TEMP0,TEMLAP,X3O.R,GSRTLtG(3) ,P0,X33 

COMMON/TtRMS/FDRAG(3),FAPM3),FAP2(3),F0UOY(3l,Ft»ASS(3), 
•MA,NU,MU,RHQ»TeMP.X<3>. FACTOR 
C 

real M,L,K 
real ma,nu,mu 

REAL MROSJG 

MR OSIG«M«(1,-RhO /SIGMA) 

DO 10 1*1.3 

10 FBUOT ( I ) »G II ) *MR0SIG 
RE TURN 
end 
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subroutine draci(vt) 

COMMON /C0NS/PI,D.SIGMA,AP,M, VOL. 

•MUREF . TEMREF , T£MP0.TEMLAP,x3O.R,GSRTL.e,l3).P0.X33 

COMMON/ T ERMS/FuR A G (3).FAPW3 >.FaP2<3>»FBUOY(3).FBaSS13), 

• M A * NO • MU . RHO . TFMP , X ( 3 ) .FACTOR »RE 
COMMON/ VElOC/U 1 3 > .V ( 3 ) .UMAG . VMAG .MAGO IF 
DIMENSION RT ABL ( 8 ) .REN ( fl ) 

DIMENSION V I ( 3 ) 

REAL M.L.K 

real ma »nu »mu .muref 

REAL MAGOIF 

Data RTAaL/7.74066.5.29a3.3.13b49,l,5,.n9S3|.-..693l47» 

* - 1. 23787.-. 544727/. nREn/8/, EPS/. 0J/.REN/-4. 6 Og,^, - 2. 3026,0.. 
•2.3026»4.6US2.6.9n7R,9.2I03» 1 1.5129/ 

RE* MAGDIF*r>/Nu 
RE Z » RE ♦ 1.00E-3O 
CO * 24./REZ ♦ 0*50 
VU-VKl) 

V I 2-V l ( 2 | 

V I3-VI (3) 

HALFRO*. 5*RH0 

coapmo»co*ap*magdif 
HACO*HALFRO*CDaPMq 
00 10 1*1.3 

FORA<j(I)*HACD*UMI»-VI(T>» 

10 CuNTlNUE 
RETURN 

end 


Function dvalia.z.m 

COMMON /DERV/PDU OX 13,3) .PDU0TC3) 
DIMENSION X(3) ,2(3) 

o v al*o • 

DO 10 J«1 ,3 

1C DyAL*DVAL6(M(J)-ZtJ))«PDUDX(I.j) 
RETURN 
EnO 
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subroutine geomet 

C This ROUTINE Calculates time independent GEOMETRIC values 
COMHGN/COnS/PI .D.SIGMA.aP.M, VOL, 

•MUREF .TEMkEF • T EMPO • TEML AP • X3n • R ., GSR TL • & ( 3 > ., PO • X33 

REAL M.LiK. MUREF 
VQL*P I *D**3/6 . 

H»VOL«S1GmA 
AP*PI *D**2/4,3 
IF(TLMLAP) 10 , 20 * 10 
10 CONTINUE 

GSRTL*ABS(G(3) ) / ( R «TEML AP ) /H • I8<4E + 7 
GO Td 30 
20 CONTINUE 

GSRTL*ABS(G(3)»/(K*TEMP0 ) / ** . I bME + 7 
30 CONTINUE 
RETURN 
ENU 


SUBROUTINE INTEG(T,V»X.r>ELTAT,lDTl 
0 I MENS ION X ( 3 ) . V ( 3 ) » D V 1 (3) ,DV2(3) , 0 V H ( 3 ) * D V 5 ( 3 ) , 

• VI (3) • X I ( 3 ) , DVoT ( 3 ) «DV 3 ( 3 ) 

DIMENSION DXI(3)iDX2(3> ,0X3(3), DXH(3> • 0X5(3) *£<3> ,EX(3) 
C0 MMUN/T I ME/N 
DATA ACCUR/1 */ 

DO 2uQ0 L 1 * I , I OT 

3 CALL ACCEL(T,V.X,DvOT,DFLTAT,1*) 

M DO S 1*1 ,3 

DVl ( I )«DELTAT*OVD t (1 >/3. 

0X1(1 )*D£lTAT*V ( I ) / 3 • 

V I ( I )*V ( I ) *0V 1(1) 

5 X I ( I ) *X ( I ) ♦DX l ( I ) 

CALL ACCEL(T*0ELTAT/3..vI,XI,DVDT. DEL TAT, ,3333. 3 33333) 

DO 6 1*1,3 

DV2( I )*DEL TAT •OVDT( I )/3. 

0X2 < I )*D£LTAT*vI ( I >/3. 

»iu)«v(i)**5*ovim**5*ov2m 

6 XI (I )*X(I ) + , 5 *l> X 1 l I )+.S*Dx2(I ) 

call accel(t+oelTat/3..vi,xi,dvbt*deltat,. 3333333333 ) 

DO 7 1 * 1 ,3 

D y 3 ( 1 )*DElTAT*OVDt ( I )/3. 

0x3 ( I ) *D£l T A T* v I ( M/3. 

VI ( I ) * V ( 1 )+0,375*0vl ( I ) + l«125»Uv3( I ) 

7 XI ( I )«X( I )*0,375*Ox 1 ( I ) ♦ 1 • 1 2 5 • D X 3 ( I > 

CALL ACCEL(T^DELTAT/2, ,VI ,XI , 0 V D T , DEL T A T • • 5 > 

DO 8 1*1,3 

D V <4 ( 1 )»0ELTAT*OVDT ( I ) /3 • 

DxR( I )*DE lTAT*v1 ( I )/3. 

VI ( I ) * V { I ) + l,5*DVl ( I )~M.5*DV3(i ) •6.0*0 VRi I ? 
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8 XI { I )*X( I I*.! «5*DXl < I )-M.5*0X3(I >«-6«r»*DX4(I ) 

Call a c c e l ( t ♦ d eltat*vi , x i , o v d r * ot ltaT#i«) 

TEST *0*0 

Do 9 1*1 .3 

OVS< 1 )*DElTAT*OVDT I I )/3. 

0X5(1 )*DElTAT*VI C I » / 3 • 

£( I ) *C • 2 *0 V 1 ( I J - • 9 * D V 3 C T ) ♦ • B * D V H ( I )-• 1 *‘D V 5 ( 1 ) 
Ex(I)*0.2*Dxin)-.9*DX3(I)^*8*DXH(i)-.l*DXS(I> 

9 TEST*AHAX1 (TEST.AB<;(£(I ) ) tABS (tX ( 1 ) ) ) 
IFITEST.LT*ACCUR> GO To 10 
WRIT£(6,l 000 ) 

100C FORMAT ( 10X *38H#**ACCURaCY TEST FAJLpO IN INTEGRATION) 

10 DO 11 1*1 ,3 

V(I}«V(I)+«5*DvM(I>+2.*DV‘HI)+.5*0V&(I> 

11 x(i)*x(i)4-.5*oxi<i)+2**ox < nn^*s*oxs(n 
T»T ♦DELT A T 

2000 CONTINUE 
RETURN 
END 


REAL F UNCTION MR E F ( T E MP , MUR E F , T E rt R E F ) 
REAL NUREF 

A»MUK£.F«(TEMP/TEMKEF>**1»S 

8*1 .505*TEMREF 

C*T E MP ♦ ■ 505 ♦TEMREF 

MHEF*A*B/C 

RE TURN 

End 


SUBROUTINE POTFLU ( t . x , Y ,Z ) 

COMMON/ INTL/K 1 .K2 .K3 .THETA 1 *THETA2 . THETA3 *0ME(5A ,Utt , va , WB , A »B ,C 

COMMON/OERV/DUDX »DV0X . DrtDX ,DUDY ,UVDY *DwDY »OUOz .pVOE ,OflUZ»QUOT > 
♦ OVOT »OWOT 

common/veloc/u.v»«.dumy(6j 
Real L l ,L2 tL3 »« 1 *K3 

XI*K1*X+K2*Y*K3*Z 
X1M0MT*XI-0MEGa«T 
COS l*CGS ( X IM0MT*THFTA n 
C0S2*C0S < X I M0MT + THFTA2 ) 

C0S3*C0S ( X I M0MT + THETA3 > 

U*U8-*-A*S 1 N ( X I MOMT + THETA 1 ) 

V*V8 + 8*S 1 N ( X I M0MT + THETA2 ) 

W**B + C*5 I N ( X ! M0MT + THETA3 1 
DU0X*A*K I *COS 1 
0U0Y*A«K2«C0S 1 
DuOZ*A*K3*C05 1 
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D VP X * B * K 1 *C0S2 
DVD Y*B*K2»C0S2 
DVUZ*B*K3*C0S2 
D#OX«C*K 1 • C 0 S 3 
DwO Y *C*K 2 *COS3 
D*iOZ*C*K3*COS3 
puDT»-A IOMEGA *COS 1 
pyi)T»-d*OMEGA*COS 2 
PWOT*-C»OMEGA*COS3 
2d CONTINUE , 

RETURN 

end 


SUBROUTINE PROPTY(XI) 

c this routine calculates temperature and density from 

C TABULATED FUNCTIONS OF ALTITUDE (X<j) OR ALT), DYNAMIC 

C VISCOSITY , MU . AND KINEMATIC V | S C Q S 1 T Y , nU , A R E CALCULATED. 

C PRESSURE ■ PO * DYnES/cM**2 

C DENSITY » RHO * GMS/C«*3 

DIMENSION X I ( 3 > 

COMMUN/CONS/PI,DiSIGMA,aP,M, VOL* 

•MUREFtTEMNEF. TEMPO .TEMLAP ,X3P.R.rtSRTLi6<3>.PO,X33 

C0MMDN/TERMS/FURAG(3),FAP1(3),FAP2(3) ,FRU0Y(3| ,FBaSS<3) , 

*MA »NU ,MU .RHO . TEMP , X I 3 > ,F ACTOR 

real m*l«k 

REAL MA ,NU »MU , mUREF 
REAL MREF 

TEMP*T£MPO^TEMl AP*(Xl ! 3 ) -X3H ) 

MU«-iK£F(TEMP,MuREF.TEMREF) 

1000 FORMAT! IOX .5HTEMP*,ElO.<4) 

2ooc format ( tox t3HMu- •£ i o*h > 

iFtTEMLAP) 1C, 20. in 

10 continue 

P*PC* (TEMP /TEMPO )**GSRTL 
GO TO 30 
20 CONTINUE 

P ■ PO*EXP I -GSRTL* 1 XI (3>-X3P) ) 

30 CONTINUE 

RH0«P/!R*T£MP ) / 4 • 1 fl 4E + 7 

NU*MU/RHO 

MA»V0L*Rh0*.5 

l do format (// 1 ox • i bhpropty variables) 

2QQ FORMAT ( ! OX • 5HT* MP«.E 1C. *f / lox, 3H X i«,,El 0 .<l/ 10 X,4HA30« it 10* S/ 
•10X,3 hMU*,E 10.4/1CX ,2HP«t£10. S/lUx ,&HGSRTL*,Ein.^/ 
*lOX,3HNU«,ElC.4/IOX.3HMA»,E10* i l/lOX,2HM*t£l0,4/l0X,4HV0L a ,tl0.H//) 
RETURN 
END 
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APPENDIX L 


Samples for Blocks #1, #2 fc #3 Outputs 

The sample outputs which follow correspond to Blocks #1 
#2 and #3 Outputs which are described in Section 4.2. These 
outputs are for the most part labeled and thus are self-explana 
tory. Where three components of a vector are printed out, the 
sequence always corresponds to the order , x^, x^ . 

L.l Block #\ Output (Sample) 

! 

1 •0000*02 1 #0000+02 I.C0PO+Q2 


L.2 Block #2 Output (Sample) 

TEST CASEtEST <1 . 1 > 

DIMENSIONLESS GROUPS 

N7« 3.5914-04 

NB* 8 • 168 1 —OH 
N9a 1 #7317-03 
N 1 0 * 1 *0000+01 

Nil* 8# 1681-03 
S INPUT 


G * 

• OOOOOOOOE + OG • 

• oooooftoot + nu « 

-• 98Q6650nE#.03 

X * 

• OOOOOOOOE + UO » 

•OOOOOoOOE+OO • 

• 1 lOOGC'Goe + 07 

K a 

•62831S38E-0& > 

i 6283 l • 

•62831536E-0S 

THETA * 

• OQOOOOOOE + OO > 

. 12OG0PC0E+O3 , 

#24000 C00£+03 

UH * 

• lGaOOQCCE + uH • 

•OOOOCOOOE+nO , 

•ooooonooE+oo 

V a 

• i ooooooae + o4 * 

•uOOOOnOOE+OO# 

•OOOOOOOOE+OO 

A a 

• 1 0000000E + 03 



OMEGA a 

• 6283 1 S38E-02 



SIGMA a 

.361494GQE-03 



0 a 

• 1 3000000E + 03 



MUREF « 

• *1 7080000E-03 



TEMKEF a 

•2731S000E+G3 



T a 

• aooaaoooE+uo 



T L I M a 

• lOOCOOOOE+OS 



OELTaT « 

• 2SOOOOOOE+U 1 



TEMPO » 

• 2 1 800000E + 03 



P0 * 

#2261 VOOOE+06 



TEMLAP ■ 

•OOOCOOOOL+OO 



H a 

♦ 6860000QE-0 1 



X 30 * 

• 1 1000000E + G7 
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L.3 Block #3 Output (Sample - next four pages) 
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.j&S. 


T|MC SCAlEO time 

SCALtO ATMOSPHERIC 

* I No 

SCALED 

BALLOU* VtlOClT lEa 

vlloci tt 

R A T 1 OS 


SO. 00 

3. 1*2-02 

1 .003+01 

6.C87-31 

9. *62-01 

1 .OOH+J1 

5.920-01 

5.837-01 

1.000+00 

9.733-01 

6* 100-01' 

100.00 

6.263*02 

1 *007+01 

6 » 37 5-0 1 

9 . 669-C 1 

i »007+u l 

6 « 3«9-0 1 

*. 3 19-01 

1*000+00 

1 .002+00 

*.*69-01 

1SO.OO 

9. *25-02 

1.0114-01 

6.6*5-01 

9.799-01 

1 .01 l+Ll 

6*661-01 

3, *76-01 

1.000+00 

1.002+00 

3.566-01 

200.00 

1.257-01 

1 *01 *+01 

6.93*t-C 1 

9.821-01 

1 .019+01 

6*919-01 

2.887-01 

1.000+00 

1.002+00 

2,900-Ql 

2S0.00 

1.571-01 

1 .018+01 

7.157-01 

9. 082-Pl 

1 .018+01 

7.J72-01 

2. **0-01 

1.000+00 

1.002+00 

2.069+01 

300.00 

1 • 805-0 1 

1.021+01 

7.005-31 

9.931-01 

1.021+01 

7. *)9-0l 

2.083-01 

1 .QOO+OO 

1 .002+00 

2*098-01 

350*1)0 

2.199-01 

1 .025+01 

7.6*9-01 

9.968-01 

1 .025+01 

7.662-01 

1 .789-01 

1.000+00 

1*002+00 

1.790+01 

*00.00 

2.513-01 

1 *029+0 1 

7 . 889-0 1 

9.991-01 

l .029+01 

7 . 9C*2-u 1 

1 .538-01 

1 .000+00 

1*002+00 

1.539-01 

*50.00 

2. 82/-01 

1.032+01 

8.129-0 1 

1 .030+00 

1 .032+o 1 

8. 1 36-01 

1.318-01 

1 .000+00 

1.002+00 

1 *314=01 

50o.no 

3.1*2-01 

1 *036+01 

8 .35*4-0 l 

9.993-01 

1 .036 + 01 

8.365-01 

1 • 119-01 

1.000+00 

1.001+00 

1*120-01 

550*00 

3 . *56-0 1 

1.0*0*01 

8. 576-01 

9.969-01 

1 .0*0 + Ul 

8 *557-0 1 

9.366-02 

1.000+00 

1.001+00 

9.396+02 

Aoo.no 

3.770-01 

1 .0***’01 

8.7 9 n-0 1 

9. 92S-Q1 

1 .o**+ui 

0 • 8n 1 *li l 

7.638-02 

1.000+00 

1 .ooi+uo 

7.696-02 

650*00 

* • Ob * -G 1 

1 .0*8 + 01 

8. 99*4-0 l 

9.861-01 

1 .0*8+01 

9 • On 3-0 1 

5.969-02 

1.000+00 

1*001+00 

6.050+02 

7oo.no 

*.398-01 

l .U52 + 01 

9. 1 85-0 1 

9.779-01 

1.052+01 

9 • l 9*-C 1 

*.329-02 

1 .OJO+OO 

1 .001+00 

*.*29-02 

75o.no 

*.712-01 

1.056+01 

9.362-01 

9.663-01 

1 .056 + 01 

9.369-01 

2.693-02 

1.000+00 

t.ooi+oo 

2.787+02 

000.00 

5.027-01 

1 .060 + 0 1 

9.52I-Q1 

9.527-31 

1 ,060 + Jl 

9.520-01 

1 *0**-Q2 

1*000+00 

1*001+00 

1*096-02 

aSQ.no 

5.3*1-01 

1 • 06* + 0 1 

9,661-01 

9.365-01 

1 • 06 * +0 1 

9*667-01 

-6.300-03 

1*000+00 

1*001+00 

-6*728-03 

903.00 

5.655-01 

1.068+ut 

9.70n-Ol 

9. 175-01 

1.068+01 

9 . 7R*“Q 1 

-2.336-02 

1.000+00 

i.ooo+oo 

-2*5*6-02 

950.00 

5.969-01 

1.071+01 

9.879-01 

8.957-01 

1.071+01 

9.070-01 

-*.070-02 

1 .QOO + OO 

1.000+00 

-*•6*0-02 

i ooo.no 

6.2B3-01 

1.075+01 

9.999-01 

8.710-01 

l .075+01 

9.9*46-01 

-5.828-02 

1 * 000+00 

1 . 000+00 

-6*691-02 

ioso.no 

6.597-01 

1.Q78+01 

9.906-3 l 

8.936-01 

l .070 + 01 

9.9*7-01 

-7.599-02 

1 . 000+00 

i.ooo+oo 

-9.007-02 

iloo.no 

6.911-01 

1 .08 1 + 01 

1 . OOp + C? 

0.135-01 

1.002+01 

1 .QnO + UO 

-9,367-02 

1 .000 + 00 

1 . 000+00 

-1.151-01 

1 150.00 

7.226-01 

1.085+01 

9.985-31 

7.800-01 

1 .085 + 01 

9 * VH 0 - Q 1 

- 1 . 111-01 

1 * 000+00 

9.999-01 

-1*023-01 

1200.00 

7.5*0-01 

1 *087+01 

9.990-31 

7.055-01 

1 ,087+ul 

9.937-0 l 

-1.281-Q1 

l *000 + 00 

9,998-01 

-1*719-01 

1 250.00 

7 • 85H-0 1 

1*090+01 

9.86*4-01 

7.080-01 

1.090+01 

9.061-01 

- l .***-01 

1 . 000+00 

9,997-01 

-2*0*0-01 

1300.00 

8. 168-01 

1 .092+01 

9.76*3-31 

6.603*01 

1 .092+01 

9.7S6-J1 

-1.597-01 

1 .000 + 00 

9.990-01 

-2*389-01 

I350.no 

8 . *82-01 

1.094+01 

9.627-0 1 

6.268-01 

1 *09 * + 0 1 

9.622-01 

-1.736-01 

1*000+00 

9.99H-01 

-2*770-01 

i*oo.no 

8 . 796-0 1 

1 .096 + 01 

9.967-01. 

5.837-01 

1 .096+01 

9. *6Q“01 

-l .059-01 

1 .000 + 00 

9.993-01 

-3*185-01 

1*50.00 

9.1 11-01 

1.097+01 

9. 20*1-01' 

6.392-01 

1 .097+ol 

9.272-01 

-1.9*3-01 

1*000+00 

9.992-01 

-3*600-01 

1500.00 

9. *25-01 

1 .090 + 01 

9.069*31 

*.937-01 

1 .098 + 01 

9.06C-0 1 

-2.0*3-01 

1 * 000+00 

9,990-01 

-*.139-01 

1 S5n.no 

9.739-01 

1 .099 + 01 

8 . 839-01 

*.*72-01 

1 .099 + 01 

0*825-01 

-2.097-01 

1*000+00 

9.989-Ul 

-0*690-01 

1600.00 

1. 005*00 

1 .100*01 

8.579-31 

*.002-01 

1 . 1 00+01 

0.569-U1 

-2. 122-01 

1 * 000+00 

9,988-01 

-5*303-01 

1650.00 

1.0374-00 

1.100*01 

8.3CS-CI 

3.520-01 

1 * 1 00 + 01 

0.29H-U1 

-2.115-01 

1 * 000+00 

9.986-01 

-5*996-01 

1700*00 

1.0684-00 

1 *100*01 

8 . Cl 9-01 

3.052-01 

1 . 1 00+01 

8.001-01 

-2.072-01 

1 . 000+00 

9 . 98*-G 1 

-6*789-01 

1750.00 

1 .1004-00 

1 .100*01 

7.707-3 1 

2.575-01 

1 . 100 +JI 

7 * 69*-Q 1 

-1.990-01 

1*000+00 

9,983-01 

-7.726-01 

1 000.00 

1 .131 4-00 

1.099*01 

7.387-0 1 

2. 100-01 

1.099+01 

7.372-01 

-1 .860-01 

1 .000 + 00 

9.981-01 

-8*873-01. 

i S5o.no 

. 1#162>00 

1.099+01 

7.059-01 

1 .628-01 

1.099+01 

7.039-01 

- 1.689-01 

1 . 000+00 

9,979-01 

- 1 * 038*00 

1900.00 

1. 19*+Q0 

1.090+01 

6.710-01 

1 .150-01 

l *098+01 

6*690-0 1 

-1**58-01 

1*000+00 

9.977-01 

-1*259+00 

1950.00 

1.225+00 

1 .097+01 

6.355-0 1 

6.9 1 6-02 

1 .097+0 1 

6.339-0 1 

-1.157-01 

1*000+00 

9.970-01 

-1*673+00 

2000.00 

1.257+00 

1.095+01 

5.990-31 

2.279-02 

t .095+01 

5.973*01 

-7 . 555-02 

9.999-01 

9.972-01 

-3*316+00 

2050*00 

1 .288 + 00 

1 »09* + 0 j 

5.61J-31 

-2.35*-02 

1 . 09* + J 1 

5.595-01 

-l .591-02 

9.999-01 

9,y6V-01 

6*760-01 

2100.00 

1.319+00 

1 .092 + 01 

5.222-31 

-6.990-02 

1.092+01 

5.273-01 

1 .837-02 

1.000+00 

9,996-01 

-2.625-01 

2150.00 

1.351+00 

1 .090*01 

*.830-31 

-1 . 1*2-01 

1 *090+0 1 

*• 0 * 8-01 

-1.921-02 

1*000+00 

1 .002+00 

1*682-01 

2200.00 

1 *382+00 

1 .088+01 

*. * 68-01 

-1 .556-01 

1 .088 + 01 

* . * 48*-0 1 

-3.5/7-02 

1 .000+00 

1 .Q0*+00 

2.298-01 

2250.00 

l.*l*+00 

1 *086+01 

*.109-01 

-1 .9*8-01 

1.086+01 

H* 127-01 

-3.0/1-02 

l. ooo+oo 

1 .005+00 

1*987-01 

2300.00 

1 « **5 + 00 

1 .080+01 

3.759-31 

-2.319-01 

1 .085 + 01 

3.779-01 

-*.323-02 

1.000+00 

1.005+00 

1 *860-01 

2350* nO 

1 .*77 + 00 

1 .082+01 

3. *22-01 

-2.669-01 

1.003+ot 

3 • <4 * 2-u 1 

-*,8u*-02 

1 .000+00 

1.006+00 

1*799-01 

2*00.00 

1.508+00 

1 .080*01 

3.097-C 1 

-2.999-01 

1.081+01 

3. n e-oi 

-5,213-02 

1.000+00 

1.007+00 

1*738-01 

2*50.00 

1.539+00 

1.078+01 

2.706-01 

-3.308-01 

1 . 0794QJ 

2 • a n7 -0 t 

-5.SH7-02 

1.000+00 

1.007+00 

1*677-01 

25oo.no 

1.571+00 

1 .076 + 01 

2. *87-01 

-3.598-01 

• 1 .077+01 

2.5n8*01 

-5.813-02 

1 .ooo+oo 

1.008+00 

1*615-01 



TIME SCALED TIME 

X 

Y 

DZ 

KE 

50 . 00 

3. 1 92-02 

5 . 009+09 

2.372+03 

3.093+03 

1 .290+09 

ioo.no 

6.263-02 

1 .009+05 

5,978+03 

5.580+03 

i . 770+09 

ISO. DO 

9 .925-02 

1*508+05 

6.738+03 

7.519+03 

2.071+09 

200.00 

1 .257-01 

2.015+05 

1 .213 + 39 

9.097+03 

2. 289+09 

250.00 

1.571-01 

2.523+05 

l • 5 6 a + 0 9 

l .092 + P9 

2,996+09 

3oo.no 

1 .<165-01 

3.0J3+05 

1 .930 + 09 

1 . 1 55 + 09 

2.575 + 1*9 

350.00 

2.199-01 

3.599+05 

2.307+09 

1.252*09 

2.680+09 

HOD. DO 

2.513-0 1 

9.058+05 

2.697+09 

1 . 335 + D9 

2.766+09 

*♦50.00 

2.627-01 

9.573+05 

3.096+09 

1 .906+09 

2.838+09 

500.00 

3. 192-01 

5 . 090*05 

3.510+09 

I . 967*09 

2.868+09 

550.00 

3.956-01 

5.609+05 

3.939+09 

1 .518+09 

2.997+09 

600. OO 

3.770-0 1 

6. 1 30+05 

9.369+09 

1 .560*09 

2.987+09 

650*00 

9.069-01 

6.65 3 + 05 

9.819+09 

I .599+09 

3.019+09 

700*00 

9.398-01 

7 . 178+05 

5.269+09 

1 .620+09 

3.093+09 

750.00 

9.712-01 

7 . 705+05 

5.733+C9 

1.637+09 

3.059+09 

800.00 

5,027-01 

8.239+05 

6.205+09 

I .697+09 

3.068+09 

850.00 

5.391-01 

8,7 65 + 05 

6.688+09 

1 .698+09 

3.070+09 

900.00 

5*655-01 

9.298+05 

7.1 72+09 

1 .690+09 

3,069 + u 9 

950.00 

5.969-01 

9.833+05 

7*663+09 

1 .625+09 

3.090*09 

1000*00 

6.283-01 

1 .037 + 06 

8.159+09 

1 . 600+09 

3.028*09 

1053.00 

6.597-01 

1 .091+06 

8.659+09 

1.566+09 

2 .998+09 

1 1 00.00 

6.911-01 

1 • 1 95+06 

9. 157+09 

1 .529+09 

2.960+09 

1150.00 

7 .22a — 01 

1 . 199+06 

9,657+09 

l .973+09 

2.912+09 

1200.00 

7 . 590-01 

1 .253+06 

l .01a+05 

1 .913+09 

2.855+09 

1253.00 

7.659-01 

1 .308+06 

1 * 065+05 

1 .395+09 

2.789+09 

1 300.00 

6.168-01 

1.362+06 

1 . 1 1 9+CS 

1 .269+09 

2*7 1 2+09 

1 350.00 

8,982-01 

1.9 1 7 + 06 

1 , 163 + 05 

1 .185+09 

2.625+09 

i <4oo. no 

6.796-01 

1 .972 + 06 

1 ,2 10 + 05 

1 .096 + 09 

2 .528+09 

1 950.00 

9. 1 1 1-01 

1 ,526 + 06 

1 , 257 + 05 

9.999+03 

2.91 9+09 

1500*00 

9.925-01 

1 *58 1+06 

l .303 + 05 

8.997+03 

2 .300+09 

1550*00 

9,739-01 

1 ,636 + 06 

1 .399 + 05 

7.962+03 

2 • 1 66+09 

1600.00 

1 • 005 + 00 

1,891 +06 

l . 39 i +05 

6.906+03 

2 .029+09 

1650.00 

1 . 037+00 

1 .796+06 

l .933+05 

S .895+03 

1 .668 + 09 

1700.00 

1 . 068+00 

1.801 +06 

l .979 + 05 

9.798+03 

1 .699 + 09 

1750*00 

1 .100 + 00 

1 .856+06 

l .513+05 

3 . 7 8 1 ♦ G 3 

1.516+09 

1800*00 

1.131+00 

1.911 +06 

l .SS 1 +05 

2.81 6+03 

1 * 3 19 + 09 

1650*00 

1 , 162+00 

1 . 9&6+06 

1.587+05 

1 .926*03 

1 . 1 05 + 09 

1 900.00 

1 . 199+00 

2.021+06 

1.621+05 

1 . 137 + 03 

8.727+03 

1 950.00 

i . 225+00 

2*076 + tf6 

1 .659 + 05 

9.806+0 2 

6. 1 72+03 

20oo.no 

1 .257+00 

2. 131+06 

1 ,685+05 

-2.281+CO 

3 ,286+03 

2050.00 

1 • 288+00 

2 . 1 8 5 + Q6 

1.719+05 

-2.929+02 

2.629+02 

2100.00 

l . 319+00 

2 .290+06 

1.791+05 

- 1 .885 + 02 

2.953+03 

2150.00 

1.351 +00 

2.295+06 

1 .766+05 

-1 .883 + 02 

3, 1 76+03 

2200.00 

1 . 382+00 

2.399+06 

1.789+05 

-3.372+02 

9.007+03 

22 So . 00 

1.919+00 

2,903+06 

1.611 +05 

-5.232+02 

5.321 +03 

2300.00 

1 .995+00 

2.958+06 

1 .830+05 

-7. 269 + 02 

6.3 13+03 

2350.00 

1 .977+00 

2,512+06 

1 .898*05 

-9.598+02 

7 .326+03 

2900.00 

1 . 508+00 

2.586+06 

1 .865+05 

-1 .205 + 03 

8.299+03 

2950*00 

1 ,539+00 

2 , 620+06 

1 .860+05 

-1 .979*03 

9.222+03 

2500.00 

1 .571+00 

2.679+06 

1 .893+05 

- 1 .758+03 

1.011 *09 
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TIME SCALED TIME 


FOHAG/FbUOr 


50.00 

3.192-02 

- 1 *057*0 1 

7 • 2 7 b ♦ G 1 

1 .665*03 



0*000 

o.coc 

- 1.972*03 

-ULo.au 

4*252 !-02 

-1. 350 + 01 

-9. 19b*CP 

3.9Q6+03 



0.000 

o.ocr 

-3.551+03 

ISO. 00 

9.925-02 

-1 .568*01 

- 1 . i e«4*s i 

9.674*03 



O.OUG 

0.000 

-9.779*03 

-200.00 

1.257*01 

-1.719*01 

-1 .262*31 

5.698+03 



0.000 

Q.con 

-5.773*03 

.250*00 

1.571-01 

-1 .828+01 

-1 .303*01 

6.599+03 



0.000 

0.000 

-6.608+03 

-200*00 

1.885-01 

-1.908+01 

•1 .318*01 

7.269+03 



O.OCO 

0.000 

-7.314*03 

-350.00 

2.199-01 

-1.965+01 

-1.313*3 1 

7.883+03 



O.OOC 

O.OCO 

-7.922+03 

. *<00.-00 

2.5 13-01 

-2 • Qu9*0 1 

-1 .291*01 

8.909+03 



0.000 

o.ocr 

-8.991*03 

.850.00 

2.827-01 

-2.027+01 

-1.254*01 

8.84 1 +03 



0.000 

O.COP 

-8.887*03 

-500.00 

3,192-01 

-2.037+01 

- 1 .207*01 

9.294*03 



0.000 

0 » OOP 

-9.247*03 

-sso+oo 

2.454-01 

-2.033+01 

-1 • 198*01 

9.571 *03 



O.OOC 

a.ccr 

-9.587*03 

-AO 0*00 

3*770-01 

-2.01 B+01 

-1 «07tt*C l 

9.890+03 



0.000 

O.CCO 

-9.852*03 

-450.00 

H. 084-01 

-1.991+01 

-9.986*00 

1.004+09 



0.000 

0.00c 

-1 .004*09 

_ 2 aix»oa 

9 . 398-01 

-1.952+01 

-9. 1 l 2 * C 0 

1.022+09 



0.000 

O.OCn 

-1 .022*09 

-250*00 

H. 712-01 

- 1.903 + 0 1 

-8.151 *00 

l .033*09 



0.000 

0.C0O 

-1 .033+09 

AO 0*00 

5.027*01 

*“ 1 .89 3 + 0 1 

-7 . 1 18*00 

1.090+09 



0.000 

O.OCO 

-1.039*09 

•So. 00 

-5-*2H L»o l 

-*1. 772+01 

-4.031*00 

1.091*09 



0.000 

0 . C 0 O 

-1 .090*09 

- 200*00 

5 .655*0 1 

*1.692+01 

-9.894+00 

1 .034+09 



0.000 

O.GOo 

-1.035+09 

— * 50.00 

5 .969*0 1 

-1*603+01 

-3 . 7 i 4+Oc 

1.027*09 



0.000 

O.COo 

-1 .025+09 

1 000.00 

5*253*01 

*1 .506 + 01 

-2 • 52o*0C 

1 .012*09 



O.OCO 

O.OCO 

-1.010*09 

1 05 O* op 

6*597*01 

— l . HQC+Cl 

-1.31 3*00 

9.915+03 



0.000 

0.000 

-9.889+03 

44 - 00*00 

5*911-01 

-1.289+01 

-1.125-01 

9.655+03 



o.aoo 

0 . 00 O 

-9.625*03 

4 - 150 * 00 - 

-7*226*01 

-4*172+01 

1 .064+00 

9.390+03 



0*000 

O.COO 

-9.305*03 

4200-. OO 

7*540*01 

-1*051+01 

2. 200+00 

8.970*03 



0.000 

0.000 

-•.932*03 

1250.00 

7.859*01 

-9.275+00 

3.279+00 

8.598+03 



0.000 

0.000 

-8.505*03 

4 >00.00 

5 * 465*04 

-^5.034 + 00 

9 .28o*00 

8. 077+03 



0.000 

0.000 

-8.029*03 

1 >50*00 

8.982-Ql 

-6*794*00 

5. 1 89+00 

7.558+03 


F*P1 / F m A S 5 


FAP2 


6.234*31 

9.951*01 

1 .824*01 

-5.021+00 

-3.986+00 

-1.970*00 

0 .000 

0.CO" 

0 .000 




6.993+01 

5*015*01 

1.679+01 

-6.858+00 

-5.297+00 

*1.768*00 

u .ooe 

O.OnO 

0.000 




4.726*01 

5.rm5*o t 

1 .506 + 01 

-8.005+00 

-6.017+00 

-1.793*00 

0.000 

O.onc 1 

0.000 



6.990*01 

5*072*01 

1.319+01 

-8.7 9 9 + 00 

-6.927+00 

-1.671*00 

0.000 

O.OnC 

O.OOO 




7.138+01 

5.044*01 

1.111+01 

-9.369+00 

-6. 699+00 

-1.958*00 

0.000 

0 . Cno 

U.UOL) 




7.320*01 

5.039+Cl 

B . '7a 4 + 00 

-9.799+00 

-6.737+00 

-1.175*00 

0.000 

O.OnO 

0.000 




7 .483 + 0 1 

9.977*01 

6.210+00 

-1.011*01 

-6.721+00 

-8.387-01 

u.OOQ 

O.OnO 

O.QOCl 




7.627+01 

9.691*01 

3 . 369+00 

-1.032+01 

-6.420+00 

-9.561-01 

0.000 

O.OnO 

O.OQU 




7.748+01 

9.773*01 

2.510-01 

- 1 .096 + 01 

-6.996+00 

-3.390-02 

0.000 

0 » OnO 

0.000 




7 .843 + 01 

9.473*01 

-3. 151+00 

-1.053+01 

-6.206+00 

4. 230-01 

0.000 

O.OnO 

0.000 




/ . 908+0 l 

9.434*01 

-6 . B 37 + 00 

- 1.053 + 01 

-5.906+00 

9.101-01 

U.OOQ 

O.OnO 

O.OOU 




7 # 94Q + U 1 

9. 213*01 

-1 . UBU+O 1 

-1 .096+01 

— 5 . 550+00 

1.423*00 

0.000 

O.OnO 

O.OOU 




7.935*01 

3.950*01 

-1*503+01 

-1 * 033 + 01 

-5. 199+00 

1.958*00 

0.000 

O.OnO 

0.000 




7 .889 + 01 

3.697*01 

-1 .950*01 

-1 .015*01 

-9.690+00 

2.508*00 

0.000 

0 * OnO 

0.000 




7.799+01 

3 • 3n9*0 l 

-2.919+01 

-9.899+00 

-9. 199+00 

3.070*00 

O.OOtl 

O.OnO 

O.OUO 




7.663+0 1 

2.9?2*Cl 

-2.909+01 

-9.599+00 

-3.65 V + 00 

3.636*00 

0.000 

O.OnO 

O.OOO 




7.478+01 

2 » 5n3*0 1 

-3.902+01 

-9.239+00 

-3. 091+00 

4.200*00 

O '. 00 0 

0 • OnC 

0.000 




7.299+Jl 

2 * 040*0 1 

-3.906+01 

-8,819+00 

-2.996+00 

4.756*00 

0 .noo 

o.onc 

0.000 




6.960+01 

1 .566*0 1 

-9.91U+01 

-8.359+00 

-1 . BBU+OO 

5.293*00 

0.000 

O.OnO 

0.000 




6.428+01 

1 *057*01 

-9 . Vq 7 + o l 

-7.892+00 

-1.251+00 

5.806*00 

0.000 

O.OnO 

0.000 




6. 250*01 

5*289*00 

-5.369+01 

-7 ■ 269+00 

-6. 146-01 

6.284*00 

0.000 

O.OnO 

0.000 




5.830+01 

-1 » 160“01 

-5.B9V+01 

-6.699+00 

1 • 333-02 

6.721*00 

0.000 

O.OnO 

u.OGO 




5.373+01 

-5*570*00 

-6 » 2B0 + 0 1 

-6*080+00 

6.303-01 

7.106*00 

0.000 

o.ono 

0.000 




9.889*01 

-1 » ln0*01 

-6.675+01 

-5.939+00 

1.225+00 

7.433*00 

0.000 

0 . OnO 

0.000 




9.371+01 

-1 *631*01 

-7.027+01 

-9.786+00 

1.786+00 

7.695*00 

0.000 

O.OnD 

O.OOU 




3.839*0! 

-2. 145*01 

-7.332+01 

-9.129+00 

2.306+00 

7.885*00 

o.ooo 

O.OnO 

O.OOU 




3.298+01 

-2*633*01 

-7.585+01 

-3.977+00 

2.775+00 

7.996*00 
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o.ooc 

o.ocr 

-7 .507*03 

O.onO 

D.ono 

0.000 




1400.0Q 

8.796-Ql 

-5. 588+00 

5.99C+0C 

6.998+03 

2,754+01 

-3.080*01 

-/, 783 + 01 

-2.840*00 

3.184+00 

• a 024*00 



0.000 

U.CCn 

-6,942+03 

0*000 

0*000 

0.000 




1450.00 

9.111-01 

”4.423+00 

6. 662+00 

6.401+03 

2.216+01 

-3,bn6*Gl 

-7,925+01 

-2.228+00 

3.526+00 

7*9-69*00- 



0*000 

O.CCO 

-6.341+03 

(1.000 

O.OnO 

U.OOu 




15DC.OQ 

9 » H 25*0 1 

-3.316+00 

7 .196 + 00 

5.774+03 

1 .689+01 

-3. 882*0 1 

-8.010+01 

-1.650+00 

3.Z9 l+.OO 

J-.4U2-H.+O0 



0.00 0 

o.oco 

-5.710*03 

0.000 

0 * ono 

0 * 000 




1550.00 

9. 739-Ql 

-2.289+00 

7 .575+00 

5.125+03 

1 « 182+01 

-4.21 l+Ct 

-8.Q3V+01 

-l. 115 + 00 

3.9/6+00 

7 *5.8-9*00 



O.OOC 

O.OOn 

-5.057+03 

0.000 

O.onO 

0.000 




1600.00 

l . 005*00 

”1.357+00 

7 .788 + 00 

4.462+03 

6.979+00 

-4.492*01 

-8.012+01 

-6. 327-Ql 

4 .-072+00 

7*24-4 *00- 



0.000 

0. OOP 

-4.390+03 

0.000 

o.ono 

0.000 




1650*00 

1 * 037+00 

-5.366-01 

7 ■ 82 1 + 30 

3.795+03 

2.434+00 

-4*722*01 

-7.932+01 

-2.101-01 

4.076+00 

6.84-8*00 



0.000 

O.OOn 

-3.719+03 

U.000 

0 * OnO 

0.000 




1700.00 

1 .068+00 

1.484-01 

7 • 6 7 p + QO 

3. 134+03 

- 1 .779 + 00 

-4.901*01 

-7 .802 + 01 

1.446-01 

3.984+00 

6*342+00 



0 . OuO 

O.OCO 

-3.055+03 

0.000 

O.ono 

0.000 




1 750.00 

1 • 100+00 

6.991-01 

7. 321+30 

2.493+03 

.5.627+00 

-5.0?9*0l 

-7 «62» + 0l 

4.240-01 

3.790+00 

5*747+00 



o.ooc 

O.OOn 

-2.409+03 

0.000 

O.onO 

0.000 




1800.00 

1 .131+00 

1 .090+00 

6.761+00 

1.884+03 

•9 . 086+00 

-5.ln6*0 t 

-7.407+01 

6.210-01 

3.490+00 

5*062+00 



0.000 

0 .000 

-1 .796+03 

0.000 

O.OnO 

0.000 




1850.00 

1 . 162+00 

1 .309+00 

5 . 9B2+0C 

1 .322+03 

-1.214+01 

-5. 135*01 

-7.148+01 

7.279-01 

3*079+00 

4*286+00 



O.ouu 

0.000 

-1 .229+03 

0.000 

O.onO 

0.000 




1900.00 

l .194 + 00 

1.336+00 

4. 957 + 00 

8.242+02 

-1.477+01 

-5. i | 6 * U 1 

-6.853+01 

7.361-01 

2.549+00 

3.415*00 



0.000 

O.OOn 

-7.263+02 

0.000 

O.onO 

0.000 




1950.00 

1 .225+00 

1,148+00 

3 .668+00 

4 . 128+02 

-1 .698+01 

-5.051+01 

-6.526+01 

6.338-01 

1.885+00 

2.435*00 



0.0 UC> 

Q.COO 

-3.070+02 

0.000 

O.OnC 

0.000 




2000.00 

1,257+00 

7 .q 98-C l 

2.039+CC 

1 .177+02 

-1*877+01 

-4.941*01 

-6, 169+01 

3.998-01 

1.053+QQ 

1.314*00 



0.000 

u . OOP 

1 .495+00 

U.oOO 

O.OnO 

0.000 




205o.no 

1 . 288 + 00 

7.Q58-02 

1 .938-01 

-8.500-01 

-2. Q13+01 

-4.787*01 

-5. 782+01 

-2.390-02 

-5.684-02 

-6*866-02 



0.000 

O.OOo 

1 .553+02 

0.000 

O.ono 

c » uou 




2100*00 

1.319+00 

2*154-02 

2.5C1-01 

-9.520+01 

•2. 105+01 

-4 »5h9*C 1 

-5.368+01 

-4.503-01 

-9,819-01 

-1*1 4.9+00 



0.000 

o.oce 

I .205+02 

0.000 

0 . 0r»0 

O.OQU 




2 150*00 

1.351+00 

— 6 * q7 8*0 1 

-1 . I 7 1+00 

-l . 100+02 

-2. 1 5 2 + u 1 

-4.363*01 

-4.953+01 

-5.447-01 

-1.104+00 

-1.254*00 



0.030 

o.ocr 

1.204+02 

O.OOC 

o.ono 

O.OOu 




2200.00 

1.382+00 

-1.202+00 

-2 . 332+CG 

-1 .745+02 

-2. 160+01 

-4 . l 20*01 

-4.549+01 

-7.496-01 

-1.430+00 

-1.579*00 



0.000 

O.COr 

2.156+02 

0.000 

- O.OnO 

0.000 




225o.oo 

1.414+03 

-1 .910+00 

-3.553 + 0 o 

-2.953+02 

-2. I 34 + 01 

-3.862*01 

-4,155+01 

“1 .047 + 00 

-1*895+00 

-2.019*00 



0.000 

G.OCn 

3.346+02 

O.OOO 

0 « One 

0.000 




2300.00 

1.445+00 

-2.561+00 

-4.533+CD 

-4.31 0+02 

-2.079+01 
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APPENDIX M 


Initial Testing of Balloon Program 

In general according to first-order theory, the equation 
for a balloon in a stationary atmosphere can be written as 




dv 


dt 


- + 2a v 3 + N 2 x£ = G 1 


+ F 


B 


(M-l ) 


In this equation, the vertical displacement from equilibrium is 


x 


< =x 3 


The mass of the balloon is m, and the apparent mass 


of the displaced fluid is m, 


In general, for a spherical balloon 


displaced a small distance from equilibrium, (m + m )/m=3/2. 

ci 

In Equation (M-l) a is the linear drag and is the non- 
linear contribution to the drag term based on the following approxi- 
mation for drag coefficient 


C D = 


24 

Re 


'm 


( M-2 ) 


where C^is a constant and Re is the Reynolds number ( Re = I v 3 1 ) 

In Equation (M-l) the term Fg represents the Basset force 
which arises in a transient flow at low Reynolds numbers. 

2 

The buoyancy term N has been reduced to its linear com- 
ponent. Perturbation techniques used to evaluate the nonlinear 
contributions show that the first-order term contributes a modifi- 
cation of about 1 cm in 10 m to the amplitude of the oscillation 
and a subharmonic term of similar amplitude. These nonlinear 
contributions should be observed using a small time step and careful 
analysis. The perturbation technique cannot be used to solve the 
linear oscillator in general, as, the solution does not converge 
due to secular terms which arise in solutions higher than first- 
order. However, the technique seems applicable here to determine, 
the order of magnitude of the nonlinear contribution since no 
resonance phenomena are expected in the unforced system. 
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In the first foul* tests respectively the following simplifications 
were made : _ 

(1) (m + m a )/m = 1, a - 0, Fg * 0, = 0 

(2) (m + m a )/m =1, a f 0, Fg =0, G ± = 0 

(3) (m + m a )/m - 3/2, a = 0 / Fg = 0, G 1 = 0 

(4) (m + m a )/m - 1 , a = 0, Fg f 0, G 1 = 0 

For test case the governing equation is the same as that for a 
parcel of air displaced from its equilibrium position in the absence 
drag forces, Basset forces, apparent mass forces and fluid accel- 
eration forces. With these simplifications the governing equation 
is 




(M-3 ) 


where the terms have their usual meaning. 


For the atmosphere, by means of the hydrostatic equation, 

0 - - :£H 3 - g - (M - 4) 

The atmospheric temperature is assumed to vary linearly with alti- 
tude according to the relation. 

T = T - y x' (M-5 ) 

o o 

The temperature of the parcel is also assumed to vary linearly with 
altitude according to the relation 


T = T 
P o 



x. 


(M-6) 


As shown by Hess [181] , based on the four preceding equations 
combined with the ideal gas law, by first-order theory the oscilla- 
tion of the air parcel can be predicted. The period of oscillation is 


2t r 


T “ 


/ g/T Q (y p -y) 


(M-7 ) 


Normally y is assumed to be equal to the adiabatic lapse rate, r , 
P s 

and the resulting relation is 

2 TT 


B - V - A/t 0 (?.“ 


(M-8 ) 


where the subscript B.V. refers to Brunt-Vaisala . 
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As pointed out by Angell and Pack [131] , if y is set equal 

P 

to the constant density lapse rate, r (instead of the adiabatic 
lapse rate), the motion of the air parcel is the same as that of a 
constant volume balloon with a period of oscillation of 


CVB / g/T 0 (r p - y) 

Now 

r p = -.0341 °K/m 


(M-9 ) 


(M-10 ) 


which is about six times larger than the lapse rate normally 
observed in the lower atmosphere. For the case of an isothermal 
atmosphere. 


X CVB 


/ ^ T o F P 


= 10.85 /T 

o 


(M-ll ) 


Equation (M-ll) closely agrees with 
Now for the test case #1, 
y = 0 

and 

T = 218°K 
o 

the time period is 

t Cvb = 160.20 sec. 


the relation given by Lally [58] . 
with an isothermal atmosphere, 

(M-12) 

(M-13 ) 

(M-14 ) 


The balloon displacement as a function of time, based on the computer 
program is presented in Table (M-l). Inspection of this table 
reveals a time period of 160 seconds, in close agreement with 
Equation (M-12). 

For the^ case #2, the solution to Equation (M-l) is 

= Ae -at cos At (M-15 ) 

where 

A - (N 2 - a 2 )* (M-16) 
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* Note: An error in the balloon density places the equilibrium 

position at x Q - Xort ~ -3m. 

«j 30 
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2 2 

In the case studied, it was found that a << N . 
period 


Consequently the 


x = 


\l 1 - T 2 /T 2 
* o ' a 


where 


- T 


T q = 27T/N 


(M-17 ) 
(M-18 ) 


T a = 2Tr/a 
N = \/ T ! 


(M-19) 

(M-20) 


Hence, no change in this period was observed in the numerical 
solution due to the dissipation term. However, the amplitude 
damping in the numerical solution was found to correspond almost 
exactly to the theoretical value. 


In the Case #3 test, it was observed that the period f 
given by the numerical solution corresponded closely to the first- 
order solution which required 


T 



(M-21) 


In the fourth test, the Basset force was non-zero, 


F 


B 



dv 3 /dt ’ 

^ t - t r 


-dt 1 


(M-22 ) 


No attempt was made to obtain a solution to the integral equation. 
The results of the test suggested that the Basset term behaves 
like a term, of the form, 


F fi = 2a v 3 - b 2 x^ 


leading to an equation of the form 

^ v 3 2 2 

^ + 2a v 3 + (« - b z ) Xg = 0 


(M-23 ) 


(M-24 ) 
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where the fundamental period was altered and the motion dampened. 
The oscillation was regular over two periods leading to the sup- 
position that 


a, b f f (x£) (M-25 ) 

The value of "a" was of the same order of magnitude as the coeffi- 

o 

cient of the integral. However, ”b ” was found to be an order of 
magnitude greater; 

b 2 ~ 10 a ~ 10" 2 (M-26 ) 

All of the preceding tests were conducted assuming an 
isothermal atmosphere with the natural period given by the relation 

T o = T o (M - 27) 

instead of the more general relation 

t = T / % ( M-28 ) 

° V Y o - Y 

One additional test case (#5) was run to ensure that 7=0 was not 
a special case. The period of this test case and the period obtain- 
ed from the analytic solution of the equation were again nearly 
identical . 
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APPENDIX N 


Results of Test Runs #6 - #11 


Based on the numerical solution produced by BALLOON, 
the variation of the vertical component of velocity (for the wind 
and the balloon) versus time is presented in Figures N-l through 
N~5 for test cases #6 - #li. 
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Figure N-2 . Vertical Component of Velocity Versus Time, Test Case #8 
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Figure N— 3 . Vertical Component of Velocity Versus Time, Test Case #9 





APPENDIX O 


Analysis of Numerical Results, Phase I (16 Origi nal Runs ) 


For each of the sixteen runs which were initially performed 
as part of Phase I a set of 6 position and velocity plots was 
generated. These plots are not included in the report because of 
the considerable bulk which they represent. They have been 
collected and bound in a separate document [l] for reference purposes. 

For each of the sixteen runs special attention was given 
to the length of the time interval over which the numerical inte- 
gration was carried out. The time interval used along with the 
time period of the associated flow field are presented in Figure 0-1. 

Figure 0-2 presents the linear period, the observed period 
and the observed phase lag in the balloon velocity for the 16 runs. 

The wind is given by: 

u^ = u^ + A sin (kjXj - cot + 0.^) (0-1) 

Based on the assumption 

x. = u.t (0-2) 

l l 


then the linear period is given by 
- 2tt 

T * i^i - °i 

since k. = k and u. = u . 5 ■ , and 

i l i li 


(0-3) 


c = 


0J_ 

k 


(0-4) 


is the phase velocity. The observed periods are taken directly 
from the data represented in Figures 5-4 through 5-6. 

In Figure 0-2 values to the right of the principal 
diagonal represent cases where u^ > c and to the left where u^ < c. 
Of course, on the diagonal* = c, corresponding to Taylor’s 
hypothesis. Off the diagonal the equation for vertical motion 
has the form 
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NOTE: Associated with each element in the array there are two numbers. The top number ., 

is the total time interval over which the particular case was integrated, and 
the bottom number is the linear period due to the forcing function. Both time 
periods have units of seconds. 


Figure 0-1. Total Time Intervals and Linear Periods 
for 16 Original Runs of Phase I 
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NOTE: Associated with each element in the array there are three numbers. The upper 


number is the theoretical (linear) period, the middle number is the observed . 
period, and the lower number is the phase lag. Positive values of the phase 
lag indicate the balloon lags of the wind. All values are in secpnds. 


Figure 0-2. Periods of Oscillation and Phase Lag 
for 16 Original Runs of Phase I 



(0-5) 


dv 3 9 

. + 2a v 3 + ITx^ = C cos ( k jXj - ait + e 3 ) + G 

where the left-hand-side of the equation is the usual equation for 
a damped oscillator. The first term on the right-hand-side is the 
forcing function due to the pressure forces where C'is assumed to 
be a constant, and G is a function which is nonlinear in‘ dx 3 /dt. 

If the function G is neglected, the balloon would be expected to 
oscillate with the period of the forcing function given by Equation 
(0-3). This expectation is borne out by the test case results. 

The largest observed deviation from the linear period occurs along 
the fourth column in Figure 0-2 where the wave number is a maximum 
and not necessarily where |u - c | approaches a minimum. 

Along the diagonal of the experiment matrix, the cosine 
of .Equation (0-5) becomes a constant. For these cases the oscill- 
ations should be expected to damp out as was observed in the A n 
case. In these cases, the balloon approximately matched the 
flow after the initial transient motion was suppressed and the 
forces balanced to hold it in this position. Thus for case 

, the balloon was observed to reach equilibrium about 120 m 
below its equilibrium level and, was observed to continue to 
have a lateral motion. This lateral mot ion . persists for longer 
time intervals, causing an ever-increasing lateral displacement 
of the balloon. Therefore, one would observe no vertical motion 
but a mean lateral motion. ' 

The vertical displacement typically displayed two modes 
of oscillation. The short -period oscillation corresponded to the 
forced oscillation. The long-period oscillation resulted from the 
initial conditions and was found to be damped as expected. 

Figure 0-3 shows case A 31 where a ^ery long total time 
interval for numerical integration was employed. This case clearly 
shows the damping in the transient mode of oscillation. 
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APPENDIX P 


Analysis of Numerical Results Phase I - Cruciform Runs 
P.l Cruciform C 2 2 

Twenty-four runs were carried out in the Cruciform C^ 2 
arrangement shown in Figure P-1. The values of k and w for each 
run are given in Figure P-1. For each run, a set of six position 
and velocity histories are generated. These plots are not included 
in the report , because of the considerable bulk which they represent 
Such plots have been collected and bound in a separate document 
tl] for reference pusposes. 

Analysis of the numerical results of the 24 runs was 
primarily concerned with comparing the first-order theoretical 
time period with the observed value. The results of this compari- 
son are presented in Figure P-2. In general, as the conditions 
corresponding to A^ are approached along any one of the four 
branches of the cruciform the first-order theoretical period 
and the observed period agree less and less, with the observed 
value increasing more rapidly than the first-order theoretical 
value. This - indicates that first-order theory, which is linear, 
is not sufficient to predict the balloon motion under conditions 
where Taylor’s hypothesis is (nearly) satisfied. 

P.2 Cruciform 

Twelve runs were carried in the Cruciform Cgg arrangement 
as shown in Figure P-3 . The values of k and <*> for these runs are 
presented in Figure P-3. As with the runs, a set of six 

position and velocity histories was generated for each run. These 
plots have been collected and bound in separate documents 1 for 
reference purposes. 

As in the- case of the C^ 2 runs, analysis of the numerical 
results of the twelve runs was primarily concerned with comparing 
the first-order theoretical time period with the observed values. 
Figure P-4 provides a summary of this comparison. Inspection 
of this figure reveals the same trend as observed in Cruciform C 22 • 
Again, it would appear that first-order theory is not adequate to 
predict the balloon motion under conditions where Taylor’s hypothe- 
sis is (nearly) satisfied. 
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Figure P~2. Dimensionless Time Periods for Runs in the Cruciform C, 
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